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LXVI. On the Inhomogeneous Deformation of a Plastic Lamina in a 
: Compression Test. 


By R. Hit, M.A.Ph.D., 
Metal Flow Research Laboratory, Sheffield *, 


[Received April 14, 1950.] 


SUMMARY. 


A theory is formulated for the calculation of stresses and displacements 
in a thin sheet of arbitrary shape which is being plastically compressed 
between parallel plates. Account is taken of inhomogeneous distortion 
caused by friction or edge constraints. By a transformation the solution 
is made to depend on that for a certain elastic problem of plane strain. 
Experimental data in support of the theory is provided by compression 
tests on long rectangular strips with both fluid and solid lubricants. A new 
method is suggested for estimating coefficients of friction relevant. for 
cold-working processes. 


$1. InTRODUCTION. 


‘THE paper is concerned with the behaviour of a metal lamina compressed 
plastically between parallel rigid platenst. The term lamina is to be 
understood as referring to a uniform plane sheet of arbitrary shape, 
having a thickness small compared with other dimensions. This restric- 
tion is made for the sake of simplicity, since the variations of stress and 
strain through the thickness can then be disregarded provided that the 
friction on the platens is not too great. 

So long as the platens overlap the specimen all parts of the latter are 
compressed equally and, if there were no friction, the deformation would be 
uniform and the lateral expansion equal in all directions. The slightest 
trace of friction, however, causes pronounced inhomogeneous deformation 
if some dimension of the lamina is very large compared with the thickness. 
For example, in a long rectangular lamina the compressive stress in the 
plane needed to overcome the frictional resistance to longitudinal expansion 
is of order 24,Yb/h, where h is the thickness, 2b the length, Y the compres- 
sive yield stress, and p the coefficient of friction. Thus, even when the 
lubricant is so good that p is 0-05, the compressive stress in the plane is of 
order Y/2 at the centre of a specimen whose length is ten times the thick- 
ness. In a plastically deforming metal this implies that local expansion 
in the length direction is almost or totally inhibited (the third principal 


stress being small). In general, a lamina changes its shape appreciably 


* Communicated by the Author. 
+ Platen is the technical term for a compression plate in a press. 
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during compression unless it is initially circular (the friction being 
isotropic). Moreover, the distribution of pressure over the platens is 
non-uniform, so that the test cannot give a direct measure of the uniaxial 
stress-strain curve. It is necessary to know the relation between the 
applied mean pressure p and the current yield stress Y in compression. An 
estimate of p for a circular disk, due to Siebel (1923), is 


= 2a 


where a and hf are respectively the current radius and thickness. This. 
formula should be a good approximation for a lamina. when ya/h is not too: 
large, say less than $. No similar formula has been suggested for a lamina 
of arbitrary shape, and it is the purpose of the present paper to show how 
one may be obtained. The change of shape of a non-circular lamina can. 
be regarded as a sensitive indication of the presence of friction. It is. 
certainly a more delicate indication than the barrelling of a solid cylinder, 
which is not known invariably to accompany friction. Furthermore, by 
the help of the theory to be described, the coefficient of friction can be 
deduced from measurements of the dimensional changes during compres- 
sion. 

When, on the other hand, the lamina overlaps the platens, the material 
outside the platens is not compressed, although a part of it will usually be 
deformed plastically by the stresses in its plane. A constraint is thereby 
exerted on the material between the platens, tending to produce an. 
inhomogeneous distortion which may overshadow that caused by friction. 
In a test suggested by Orowan and described by Ford. (1948), advantage is. 
taken of this constraint to secure compression under approximate conditions 
of plane strain. A rectangular strip is indented by platens completely 
spanning its width. The material outside the platens is almost entirely 
rigid and, except near the edges, prevents all strain in the width direction 
when the compressed area is sufficiently long and narrow. The falling-off 
in pressure towards the ends of the platens is investigated here. The 
situation resembles that in the roll-gap in strip-rolling, and the results of 
the investigation should furnish a method for estimating the state of stress. 
near the edges of the strip. 


§2. FRICTIONLESS PLATENS witH EDGE CONSTRAINT. 


We suppose, to begin with, that the platens are perfectly smooth, the 
non-homogeneous distortion being due solely to constraints applied to the 
edge of the material between the platens. Although the constraints would. 
in general produce a varia ion in stress and strain through the thickness 
this is confined to the part of the lamina lying within a distance of aon 
h from the edges of the platens*. It can therefore be neglected since we 
are assuming the width to be large compared with the thickness. 


* See, for instance, the solution by Hill and Lee (1946) to th bl 
plane indentation of a block by perfectly rough en ee eet a 
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Let Cartesian axes of reference (x, y, 2) be taken, with z normal to the 
lamina. The components of the mean stress (averaged through the thick- 
ness) satisfy the equations of equilibrium : 
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Let S denote the speed of approach of the plates at any moment ; the 
corresponding rate of compression is S/h. The Lévy-Mises plastic 
stress-strain relations are then 


kh _o,—hMo,+o,) ) 
S dx o,—4(0,+0,)’ 
h OV _ oy—4(2,+¢5) 2 
S oy O72 o4-0,) | 
‘y h dU ov = ST, 
S\oy ax} o,-Ha,+0,)’ | 


where (U, V) are the (x, y) components of velocity averaged through the 
thickness. We have, finally, the criterion of plastic yielding : 

ease Atel) ee ee ee cares Gt 
where o, and o, are the principal stresses in the plane of the lamina. 

Since the lateral strain is non-homogeneous, different elements receive 
different amounts of work-hardening. However, when the variation in 
the lateral component of strain is not too great (as we shall assume), all 
elements are in virtually the same state of work-hardening at any moment. 
Let Y be the current value of the uniaxial yield stress corresponding 
to this state. We now take advantage of the fact that the stress in any 
element will normally not be far from uniaxial compression in the z direction. 
Then, the hydrostatic component of stress being without influence on 
yielding, the criterion (3) can be approximated by the linear expression 

os (a5 0, )—=— ry, (hepa e's | SCY) Pee nee (4) 
where m is a constant which is to be chosen to give the best overall approxi- 
mation, and we have used the fact that o,+0, =0,+0¢,. When the 
values of the principal stresses are plotted along rectangular axes, the 
equation (4) represents one face of a regular hexagonal prism whose axis 
is in the direction (1, 1, 1). When m=1, the prism encloses the circular 
cylinder representing von Mises’ criterion and touches it along the gener- 
ators corresponding to uniaxial tension or compression. By suitably 
choosing m, equation (4) can therefore be made to approximate von Mises’, 
or indeed any, criterion for states of stress which do not depart far from 
uniaxial compression*. If m is assigned the value 1/3/2, (4) agrees with 


* Similarly, the maximum shear-stress criterion (represented by a hexagonal 
prism inscribed in the Mises’ cylinder) can sometimes be used as a convenient 
linearizing approximation to any actual criterion for states of stress near a 
pure shear. 

Be We 
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von Mises’ criterion in a state of pure shear. Hence by choosing m between 
4/3/2 and 1, we ensure that (4) is never in error by more than about 8 per 
cent for any stress in the range |o,—o,|<2mY/3. If the shear stress in 
the plane exceeds mY/3, the error can only be kept below 8 per cent by 
replacing (4) by another face of the hexagonal prism. This would not be 
analytically convenient and we shall therefore take (4) as an approxima- 
tion for all stress states. However, since (4) is introduced mainly to 
simplify the stress-strain relations (2), the accurate criterion (3) can still be 
used for an improved estimate of ¢, when approximate values of o, and 
a, have been calculated. 
On inserting (4) and (2), we find that 


o dU eae 3(9,— Fy) 
Sor * 4mY ° 
hii “Cp Oy) 
Say 2 dmY 6) 
x Aap BENE be etey 
S\ag * az) ~ my * 


We now substitute 


mY (hU 5 mY (hV_ , 
art Sa ie eek 3 v= “K (= —ly), . * (6) 


where E denotes Young’s modulus, and the new variables (w, v) have the 
dimensions of length. We obtain 


These will be RHE: as the elastic equations for a material which has a 
Young’s modulus E and Poisson’s ratio $, and which is deformed in plane 
strain. The distribution of stress and strain in the plastic lamina can 
therefore be derived from that in a geometrically similar lamina of 
incompressible material deformed elastically in plane strain. The (a, y) 
stresses are indentical in both, while according to (6), the elastic displacement 
is proportional to the plastic velocity representing the non-homogeneous 
part of the strain (7. e. over and above uniform compression and expansion). 
The z component of stress needed to maintain plane elastic strain in an 
incompressible material is $(7,,++-¢,,) ; thus, from (4), the pressure exerted 
by the platens on the plastic lamina is equal to mY plus the pressure on the 
elastic lamina. 

Consider, for example, a long rectangular lamina which is being com- 
pressed between smooth platens completely spanning its width (tig. 1). 
Let the width of the lamina be 26 and the width of the platens 2a. Let thex 
and y axes be taken parallel to the sides of the lamina, with the origin at the 
geometric centre of the area under compression. If the material lying 
outside the platens is assumed to be rigid, the conditions along the 
boundary of the rectangular compressed area are 


U=constant=-+ASa/h, say ; Viz=0;; onz=-+a; 
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together with zero applied stress on the surfaces y=-b. The unknown 
constant A must lie between } (corresponding to free expansion) and 1 
(corresponding to complete prevention of lateral strain). The precise 
value is determined by the condition that the resultant force over a 
lateral section, acting in the x direction, is zero (the ends of the lamina 
being unconstrained). The boundary-conditions in the associated elastic 
problem for a rectangular lamina (with sides 2a and 2b) are 


Fig. 1. 


ae PLATENS 
feet] 


¥ 
“LAMINA 


Lamina compressed between platens which overlap 
across the width only. 


| oe RIGID BLOCK 
a7 x 
Ui _~- CEMENTED EDGE 


Associated elastic problem. 


together with zero applied stress along y= +6 and zero resultant force on 
the surfaces x —-+ a. These conditions are satisfied by superposing a 
uniform tensile stress 2mY/3 in the x direction on the stress distribution 
due to a compressive load 4mY6/3 applied through inextensible rigid euus 
cemented along x= -La (fig. 2). According to (7) the tensile stress produces 
displacements 


+mYa/2E, —mYy/2H, on v= +4. 
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The compressive load balances this tension (so satisfying the condition for 
zero resultant force), and produces along « =-+ a a zero displacement in 
the y direction (since the bars are cemented and inextensible) and a certain 
uniform displacement in the 2 direction (since the bars are rigid). The 
magnitude of the latter has to equal (1—A) mYa/E, and this determines the 
value of A. 

The problem of compressing a rectangular elastic block between massive 
cemented bars does not appear to have been solved*. A photo-elastic 
investigation is projected, though the achievement of plane strain presents 
difficulties. We may, however, safely make certain generalizations about 
the distribution of stress, without a detailed solution. In the first place, 
by Saint-Venant’s principle, the stress must be virtually uniform except 
within a distance of order 2a from the free surfaces. Secondly, the pressure 
transmitted through the bars must fall off at the edges, since the material 
there has greater freedom to expand and can therefore be compressed. by 
a given amount under a smaller pressure. Thirdly, since the bars are 
cemented, they exert a shearing stress and this, by Saint-Venant’s principle, 
must be concentrated mainly within a distance of order 2a from the free 
surfaces. Translating these conclusions in terms of the plastic lamina, we 
can expect a state of plane strain except in a region near the edge, where 
there is both a shearing stress and a ten_ion parallel to the edge. The 
pressure applied through the platens therefore diminishes towards the 
edge, the fall-off being confined to a region whose width is roughly equal to 
that of the platens. 


$3. Roueu PLatens WirHout EDGE ConstTRAINT. 


We consider, now, a lamina free from edge constraints and compressed 
by rough platens which overlap it all round its circumference. The 
coefficient of friction uw is assumed so small that terms containing second or 
higher powers of . can be neglected. For example, c, is a principal stress 
to first-order, and so the yield criterion is still (3) or (4). The theory is only 
a valid account of a first-order departure from pure uniform compression. 
This is the reason for excluding edge constraints. 

It is necessary to determine which element of the lamina remains 
momentarily at rest at each stage of a slow compression. In a fast 
compression, where inertial forces are large compared with frictional 
forces, it is clearly the element momentarily coinciding with the centroid, 
since there are no other lateral external forces. The situation is different, 
however, when the compression is so slow that inertial forces are 
negligible compared with frictional forces. For a first-order calculation 
of the distribution of the frictional shear stress over the faces of the lamina, 
the deformation may be taken to be uniform expansion. We must also 
allow for the possibility of a small rotation of the lamina in its plane. Then, 


* An approximate method, yielding limited information and applied only to 
a square block, has been described by H. J. Greenberg and R. Truell (1948) 
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‘to first order, the relative velocity of an element distant r from the 
momentarily fixed point is compounded of a radial velocity proportional 
to r (due to the uniform expansion) and a velocity which is perpendicular to 
this and also proportional to 7 (due to the rigid-body rotation). The 
resultant velocity is therefore proportional to r and inclined at a fixed angle 
4 to the radius. Hence the frictional stress, opposing the relative motion, 
is directed at the same constant angle @ to the radius. To first order, the 
magnitude of the frictional stresses is ~Y, and they therefore exert no 
resultant torque on the lamina if 


pwY sin || rdx dy =0, 


taken over the lamina. Since the integral is essentially positive, it follows 
that 6=0, or that there is no rotation of the lamina as a whole. The 
friction is therefore directed toward the po'nt at rest, and the condition for 
zero resultant frictional force is 


| | tx dy=o, jetties, *d'T 68) 


where t is the unit vector in the radial direction. This is the equation 
determining the position of the element momentarily at rest. If the 
Jamina has two axes of symmetry, the centroidal element always remains 
at rest. For a lamina of arbitrary shape, however, the point defined by 
(8) is not the centroid. 

Let (x, y) axes be taken with the origin coinciding with the point at rest. 
The frictional stress contributes to the equations (1) additional terms which 
are equivalent to a body stress with potential 2uYr/h. Such a body 
stress can be balanced by a variable hydrostatic tension whose value 
is equal to the potential. The remaining part of the actual stress then 
satisfies the equations (1) of equilibrium. Furthermore, the remaining 
part of the actual stress still satisfies the plastic stress-strain relations and 
yield criterion, since these are unaffected by the addition of a hydrostatic 
stress. The corresponding external stress on the edge of the lamina must 
be equal to a normal pressure of amount 2uYr/h at each point in order to 
produce, when combined with the hydrostatic tension, zero resultant stress 
on the edge. 

Thus, to obtain the distribution of stress in a lamina compressed by 
rough platens, we combine a varying hydrostatic tension 2uYr/h with the 
distribution of stress in a similar lamina compressed along its faces by 
smooth platens and along its edge by a varying normal pressure 2uYr/h. 
The deformation is identical in the two lamine. Since the edge pressure 
is proportional to the distance from the point at rest, it offers the least 
resistance to lateral expansion of the parts of the edge nearest to this 
point. In other words, as one would naturally expect, friction on the 
platens causes the lamina to expand most in the direction of its least 


dimension. 
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The elastic analogy described in §2 is, of course, still applicable. 
Let us, for example, apply it to determine the stress in a circular lamina of 
radius a. The corresponding edge pressure produces a uniform hydro- 
static pressure 2uYa/h in the plane of the elastic lamina and hence also in 
the plastic lamina between smooth platens. Combining this with the 
varying hydrostatic tension 2uYr/h, we obtain equal radial and cireum- 
ferential compressive stresses 2uY (a—r)/h in the plane of the lamina 
between rough platens. From (4), taking m=1, the pressure exerted by the 
platens is 


Y+2unY(a—r)/h 


and the mean pressure is Y (1+ 2ua/3h), in agreement with Siebel’s formula. 


§4, EXPERIMENTAL DETERMINATION OF THE COEFFICIENT OF FRICTION. 


The preceding theory provides a means of estimating the coefficient of 
friction in a compression test from observations of the change in shape of 
a rectangular lamina. 

The coefficient of friction, so determined, corresponds to the conditions 
in a compression test, viz: where one of the sliding bodies is distorted 
elastically in bulk and the other has yielded plastically. Such coefficients 
of friction are relevant for all processes for cold-working metals . 

The method is obviously more sensitive when the change of shape is 
large, but then the first-order approximation of the theory is less accurate. 
However, when the lamina is long and narrow, a more exact calculation 
can be made. If the length 26 is greater than about 10 times the width 
2a, and if a/h is of order unity, it is probably legitimate to neglect the 
variation in stress across the width in comparison with the variations 
along the length. The only stress components are then the longitudinal 
compressive stress q, the pressure p, and the friction wp. To the same order 
_ of approximation we may disregard the change of shape when setting up. 
the equation of equilibrium. This then reduces to 

d 
lege ORD A Sawn Eg A eM) 
where y is measured along the lamina from its shorter side 2a. For any 
yield criterion near that of von Mises, it is sufficiently accurate to use 
equation (4) : 
oS tg tr Vous: Wy Gal aya, eee ade 
Eliminating ¢ from (9), and integrating, we obtain 


prem exp (kyl) aan. x dot oe te LT 


after using the boundary-condition p=mY on y=0. 


The longitudinal strain-increment vanishes at a section distant y’ from 
the shorter side, such that g=4p or 


py'Jh= In (4/3) 220-288, . -. « . . . (12) 


} 
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Hence, if ub/h <0-288, all parts of the specimen continue to elongate. It 
follows from equations (2), (10) and (11), that the increase in 6 during a 
further small decrease dh in thickness is 
-b 
db= wi ; (2—3 eM*/*) da—dh E mat i (enn | , 
This may be written to the present order of approximation as 
db (1 ) dh 
| On 2h J 2h 
Integration gives 
_# 


ci ee fo) HO <o.ogs- | 
aalloeeen g) SX Sata myetien) (13) 


2h 
where b, and hy are initial values. The first factor in this equation becomes 
increasingly less than unity as the compression proceeds, and represents 
the effect of friction in hindering changes of length. The corresponding 
formula for the mean pressure is, with m~1, 


Box (14 ) Y; a LAU e Mn | Pe ele) 
On the other hand, when pb/h>0-288, a central region of the lamina, of 
length 2(b—y’), momentarily undergoes no extension. Since. there is no 
relative longitudinal slip in this region, p has the uniform value 4mY/3 and 
q is equal to 2mY/3 (m should be taken nearer to +/3/2 as y’/b becomes 
progressively less). The mean pressure is easily found to be 
pa (j —0-001 5 ) TN ue 0:2 88a 24 age sal ae (Lo) 
The increase in length during a small decrease dh in thickness is given by 
db h \ dh 
if ae 
y' ( A 2h 
db 0-076. 
aie was 

In the present experiments the nominal specimen dimensions were 
1-0 0-10-05 in., so that b,/hj>=10. Fig. 3 shows the corresponding 
theoretical relation between the fractional change in length (b/b)—1) and 
the fractional reduction in height (1—h/h,) for various values of. At the 
critical value of b/h, for which the centre of the lamina undergoes no ex- 
tension, the curve (13) is replaced by a line of slope 0-076 ho/ubo, according 
to (16). When pis greater than 0-0288 the relation is linear from the start 
when 6,/hj>=10. 

Also shown in fig. 3 are experimental points for two lubricants : 
powdered calcium oleate and vacuum solvac. Copper specimens were used, 
the respective dimensions being 

1-000 x 0-100 X 0-050 in. and 1-002 x 0-103 x 0-049 in. 


or 
Ho 0-288, ig, ape, 24 Us7(16) 
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‘The copper was rendered isotropic by suitable rolling and heat 
treatment ; this was confirmed by compression tests on circular disks, 
which femained effectively circular The lubrication was renewed when- 
ever the test was interrupted to take measurements of length and thickness. 
The experimental points follow the trend of the calculated curves, thus 
affording indirect evidence that the theory is sound. The deduced 


Fig. 3. 


O15 


0.10 


0.0$ 


Comparison of calculated and observed values of the extension of lamine of 
dimensions 1-0 x 0-1 x 0-05 ins. 


coefficients of friction are 0-015 and 0-045, respectively. The latter is 
comparable with values inferred from measurements of roll-force when 
rolling short lengths of copper strip with vacuum solvacasa flood lubricant *. 


a ee ee eee 
* Private communication from Mr. R. B. Sims, Head of the Rolling Section, 
‘British Iron and Steel Research Association. 
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‘The low coefficient of friction obtained with calcium oleate indicates that 
‘true boundary-lubrication was not operative, and this was confirmed by 
the appearance of the copper after compression. The surface was still 
mainly in the ‘ orange-peel ’ condition, in marked contrast to the smoothed 
and burnished surface of the specimen lubricated with vacuum solvac. 

For a direct check of the theory it is necessary to use a lubricant whose 
frictional properties are known independently. 

A sufficiently thin foil of some softer metal acts as an efficient lubricating 
agent, the frictional stress being equal to the limiting yield stress ¢ in 
shear of the foil at large strains. The foil must not be too thick, or it will 
extrude in the width direction before (or while) the specimen itself is 
compressed, and so exert a negative frictional drag in that direction. A 
rough criterion, ensuring non-extrusion, is that Y/s should be less than w/t, 
where ¢ is the foil thickness and w is the specimen width. This follows 
from the condition that the hydrostatic component of the stress at the 
middle of the foil, necessary to cause extrusion in the width direction, 
should be greater than Y. For lead foil on copper, Y/s is less than about 
25/0-5 (units of tn./in?.) at all strains, and so ¢/w should be less than 1/50. 
In the present experiment ¢ was not more than 0-002 in., while w was 
initially 0-1 in. The foil was removed and replaced by a new piece each 
time measurements were taken. The results are shown in fig. 3 and corre- 
spond to an apparent coefficient of friction of about 0-017 (the changes 
needed in the theory when the frictional stress is constant are insigni- 
ficant here). By comparison s/Y is about 0-4/20 or 0-02 for large 
strains. The agreement is excellent in view of the various uncertainties. 
Similar confirmation was obtained with pure lead foil on aluminium 
where the apparent coefficient of friction was found to be 0-042 
(s/Y~0-4/10 = 0-04), and with indium foil on a lead alloy where the co- 
efficient of friction was 0-033 (s/Y~0-1/3 = 0-033). The limiting values of Y 
_at large strains (greater than 20 per cent compression) were obtained from 
measurements of the load and thickness during the test (the current area 
being calculated as initial volume divided by current thickness). The 
limiting shear stress of the indium foil was inferred from hardness tests on 
the bulk metal. An 4 in. diameter ball under a 4 lb. load of 20 seconds 
duration gave an impression diameter of 0-066 in. (average of three tests). 
‘This corresponds to a Meyer hardness of 0-52 tn./in®. and hence to a 
mean uniaxial yield stress of about 0-17 tn./in®. over the range of strains 
produced around a ball indentation. 

For indium foil on copper and aluminium, or lead foil on steel, extrusion 
must have occurred since the change in length was much smaller than would 
be expected if the foil had behaved as a true lubricant ; extrusion would, 
indeed, be expected by the criterion mentioned previously. In agreement 
with this hypothesis, a solid cylinder of copper showed barrelling when 
lubricated with lead, but a concavity when lubricated with indium. 

"A further check on the theory is provided by a comparison of the 
observed and calculated values of the length of that part of the specimen 
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which is deformed throughout in plane strain. Theoretically this is; 
from (12), 
0-288 hg 


For the copper lubricated with vacuum solvac, where b)=0-501, hp =0-049,. 
u=0-045, this unextended central region is calculated to be of length 
0:37 in. This is indicated on fig. 4, which shows the measured profile of the 
lamina after a reduction in thickness of 60 per cent. The agreement is. 
excellent. The final maximum width is 0-248 in., while the width cal- 
culated on the assumption of plane strain is 0-26 in. When calcium oleate 
was the lubricant the profile was completely rounded, again in agreement 
with theory since p.b)/hy=0°15 <0-288. 

It should be carefully noted that the theory becomes less accurate as the 
lamina is progressively compressed, since the length/width ratio decreases. 
from the initial value of 10 to about 5 after a 50 per cent reduction. Tests. 
on laminz with an initial length/width ratio of 4 showed that the theory 
significantly underestimated the increase in length ; this was to be expected. 


Fig. 4. 
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from the neglect of the component of the frictional forces restraining the: 
expansion in the width direction. This effect can be recognized in fig. 3, 
where the upward trend of the experimental results at large reductions. 
is more pronounced than the theory predicts. For this reason the 
coefficients of friction were deduced from the experimental results below 
about 40 per cent, where the length/width ratio is greater than 7. 
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SUMMARY. 


The theory of plane plastic strain due to Hencky and Geiringer is 
applied to calculating the minimum dimensions that a specimen should 
have for a valid hardness test. Critical widths and thicknesses are 
determined for rectangular specimens indented by flat dies and wedges 
of any angle. The manner in which the specimen deforms during the 
indentation is considered in detail. 


$1. INTRODUCTION. 


THE definition of indentation hardness adopted in the present paper is 
the ratio of the applied load to the projected area of the impression. 
It is equal to the work done per additional unit volume of indentation, 
or, when the frictional resistance is negligible, to the mean normal pressure 
between indenter and medium. It is well known that the hardness is 
not observably dependent on the size of the specimen provided this is 
greater than a certain minimum (other things being equal). The minimum 
size depends mainly on the properties of the material under test and on 
the shape of the indenter. If the specimen is too thin or too narrow, 
the plastic zone penetrates either to the rear or lateral surfaces ; this 
is accompanied by a change in the apparent hardness (usually a reduction). 

The minimum width and thickness of a rectangular specimen are 
calculated here for indentation by a wedge. At present, this is the only 
indenter for which a reasonably accurate calculation of the surface 
distortion and stress distribution is feasible. The problem is treated as 
two-dimensional, the strain being zero in the direction made by the 
intersecting faces of the wedge. The material is supposed to be rigid 
up to the yield point, after which the maximum shear stress retains a 
constant value in plane strain. 


§2. CRITERION FOR THE CRITICAL THICKNESS. 
The slip-line field ABCDE in the region of plastically-deforming 
material around a wedge indenting the plane surface of a semi-infinite 
mass is shown in fig. 1 (Hill, Lee, and Tupper 1947). The wedge angle is 


* Communicated by the Author. par 

For an account of the theory of plane plastic strain, see the writer’s 

The Mathematical Theory of Plasticity, Chap. VI. (Oxford: Clarendon Press, 
1950). 
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denoted by 20, the depth of penetration by c, and the length of the line’ 
of contact by h. The displaced surface is straight, and, when the wedge 

is smooth (as assumed here), the coronet is of length h. The inclination. 
of the displaced surface to the horizontal is denoted by 6—; such a 

raised coronet is characteristic of a cold-worked metal (Heyer 1937). 

Fig. 2 shows the relations between h cos 6/c, 4 and @. When the specimen. 
is semi-infinite and homogeneous the configuration remains geometrically 

similar and the hardness is independent of load. This may be realized 

in practice, for example in a Vickers pyramid hardness test, provided the: 
impression is large compared with the grain-size of the specimen and the: 
radius of curvature of the tip of the indenter (Field Foster 1948). 


Big? 1. 


sah 


ABCDE does not constitute the whole plastic region. Beyond this 
field lies a further plastic zone where the material is stressed to the yield- 
point but remains rigid, being prevented from deforming by the 
constraint of the non-plastic matrix (Hill 1949). The state of stress and 
strain in the material undergoing deformation is independent of the 
distribution of stress in the rigid material. Consequently, there is no 
need to calculate the latter whenever it is eviaent that the rigid material 
is capable of supporting the stresses acting around its perimeter. This is 
necessarily so when a specimen is sufficiently large compared with the: 
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localized area of deformation. It follows that conditions in the plastically- 
deforming zone are independent of the size of the specimen, provided 
only that it is sufficiently great. In particular, the configuration around 
a wedge remains geometrically similar (though the shape of the plastic 
boundary does not). 

These conclusions are strictly valid only for a material that remains 
rigid whenever the maximum shear stress is less than k. In an actual 
solid, where Young’s modulus E is finite, the deformation in the 
plastic region is dependent on the elastic distortion elsewhere. For: 
example, the volume of the coronet is governed by the extent to which 
the elastic resilience of the whole specimen can accommodate the 
displaced material. Such departures from the behaviour of a plastic-rigid. 
body are presumably only of order k/E, and the hardness should be 
independent of specimen size over a certain minimum, within normal. 
experimental error (as is observed). 


Fig. 2. 


In a specimen whose thickness is small compared with its width, the- 
plastic boundary moves relatively further from the indenter and nearer 
to the rear surface as indentation proceeds. When the specimen. rests. 
on a plane rigid foundation, it is a natural supposition (and one which 
is suggested by observations of the plastic distortion of the rear surface) 
that the relative expansion of the plastic zone is mainly concentrated 
directly below the indenter. At some stage of the indentation the element 
situated at the midpoint O of the rear surface is stressed to the yield limit. 
At a little later stage the plastic region contains all the material within a 
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certain area which lies below the indenter and intersects the rear surface 
over a finite length (the conjectured plastic boundary at this stage is 
indicated by broken curves in fig. 1). It may be proved (Hill 1949) 
that deformation is confined to the area ABCDE up to the moment when 
the plastic region has spread so far that it just includes the two slip-lines 
OS from O. Thereupon the constraint fails, the plastic core deforms, and 
the rigid ends are thrust apart, sliding over the foundation. The 
configuration round a wedge no longer remains geometrically similar 
and subsequent hardness measurements are spurious. 

We now examine the criterion determining the critical depth of 
penetration. The slip-line field bounded by OS and AS is uniquely 
determined by the known shape of AS; the particular position of S on 
ADKEC depends, of course, on the ratio of AO to AB. The field having 
been calculated, Hencky’s theorem determines the normal stresses acting 
across AO, and hence the total lateral thrust exerted on each half of the 
specimen by the stresses along OAB. The thrust vanishes when the 
‘foundation is smooth and the ends of the specimen are unconstrained. 
Now the field OSAS could always be formally constructed at any stage 
of the indentation, but only for one particular depth of penetration 
would the corresponding thrust across OAB vanish. Before this depth 
is attained, the plastic region can therefore not have spread sufficiently 
far to enclose the entire length of the slip-lines OS. Hence, as already 
stated, the material below the indenter must have remained rigid. The 
criterion for the critical penetration is therefore that there should be no 


resultant lateral thrust over OAB corresponding to a formal construction 
of the field OSAS. 


$3. CALCULATION OF THE CRITICAL THICKNESS. / 


Slip-line fields with a 5° mesh were constructed for smooth wedges of ° 
60° and 120° total angle by the numerical method due to Hill and Lee 
(1946). The field OSAS consists of a region ADF where the slip-lines are 
radii through A and circular ares ; a region DFHS defined by the circular 
ares DF and DS of equal radius; a uniformly stressed region AFG ; 
a region FGHJ where one family of slip-lines are straight ; and a region 
GJO defined by the slip-line GJ and the axis of symmetry. 

If ¢ is the minimum specimen thickness, it is found that ¢/2h is equal 
to 1-04 for 6=30° and to 2-17 for 6=60°. The corresponding distributions 
of normal stress across OA are shown in fig. 3; it will be observed that 
this stress is tensile immediately under a narrow wedge but compressive 
under a wide wedge. The stress in the rear surface is always tensile. 
When ¢@ tends to zero, so that the indenter becomes a knife-edge, the 
pressure 2k along AB is balanced by a uniform tension zk across OA. 
t/2h is therefore equal to (7-+2)/27 or 0-82 approximately. At the other 
extreme, when @ tends to 90°, the slip-line field becomes similar to that 
for a flat die, for which the critical t/2h is 4-37 (Hill 1947). The relation 
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‘between the critical ratio t/2h and the semi-angle @ is given in fig. 4*, 
The critical ratio of load to thickness is also shown ; since this has a 
maximum value when 6~60°, the hardness of a given specimen can be 


measured over a greater range of loads with a 120° wedge than with any 
other. 
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The effect of friction on tue foundation is to delay the moment when the 
eritical penetration is reached. Calculation shows, however, that the 
effect is so small as to have little practical importance for coefficients of 
friction less than 0-1. Lateral compressive forces applied to the sides of 
the specimen would be more effective (so long as they do not cause general 
yielding). 

+ Fi mparison, the critical ratio of specimen thickness to impression 
pene 3-14 when a semi-cylindrical cavity is expanded in the surface 
(Hill 1948). 
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The speed U with which the ends of the specimen slide sideways is- 
determined by the following considerations. If the speed V of penetration 
is given, and U is assigned arbitrarily, the velocity at any point in the 
plastic region can be calculated from Geiringer’s equations. In particular, 
the velocity in AFG, ABD, and BEC is uniform, so that the displaced 
surface BC remains straight during the next infinitesimal increment of 
penetration. (As the penetration continues, the part of the coronet 
above the element now at C remains straight, but the additional lower 
part is curved.) Corresponding to a slightly altered position of the 
surface, it must be possible to find a new plastic state of stress throughout 
the material which was plastic prior to the increment of penetration. 
Otherwise plastic material unloads during the incremental penetration 
and this is inconsistent with the assumption of continuing plastic flow 
implicit in the application of Geiringer’s equations*. This is the condition 
which determines the true ratio U/V. A calculation of U/V has not been 
made by the writer. 


§4. CALCULATION OF THE CRITICAL WIDTH. 


In a specimen whose width is small compared with its thickness, the 
plastic region spreads out to the sides as indentation continues. At a 
certain penetration the plastic region (whose conjectured boundary is 
indicated by broken curves in fig. 5) just contains a complete slip-line 
HJKLM, extending from the axis of symmetry to the lateral surface 
(which it must meet at an angle of 45°). Up to this moment the plastic 
material below the indenter has been rigidly constrained. After this. 
moment the rigid corners of the specimen are assumed to be displaced 
sideways by sliding over the slip-line KLM; for this to be possible, 
KLM must be of a suitable shape. Now the development of a plastic 
region having a boundary cut more than once by a slip-line family is not 
uniquely determined by the state of stress alone, but proceeds in such 
a way that certain conditions for the existence of displacements are 
satisfied (Hill 1950). It is reasonable to suppose, therefore, that the 
plastic region spreads out to the sides of the specimen in such a manner 
that the slip-line KLM is of the correct shape to permit a sideways 
displacement of the rigid corners. Since a rigid-body motion in a plane 
can only be a rotation about some instantaneous centre (which is at 
infinity for a pure translation), two rigid bodies can only slide over each 
other if the line of contact is a circular are. Thus KL must be circular. 
We must also include a straight segment LM such that LMN is a region 
with straight slip-lines meeting the surface at 45°. The stress in LMN 
is a uniform compression 2k parallel to the surface, while the ensuing 
deformation caused by the rotation of the corner is a non-uniform shear 
parallel to LM. 


* These remarks are of general significance for removing an apparent lack 
of uniqueness in many problems where a plastic-rigid body is assumed. 
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Since the nett angle turned:through in going from M to H along the 
slip-line is zero, the hydrostatic components of the stress at H and M 
are equal; their value is —k, since the stress in LMN is a compression 
of amount 2k. Hence H is the point on the axis at which the hydrostatic 
stress-component is —k in the unique continuation of the field ABCDE. 
| By Hencky’s theorem, H lies on the slip-line through B which is inclined 
at an angle {m—}(0—y#) to BE. For @ less than about 50°, H coincides 
with A, and the slip-line HJKLM makes an angle }m—4(@—) with AD. 


H having been fixed in this way, there are three unknown quantities 
to be calculated. These are the length b of JK (all slip-lines of this 
family of EJKP are of equal length), the radius ji of KL, and the length 
a of LM. When these are known, the critical width of the specimen 
follows at once from simple geometry. To determine the three unknowns 
there are the three conditions for the equilibrium of the rigid corner, 
viz. : two conditions for the vanishing of the resultant force components, 
and one for the vanishing of the resultant couple. The resultant force 
and couple are calculated from the distribution of stress over the perimeter 
CPKLN, or, more conveniently, from the distribution over BAHJKLM 


a2 
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(the equilibrium of the region covered by the slip-line field is, of course, 
automatically ensured by the equations of Hencky used in constructing 
the field). Me 

When the indenter is a knife-edge an explicit analytic solution is 
possible (fig. 6). The perimeter of the rigid corner has an axis of 
symmetry bisecting the angle at the corner. It may be shown without 
difficulty that 


4 4 
afh=/2— —, Note ar a eer 


The critical width is therefore 2h(7+4)/7 or about 4-54 times the depth 

of penetration. When 6=90° (flat die) calculation gave the values 
ajh=0-95 , b/h=0-54 , rih=5-67 . 

The critical width is 8-6 times the width of the die. M is 3-67 die-widths 

below the corner. When 6~50° the point H coincides with A (@+4~90° ; 

see fig. 2). The corresponding values of a, 6 and 7, were calculated to be 
aj/h=0-205 , blh=0-51 , t= 1:99. 

The critical width is about 4:5 times the width of the impression. M is 

2:14h below the corner. 


$5. CONCLUSIONS. 
The minimum dimensions for the two-dimensional model are somewhat 
greater than those for actual indenters. For example, the minimum 
thickness recommended with a 136° pyramid is about twice the diagonal 
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of the impression (Hankins and Aldous 1934). The difference is under- 
standable since there is greater freedom of flow around a_three- 
dimensional indentation. Moreover, the plastic zone penetrates less far, 
owing to the diffusion of stress concentrations through a relatively 
larger volume. The present calculations have, therefore, a limited 
practical value. The value of the investigation lies rather in its 
contribution to a general understanding of a hardness test, especially 
in showing why the hardness value is independent of specimen size, 
down to dimensions which are comparatively small in relation to the 
indentation. The investigation is also of fundamental interest in that ~ 
it formulates methods for solving plane plastic problems of a kind different 
from any previously attempted. 
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SUMMARY. 


This paper is based on the theory of the metal semiconductor contact 
of Schottky and of Davidov, and is meant (a) to make this theory 
applicable in practice to all existing semiconductor types, and (5) to 
introduce quantitatively the influence of the electric field at the interface 
on the threshold energy. It appeared to be useful to introduce an in- 
dependent variable y, proportional to the reciprocal of the concentration 
of the conduction electrons. A differential equation of the Abel type is 
derived for the electric field in the semiconductor as a function of y. 
Taking into account the diminution of the threshold energy by the 
electric field at the interface, a relation is derived between the variable y 
and the electric field at the interface. The differential equation is given 
for the different types of semiconductors. When the differential equation 
is solved, the current-voltage characteristic can be found by numerical 
integration. Some examples are given. In the case of back currents, an 
asymptotic solution of the differential equation exists, which leads to a 
parametric representation of the current-voltage characteristic in the 
back direction. The shape of the back current characteristic is discussed. 
The differences between rectifiers of the selenium type and the crystal 
detectors are elucidated. Some remarks on the physical background of 
the electrical pretreatment of selenium rectifiers are made. 


$1. INTRODUCTION. 


Tue first theory of the blocking-layer rectifier, not based on ad hoc 
hypotheses, was given almost simultaneously by Davidov (1939) and by - 
Schottky (1939). These theories differ in their mathematical form, but 
are based on the same physical conceptions (Gisolf 1942). They account 
satisfactorily for the current dependence of the resistance for forward 
and for small back currents. (The “ forward ” direction is the direction 
for which the resistance is low.) For larger back currents considerable 
deviations from the theory occur. The resistance first passes through a 
maximum and tends to vanish for large back currents. Schottky (1941) 
did mention the presumable cause of this deviation, namely that in 
consequence of the strong electric field in the semiconductor at the 
boundary of the metal electrode, the threshold energy for the transition 


* Communicated by the Author. 
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of charge carriers from the metal to the semiconductor is diminished and 
the number of these particles in the semiconductor is increased. In the 
present paper a treatment of the above problem is given, which enables 
one to derive formule for the resistance of back currents with due 
allowance for the influence of the field on the threshold energy. The way 
by which this result is achieved may be outlined as follows. The 
concentration of the free electrons in the semiconductor is a univalent 
function of the distance « of the point under consideration from the 
surface of the electrode. Instead of a an independent variable y is 
introduced, which, except for a constant factor, is equal to the reciprocal 
of the concentration of the free electrons. From the fundamental formula 
of the electric field a differential equation is derived which gives the 
electric field strength in the semiconductor as a function of y for each 
value of the current density. The equation can be solved by numerical 
methods for the different types of semiconductors. A boundary condition 
is given by the field strength at a long distance from the electrode where 
Y=Yx.. Schottky’s formula for the diminution of the threshold energy 
at the electrode as a function of the field strength together with the 
relation between the field strength and the variable y, leads to the value y, 
of the variable y at the electrode. The resistance and the capacity of 
the metal semiconductor contact are found as integrals of completely 
known functions of y between the limits yy and y,,. 


§2. THEORETICAL PICTURE OF A SEMICONDUCTOR. 


In the following pages ordinary electrons with a charge —e, as well 
as “electron holes” with a charge +e will be spoken of as ‘ electrons ”’, 
It is possible to classify the rich variety of semiconductors on account 
of the electron energy scheme (Gisolf 1947). The following types may be 
cited :— 

A. The concentration of the free or conduction electrons n, is large 
compared to the concentration of the electron traps N,. The free 
electrons originate from impurity atoms nearly all of which are ionized 
by the crystal lattice. These impurities may be a stoichiometric excess of 
metal or metalloid atoms (for electrons and electron holes respectively), 
similar to those constituting the crystal lattice. The concentration , 
is then practically independent of the temperature. The conductivity 
slightly increases with decreasing temperature. | See 

B. Although only a fraction of the impurity atoms 1s ionized by the 
crystal lattice, the concentration of the free electrons n, is still large 
compared with the concentration of the electron traps N, The degree 
of ionization increases with increasing temperature. The concentration 
of the free electrons, and in consequence the conductivity, increases 
exponentially with increasing temperature. A definite ionization energy 
can be attributed to the impurity atoms. This case Bis the only one 
that is treated in. textbooks though it rarely occurs in practice. 


756° J. H. Gisolf on the 


C. In this case the concentration of the conduction electrons 7, is small 
compared to the concentration of the electron traps N;. One can assume, 
without loss of general applicability, that the impurity atoms are all. 
ionized. The trapping energies are supposed to be practically homo- 
geneously distributed over a range of energies JE, which is large compared. 
to the Boltzmann energy kT. In this case the concentration of the free 
electrons, and consequently the conductivity, is again an exponential 
function of the temperature. The pseudo “ionization energy ’’ deduced 
from this exponential dependence, however, now strongly depends on 
the concentration of the impurity atoms N; and in such a way that it 
becomes a universal linear function of the logarithm of the conductivity 
(the Meyer rule). 

The electrons of the semiconductor, which do not belong to the atoms 
of the normal lattice, are either bound to impurity atoms or are ‘‘ free ”. 
In any case they are supposed to have a definite energy. The energy of 
a free or conduction electron will be denoted by E, and the energy of an 
electron bound to an impurity atom by E;. The concentration of the 
conduction electrons 7, is given by the formula 


NN J{l-bexp(E,—g)/kT} > 2. ee 
with 
Inm*ET\ 3/2 
N= (=) (2) 


In these formule m* is the “ effective” electron mass. m* may be 
several times the normal electron mass but the exact value cannot be 
fixed. T is the absolute temperature; / is Boltzmann’s constant and 
his Planck’s constant. N,210!%cm.-3; g is the free energy per electron. 
In a stationary state the value of g is the same at every point. 

The concentration n, of the electrons bound to the impurity atoms with 
the energy E,, is given by the analogous formula 


m=N,/{l--exp (Be-g)/kT}. .. . . . 2. (8) 


In this formula N, is the concentration of the impurity atoms to which 
electrons can be bound with the energy E;. With the abbreviation 


= OX (HG) Ee ete. | wale) ack ta Eee) 

the formule (1) and (3) can be written 
n=N,/I+y)ZN_/y, (lopiie-< 1018) 2. ao ante} 
n =N,/{1+yexp(—E,+E,)/kKT}SN /yexp(—E,+E,)/kT, (forn;<N,) 
(6) 


When the contact surface of the metal and the semiconductor is flat, 
y only depends on the distance from the contact surface. 

The space charge density p can readily be expressed in terms of the 
parameter y for each of the semiconductor types mentioned above. It 
is assumed that the impurity atoms in the crystal are electrically neutral 
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so long as they have not lost electrons by thermal dissociation, For the: 
“traps ”’, on the other hand, it is assumed that they are electrically 
neutral so long as they have not captured an electron. The following 
formule for the space charge density p are found : 


For the case A : p=e(N,;—N,]/y). SO ee ek ee ek 
For the case B : o=eN,Y TYA—N,IN) NAN exp 4E/LT} 
{y+y(1-+exp JE/kT)+exp JE/kT} 
we Ny’ —Nolytexp 4E/kT)} 


“Wetexp ERD) ME MIT yy eH 
AE being the dissociation energy of the impurity atoms. 
For thecaseC: p=e{N,—N,1—kT/4Eln y)—N./y}, . . . . (9) 


AE now being the range of the trapping energies of the electron traps. 

The types A, B and C represent extreme cases. All practical cases 
are more or less intermediate. The bulk of the semiconductors belongs 
to type C. The discussion of this type is important for the theory of 
the crystal detectors. Cuprous oxide belongs most probably to type B. 
For selenium the situation is not clear. From measurements on single: 
crystals one might conclude that selenium belongs to type B, with a 
rather small value of JE. For polycrystalline selenium AE strongly 
differs for variously treated samples, but it is always small. Since the 
value of y is large in the blocking layer, the formule (7) and (8) show that, 
when dealing with blocking layer problems, one may reckon selenium to. 
class A, as has already been pointed out by Schottky (1939). 


§3. DERIVATION OF THE DETECTOR EQUATION. 


In the semiconductor the following two fundamental equations hold 
universally. 
(a) The Poisson equation 


dWigpen obec Ay =<. .)- ap on eee CLO) 


(b) The “ continuity equation ” for the current 
j= eur tkludn jdt, . ... - . + ,- (1) 


where F is the electric field strength, « is the dielectric constant of the 
semiconductor and ¢)=8-859 . 10-14, A sec/V cm. is the permittivity of 
free space, 7 is the current density, ¢ is the value of the electronic charge, 
u is the electron mobility in the semiconductor. 

The x-axis is perpendicular to the interface of the metal and the semi- 
conductor, the positive direction pointing to the semiconductor. By 
means of formula (5), , in equation (11) can be expressed in terms of 
the parameter y. A simple transformation results in 


dy/dx—=(F—iyjeuN,jeyfkT, - - . . » ~ (12) 
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y being a function of a only, equation (10) can be transformed in the 
following manner : 
dB /dy .dyjdx—ply)jegp. - - - > + + = (13 ) 


The function p(y) is given by one of the formule (7), (8) or (9) according 
to the type of the semiconductor under consideration. Substituting (12) 
n (13) one obtains 

dF /dy . (F—iy/euN,ey/kT=ply)lecp. - - - + (14) 


The solution of this non-linear differential equation gives F as a function 
-of the parameter y. 

Equation (14) can be shortened, by introducing an auxiliary quantity 
f(y), defined in the following way : 


fQSFY)—wileaNww 0 2 a 
Introducing (15) into (14), we have 
df[dy= kT eee, --p(y)fyf—ufeuN, | (18) 
Or 
affdg=N(Y)f—Fe | + = = ETD 
where 
x(y)=hTjeeeg. p(yly - - - - - - ~ (18) 
and * 
f.=t/euN,=constant. . . . (19) 
Equation (17) is the starting point for the BB outa probe It 
will be referred to as the “‘ detector-equation’”’. It is a so-called Abel 


equation. The solutions cannot be expired in terms of finite 
elementary functions. It is however very suitable for numerical 
integration. For each type of semiconductor the function y(y) is 
completely known and can be tabulated once and for all. When the 
quantity F(y) is found as a function of y by integration of the detector 
equation, it leads to the desired information concerning the electric 
resistance, the capacity, etc., for each value of the current density 7. 
(See § 5.) 

The quantity f(y) has an obvious physical meaning. It is that part 
of the electric field that is needed to maintain the differences in the 
concentration of the free electrons from point to point. 


§ 4. THE BoUNDARY CONDITIONS FOR THE SOLUTION OF THE 
DeEtrectToR Equation. 

At a large distance from the interface, the electron distribution in the 
semiconductor will be the same as if no metal were present. In practice, 
this “large” distance means about 1). According to formula (5) the 
‘concentration of the conduction electrons for z+0o is given by 


ACSC) ae NU ene ee ee Sauer tee (On 
, Yo is the value of y for x= oo. 
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For x+0o the electron distribution becomes homogeneous and in 
consequence f(x 00)+0. Thus 


for y>Yyn, f(yJ>0 and y(y)>0.. . . . (21) 


Substituting (21) in the detector equation (17) one sees that for y>y,,, 
df/dy becomes indefinite. We can abolish this indefiniteness by differen- 
tiation of the numerator and the denominator of the quotient x/f. This 
gives for (df/dy),,, a quadratic equation with the following solutions : 


Nis. steals Je dx 
Gs a rca Diet ee (22) 


f(y) and df/dy are always negative and in (22) the sign must be chosen 
accordingly. The integration of the detector equation can be started at 
Y=Y~. Because no great accuracy is wanted, the integration is readily 
performed. 

For the case in which the current is zero, the constant f, also becomes 
zero and the detector equation reduces to 


dfidye=xif?- @(fors=0), 2 Se ee) 


This equation can be integrated immediately : 


Yo 


f?(y)=2 : x(y) dy eet eed BO (RS) 


Yo 


§ 5. CALCULATION OF THE CHARACTERISTIC PROPERTIES OF THE 
METAL-SEMICONDUCTOR CONTACT. 
In the preceding section it was shown how, for a given semiconductor 
f can be calculated as a function of y for each value of the current 
density i. The next task will be to express the characteristic quantities 
such as the resistance, the capacity, etc., in terms of f and y. 
The conductivity is given by the following formula : 


o=eun,—euN,/yohm-tem.-7. . . . . . (28) 
The barrier-layer resistance R, is defined as the difference of the actual 


resistance per cm. and the resistance per cm.” such as it would be if no 
metal were present. Thus 


r co 


Ro= ay, | Wye) de ve (26) 
From (12) it follows that 
da=kT /e . dy/yf a Ng ete earn ook) 
and thus: 
kT (%@ (Y—Yeo) 


22 eS Se Ty ONNLICIOkey Ta A ee oma 8 
Pe ceuN yy Uh 4 ) 


y, is the value of y at ~=0. 
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Multiplying R,, by 7 one finds the barrier-layer tension U,. U, is that 
part of the applied voltage that appears across the barrier-layer. 


kTi (¥@ (Y—Yoo) ) 
eUN, y, yf 4 


The electric field strength at the surface of the metal F, follows from (15): 


U,= 


F,=f(y=4yo) tiy,feuN, voltem.-4.. . . . . (380); 
The barrier-layer capacity C, is given by 
C,=€€) dF, /dU, farad. cm. 1s ee aes ae 


and follows therefore from (29) and (30). 


§ 6. THE VALUE OF y) AT THE ELECTRODE. 


In order to calculate the characteristic quantities of the blocking layer 
with the help of the formule of the preceding section, the value y) at the 
electrode must be known. This is synonymous with the assertion, that the 
concentration of the free electrons in the semiconductor at the electrode 
must be known. 

The electrons in the semiconductor at the boundary of the metal are 
in thermodynamical equilibrium with the electrons inside the metal. It 
is assumed that this equilibrium is not disturbed even by the largest current. 
The equilibrium concentration of the free electrons at the interface is 
determined by the threshold energy e¢ for the transition of electrons from 
the metal to the semiconductor. As the top of the occupied energy 
levels in the metal coincides with the free energy per electron g introduced 
in (1), this threshold energy is kT In yp, as can be seen from (4). Thus 


eb=hMGnyys oe 6) tes ne 


Now it is known that under the influence of an electric field F, at the 
boundary of the metal, the threshold energy decreases by e/(—eF)/ee9). 
(See, for instance, Handbuch der Physik, xxiv., 2, p. 437.) 


This is a consequence of the “image force”. If Y, is the value of y 
at the interface in absence of any electric field, then : 
€4/(—eH [een )=e(d) —b)=KT In Yofyy. . . . . (33) 


Eliminating Fy from this equation with the help of formula (15) one gets 


fotiyo/euN,=—h? Tce, /e3 . In? (Yo/¥0). (34) 


The quantity kT In Yy) might be called the “true” threshold energy. 
It is a property of the combination of the particular metal and semi- 
conductor under consideration and must be known a priori in order to: 
carry through the preceding calculations. ~The equation (34) describes 
the quantity fp as a function of yy). The point of intersection of the curve 
representing f, as a function of y) (formula (34)) and the curve representing 
f as a function of y (detector equation, formula (17)) gives the values of 
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Yo and to which are needed to apply the formule (28), (29) and (30). 
With this result the mathematical solution of the blocking-layer problem, 
with due allowance for the field effect is obtained. 

In the next section the method will be demonstrated with the aid of 
Some examples. In the following section it will be shown that for large 
_ back currents an asymptotic solution of the detector equation exists 
which enables us to calculate the characteristic quantities of the blocking- 
layer without numerical integration. This asymptotic solution will then 
be used in the discussion of the properties of the blocking-layers at large 
back currents. 

§ 7. SomE EXAMPLes. 

In the following examples the values of the different constants, such as 
the dielectric constant « and the electron mobility u, are not specific for 
a particular material, but do represent plausible values and are chosen 
in such a way as to make the differen coefficients of the detector equation 
simple numbers. 

For N, we take the value 10-! cm.~* and for w the value u=6-3 em.? 
V-tsec.-1, which makes euN,=10. For « the value «—4:55 is chosen 
and hence kT/eeg=10-§ V2cm. for T=293° K. 

The following examples refer to semiconductors of the types A and C. 
For both, the ‘“ true ” threshold energy ed, is taken as 0-349 eV., which 
gives Y>=10% at T=293° K. For the semiconductor of type C the value 
AE is arbitrarily chosen as 0-582 eV. (=23 kT). The concentration of 
the electron donor impurity atoms and the concentration of the electron 
traps are chosen in such a way as to give both semiconductors at room 
temperature the same conductivity when no metal is present. This 
means that both will have the same y,,. For y,, we take the value 104. 
At a large distance from the metal, that is to say for y=y,, c=10-% 
ohm cm. (see formula (25)) and the space charge denisty p is zero. ‘Thus 
formula (7) for the semiconductor of type A gives N,;=10!? cm. and 
for the semiconductor of type C formula (9) gives N;—0-6N,—101!° cm.~*. 
If the arbitrary value N,=1018 cm.-* is chosen for N,, it follows that 
N,=6-01 10!7cm.-*. The numerical integration of the equation (17) 
has been performed by means of a very straightforward method, due to 
Horst van Sanden (1946). The integrations are carried through for a 
number of positive and negative values of ¢. 

In the Tables I. and II. some of the results are summarized with the 
purpose of demonstrating that two semiconductors, which are alike from 
the electrical point of view, may behave quite differently as rectifiers in 
contact with the same metal even when the threshold energy is also the 
same. For the type C semiconductor the blocking-layer resistance is 
much less than for the type A semiconductor and for small current 
densities the direction of the rectification is opposite in the two examples, 
as can be seen from the current-voltage characteristics in fig. 1. These 
differences are brought about by the much larger space charge density 
in the case of the type © semiconductor and its strong current dependence. 
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TABLE I. 
Typo p= 0s mN y= 10. cle 
a R Fy Po bo—¢ 
Acm.-2. Ohm cm.? Yo Vo cm-* el. charges volt 
+10 t-Oo oc." 1-00 x 108 +9-02 x 10° 9-9 1014 0 
+1 4-36 x 10-? 2-45 «10° —3-20 x 108 9-9 1014 0-036 
0 6-40 x 10-2 1-48 x 105 —5-93 x 10° 9-3x 1014 0-048. 
—] 8-10 x 10-2 8-06 « 104 —9-80 x 103 8-8 x 104 0-064 
—10 5-40 x 10-2 2:23 104 2-3 LO: 5-6 x 1014 0-095 
TABLE II. 
Type C: Y,=108; N,=6-10!" cm.*; N,=10% cm.*; AE=—0-582 eV. 
a R Fy Po by—¢ 
Acm.2 Ohmecm.? Yo V. cm.—* el. charges _-volt 
+10 2°38 x 10-8 5:70 x 104 —1-34x104 7-54x101§ 0-073 
+1 1-65x10- 2-93 x 104 —2-02'x 104 4:64101% 0-089 
0 1-60 x 10-8 2-76 x 10# —2-08 x 104 4:36x101 0-091 
—l 1°53 107% 2-60 x 104 —2-15 x 104 4-11x101 0-094 
—10 Talo 10m 1-80 x 104 —2-62 x 104 2:49x101§ 0-101 
Fig. 1. 
10 
5 
typec -O02 -OO! VOLT 
typed -OS ee -O3 = o -O) fe) o2 VOLT 
r =5 
type a 
-lIO 
amp.cm.~ 


The barrier-layer tension as a function of the current density of two crystal- 
rectifiers of the A-type and the C-type respectively. 
For the A-type: Y,=10°; N,=1015 cm.-. 
For the C-type: Y,=10°; N,—6:10'? com.*; N,—101!8 em-3- 
AEK=0-582 eV. 


§ 8. AN ASYMPTOTICAL SOLUTION OF THE ABEL EQuattion. 

Consider the case in which the constant f, in the differential equation (17) 
is negative. Physically this means that a back-current flows. When y 
is increased to infinity, the function y=p(y)/y, tends to zero as can be 
seen from the formule (7), (8) and (9). f being always negative, df/dy 
cannot tend to a finite positive value and therefore y/f cannot become zero 
for large values of y. But then f and df/dy must go to zero with 
increasing y and thus the differential equation reduces to 


S(YW=xMlfe- (35) 
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Substituting (35) into (28) and considering that y>y,,, one finds for the 
back-resistance the approximation 
pape Let Yo dy 
awe ane ),, x(y) 
The above approximation holds as well for large values of y, as for large 


negative values of 7. In this approximation the equation (34) reduces 
to the following transcendental equation for y, : 


(36), 


typ/euN,—=—k?Tce,/e® . In? (Yo/yy). . « . . . (87) 


In the next section the above approximation will be used to calculate the: 
back-current characteristic of the contact of a metal with a semiconductor 
of the A-type. 


§ 9. THE BAcK-CURRENT CHARACTERISTIC FOR A SEMICONDUCTOR 
OF THE A-TYPE. 


For a semiconductor of the A-type, equation (35) reduces to 
| p—PEUNN: (y=ya) 


—_— . 3 y J : 38 
~ Eg y? (38), 
and (29) reduces to 
= — EEO Re Bd yal Ws 2 39 
Us 2eN ,(euN,)? ’ OE Reece EM as ict aril (39), 


From equations (37) and (39) the equation for the back-current characteristic: 
can be derived, computing y, from (39) and substituting in (37). For 
practical calculations, however, it is more convenient to consider 
equations (37) and (39) as a parameter representation of the back-current 
characteristic. For each chosen value of y, these equations give the 
corresponding values of 7 and U,. 

When the back-current increases, the value of y) decreases, as can be 
seen when computing di/dy from (37). Thus in (39) the factor (yj—y%,) 
will decrease when 7? increases and the possibility of a maximum in the 
border-layer tension exists. With the aid of (37), ¢ can be eliminated 
from (39). For the occurrence of a maximum for U, dU,/dy) must be 
- gero. This leads to the following transcendental equation for yj, the 
value of y, for which the barrier-layer back tension attains its maximum : 


iri NEO OC 0 ei eee ee A) 
Substituting the value of y* found from this equation into (37) one finds. 
the corresponding value of the back-current 7*, and substitution of yj 
and i* into (39) gives the maximum U# of the barrier-layer tension. 
Equation (40) can be solved graphically in a simple way. 
Some examples : 


N; Yo Yo UF * 
(a) ulO* cem.-> 10% 2-03 x 104 —3:7l volt —30 A/em.? : 
(b) 104% cem.-3 1-410! 2-62 x 104 —100 volt —107 A/cm.? 


(c) 1015 em. 1012 3-07 104 —279 volt —157 A/cm.*: 
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The voltage maximum in the back direction cannot be observed on 
selenium or cuprox rectifiers because at current densities much smaller 
than those at which the maximum would occur, the rectifier would have 
been disturbed by heat. The maximum is observable on crystal rectifiers 
when high-frequency voltage pulses are applied (Torrey and Whitmer 
1948). The maximum is found at a voltage of some 100 or 200 volts. 
‘The value of the current density in crystal rectifiers of the whisker-contact 
type is difficult to estimate, because the contact area between the whisker 
and the crystal is not exactly known. At any rate, the back-current 
densities at high back voltages can be estimated at magnitudes of a 
hundred A/cm.’ or so. 


Fig. 2. 


{6 =P) 54 
«x!O  forad cm 


=! 2 — 3. =36/ 4 
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1/C; of a crystal rectifier of the A-type as a function of the back-voltage U,. 
Y¥o=105> N=10) em. 


The example (a) corresponds to a rather bad rectifier ; the example (b) 
corresponds to a modern high-inverse-voltage crystal detector and the 
example (c) corresponds to a rectifier of the selenium type. 

The barrier-layer capacity C, follows from (31) : 


dE (dU p= (dE y/dyo)/(dUz/dy), .« . . . « (41 ) 
d¥y/dy) follows immediately from (33) and dU,/dy, follows from (39) 
when 2 is eliminated with the aid of (37). 


As an example 1/C; is calculated as a function of U, for a semiconductor 
of the A-type with N,—1015 cm. and Y,=10’, (see fig. 2). 
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For small back-voltages 1/C? is the same linear function of the back- 
voltage as has been found by Schottky (1941-1 942). 

For large back-voltages deviations from linearity are observed which 
agree qualitatively with the calculated curve of fig. 2. The curve has a 
maximum analogous to the maximum in the current-voltage characteristic 
in the back direction. 

The maximum occurs at a back-voltage U** corresponding to a value 
yo* in the parameter representation of (37) and (39), i.e. to the value 
which is a solution of the transcendental equation 


In Yolyg*=3(Yo*—YoYeoe - - - + + (42) 


In the example of fig. 2 one finds y#*=2-35 x 104, and substituting this 
value into (37) one finds 1**——25-3 A/em.-*. Substituting y** and i** 
into (39) one finds U**=—3-67 V. Just as with the maximum in the 
current-voltage characteristic it is to be expected that the maximum 
in the (1/C;, U,) curve will only be observable with short voltage-pulses. 
As far as the author is aware, such a maximum has not yet been observed f. 


§ 10. THE RECTIFIER OF THE SELENIUM-TYPE. 


Among the crystal rectifiers used in practical applications two quite 
different types can be discerned: the selenium and the cuprous oxide 
rectifiers on the one side and crystal detectors on the other. The principal 
difference lies in the fact that the crystal detector operates at 10* or 105 
times as large a current density as the selenium or the cuprous oxide 
rectifier. That the applied current density is of preponderant influence 
on the unipolar conductivity of a metal semiconductor contact, will be 
clear when comparing, for example, fig. 1 and fig. 3. Most semiconductors 
in contact with a fine metal point, have an asymmetrical current-voltage 
characteristic, when the current density is made sutficiently high. This 
is effected by touching the crystal very loosely, so that the contact area 
is very small. Extremely few semiconductors have good rectifying 
properties at as low a current density as about a milliampere per cm.?. 
Only two are used in practice, viz. selenium and cuprous oxide. 

For a good rectifier of the selenium or cuprous oxide type, it is claimed 
that at a back-voltage of about 10 volt, the back-current will not exceed 
one milliampere per cm.”._ The contact being characterized by the two 
parameters Y, and N;, one can ask for which pairs of values of these 
parameters a metal semiconductor contact will satisfy this requirement. 
In the formulz (37) and (39) we substitute for 7 the value —10~* amp em,~* 
and for U, the value —10 volt and thereupon these equations are solved 
for Y, and N,. This is performed by putting N; arbitrarily equal to 
integral powers of ten and computing the corresponding values of Yo. 
The pairs of values thus obtained, together with some other calculated 


ee ____ ee 

+ Note added in proof.—R. N. Smith (1946, Phys. Rev., 69, 683), reports that 
the contact capacity of Si-detectors, plotted as a function of the voltage, shows 
a minimum at a backvoltage of 0-5 volt- 
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features of the contact, are tabulated in Table III. Columns | and 2 
give N; and Yy. Column 3 gives the value of 4 corresponding to Yo. 
Column 4 gives the value of yo. Column 5 gives the electric field at the: 
electrode. Incolumn 6 V, the so-called “ diffusion-potential ” is tabulated, 
as derived from 

V =AT In Y,/y.=0-058 log Yy/y., volt. - = 7 smu(43) 


Column 7 gives the conductivity o of the semiconductor. 

It must be kept in mind, that the figures in this table are not calculated 
with values of «, w and m* which hold exactly for selenium ; these figures. 
are therefore only meant as an illustration. The quantities that are open 
to experimental check are: N; (from the voltage dependence of the 
capacity); V, (from the current-voltage characteristic in the forward- 
direction) ; and o (also from the current-voltage characteristic in the 
forward-direction). The figures inside the framed part of Table ITI. agree: 
with experimental observations on actual rectifiers. 


TasieE ITI. 
S ae do in RG in Vqin o in 
N,in cm. Yo ‘volt Yo volt em.—} volt Ohm? cm.-t 

1016 4-96«1078 1-66 2-81«1019 2-81 «108 1-61 1 

1017 1-39 1029 1-17 8-88x 10° 8-88 «x 105 1-05 10 
1016 1531014 0-82 2-81 109 ZeOle< 102 0-05 102 
1015 1:-47<x1022 0-71 8-88x 108 8:-88x10* 0-48 108 
1024 1:84 102° 0-60 2-81 108 Peel SNe Oea3il 104 
1038 9-24«108 0-52 8-88x107 8-88x10% 0-17 10° 
1012 1:16<108 0-47 2:98 107 2:98x108 0-062 106 
1011 1-46«109 0:53 1:09 108 109x104 0-068 107 


In fig. 3 the back-current characteristic is drawn for a rather arbitrary 
example, as calculated according to §9. The decreasing back-resistance. 
with increasing back-voltage is a common experience on selenium rectifiers. 

Surveying Table III. the most striking feature is the small range for the 
threshold potential ¢,) in which rectifiers are actually realized. 
This is the more astonishing, as good selenium rectifiers can be made with 
a rather large variety of electrode metals. Now, it is known that every 
selenium rectifier gets its rectifying properties not until it has been 
subjected to an electrical treatment. Usually a constant back-current 
is forced through the rectifier, the back-voltage increasing in the course 
of the treatment until the desired back-voltage is reached. The question 
arises as to which of the two characteristic parameters Y, and N, is 
affected by this treatment. It has been suggested (F.I.A.T. Report), 
that during the electrical treatment under the influence of the electric 
field the impurity atoms, which are assumed to be electrically negative, 
are driven from the semiconductor to the metal. There they should either 
be stored in a very thin layer at the interface, or be taken up by the metal, 
but in all events to be put out of action. A reduction of the impurity 
atoms will increase the back resistance. 
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The curve of fig. 4 represents the back-voltage for a back-current of 
- 10% A/em.~*, as a function of the concentration of the impurity atoms 
at an arbitrarily chosen value of Y,. As can be seen from this curve, 
which would describe the falling off of the concentration N, during the 
electrical treatment, the concentration of the impurity atoms must be 
diminished by a factor ~102, to compensate a “ mismatch” in the threshold 
potential J) of ~10-1 volt. Curves as of fig. 4 have been calculated for very 
different values of Yo, leading to practically the same result. This result 
makes it rather improbable that the effect of the electrical pretreatment is 
due to a diminution of the impurity atoms. The scanty experimental data 
concerning the change of N , during the electrical treatment show a decrease 
of N; during the last stage of the electrical treatment of a factor 5 at most. 


Fig. 3. 


10 amp. em* 


i i i back- t of a crystal 
The barrier-layer tension U, as a function of the back current of 
‘ rectifier of the A-type. N,=1015cm.-?; Y)=1-47x10". This example 
is meant to represent the behaviour of a normal Selenium rectifier. 


It is conceivable that the observed influence of the electrical pretreatment 
on the concentration of the impurity atoms is only a secondary effect 
and that the primary effect is an increase of the threshold potential Po- One 
might for instance suppose that during the treatment under the influence 
of the electric field at the boundary of the metal, metal ions enter the 
selenium, building up a double layer at the interface, by which the thres- 
hold potential 4, for electron holes is increased. ue this sesump Hon ie 
curve of fig. 5 is calculated for N;=10'? cm.~ and 1=10 av : 
It appears that an increase of ¢, of 0-1 volt only, can raise the back-voltage 
to an enormous.extent. It is from present day experimental knowledge 
not possible to decide between the two suggestions. In virtue a the 
foregoing considerations the author attaches the greatest probability to 


the last. 
. BG 2 
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From Table III. it is evident that the chance that an arbitrary combina- 
tion of a metal and a semiconductor will be a good rectifier, is extremely 
small, and so is the number of the actual technically useful combinations. 
The copper-cuprous-oxide rectifier is a case of fortunate coincidence, 
very sensitive to traces of admixtures and the thermal pretreatment. 
The great technical success of the selenium rectifier is due to the possibility 
of adapting it to its purpose by the simple electrical pretreatment, and this 
possibility is linked up with the extraordinary mobility of foreign ions in 


Fig. 4. 
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The back-voltage at a current density of 10-3 Axcm.- as a function of the 


ian of the impurity atoms N; for a rectifier of the A-type, with 
Om : : 


hexagonal selenium (Torrey and Whitmer 1948). But this merit is also 
its weakness, for the electric fields during the normal operation of the 
rectifier give rise to the same effect as those during the electrical pre- 
treatment and this results in a sometimes very appreciably increasing 
resistance both in the forward and in the back direction. 

In connection with the above considerations the effect of artificial 
barrier-layers, which are sometimes applied, for instance in the form of a 
thin film of lacquer between the selenium and the metal, can easily be 
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explained. When the thin lacquer film is pierced at several spots, the 
local current density in the different parallel current paths is much higher 
than it would have been, when the same current was homogeneously 
distributed over the whole cross-section of the rectifier. It follows from 
the shape of the current-voltage characteristic that the rectifier-effect is 
strongly favoured by an increase of the current density. The effect of 
the pierced insulating film is thus essentially the same as that of a loosely 
applied whisker contact. 


Fig. 5. 
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The back-voltage at a current density of 10-? Axcm.~* as a function of the 
threshold potential ¢, for a rectifier of the A-type with a concentration of the 
impurity atoms N; of 101° cm.~%. 
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SUMMARY. 


In examining the magneto-resistance of thin fibres of bismuth, which 
are known from previous work to approximate to single crystals, it was 
found that a considerable proportion of these, about 10 per cent, showed a 
decrease in the resistance on application of a magnetic field. The effect is 
only present in fields of the order of 1000 oersteds or less and is of the order 
of 10-4p). The influence of direction of current, direction of field and orien- 
tation of the crystal axes is examined. It is considered that these observa- 
tions provide the first unequivocal evidence of the existence of a longi- 
tudinal Hall coefficient. Used in conjunction with the theory of conduction 
in anisotropic solids, it is possible to indicate the appropriate crystal 
orientation. It has not, as yet, been possible to designate this with 
accuracy. 


§ 1. INTRODUCTION. 


In the first paper of this series (Donovan and Conn 1949), the initial 
results of an investigation into the magneto-resistive behaviour of bismuth 
fibres were discussed, and particular attention was paid to the eluci- 
dation of the structure of the specimens. Comparison of the results with 
those obtained for single crystals (Stierstadt 1933) led to the conclusion that 
the fibres were, with certain exceptions, single crystals with the main 
cleavage plane parallel to the fibre axis. The exceptions were provided by 
the low-resistance (7.e. the thicker) fibres and it seemed feasible that these 
should possess a less highly developed crystal structure; their magneto- 
resistive behaviour was not radically different from that of the high- 
resistance fibres, but was sufficiently distinctive to enable positive identi- 
fication. to be made. 

In the course of this investigation certain anomalous results were 
obtained and it is the purpose of the present paper to discuss these. The 
essential feature is that, under certain circumstances, the change in resistance 
resulting from the application of the magnetic field may be negative, or, in 
other words, the change in resistance may reverse in sign when the direction 
of the magnetic field is reversed. The possibility of such an ‘“‘ Umkehreffekt”’ 
has been discussed in a number of theoretical papers, but no conclusive 
experimental evidence has hitherto been reported. The work of Kapitza 
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(1928) and others has lead to the conclusion that the change in resistance 
is proportional to the square of the field strength in weak fields and tends 
gradually to a linear relation at very high field strengths. The results 
discussed below indicate that the resistance cannot be represented simply 
as a quadratic function of the field strength and suggest that a linear term 
is also present which outweighs the quadratic term below a certain 
critical field strength. A preliminary notification of the present work has 
already appeared (Conn and Donovan 1948). 

This paper does not contain a detailed treatment of the ‘‘ Umkehreffekt”’ 
as such, but presents experimental results which may be regarded as 
evidence of its existence in bismuth; in the absence of sufficiently extensive 
experimental data, a complete interpretation cannot be attempted, but it 
is felt that a qualitative explanation can with advantage be given. 


§ 2. EXPERIMENTAL RESULTS. 


The preparation and mounting of the fibres, and the experimental 
technique employed, have been described in the previous paper (Donovan 
and Conn, loc. cit.). The same potentiometer arrangement was used in 
the present series of measurements and it was again found convenient to 
use the galvanometer as a continuous recorder of potential difference. 
In the earlier experiments, the change in resistance of the specimen was 
determined as a function of its orientation with respect to the magnetic 
field, which remained constant for a given “ cycle.’’ In the experiments 
to be described the fibre was maintained in a fixed position with respect 
to the field and its resistance was obtained as a function of the field strength 

At low field strengths the residual magnetism of the pole-pieces presents 
difficulties and it was deemed advisable to employ a source of magnetic 
field which did not involve the use of iron cores. Two large coils were 
available, each wound with 1790 turns, the outer diameter being 36 cm. and 
the inner diameter 12 cm. These were supported coaxially, as in an 
Helmholtz arrangement. With a separation of 3:8 cm. a field of 120 
oersteds per ampere was obtained at the centre of the gap. Owing to the 
small dimensions of the specimens, no difficulties arising from inhomo- 
geneities of the field were encountered. 

With this arrangement the procedure of varying the field was consider- 
ably simplified. 

It will be convenient to refer to the ¢ and 0 orientations already intro- 
duced (Donovan and Conn, Joc. cit.) and fig. 1 has been reproduced to make 
the notation clear. In both cases the axis of rotation 1s perpendicular to 
the magnetic field; the axis of the fibre is parallel to, and perpendicular 
to the axis of rotation in the ¢ and 0 cases respectively. In the latter case, 
therefore, the fibre rotates in a plane parallel to the magnetic field. 

In the ¢ case the anomalous behaviour has appeared in certain fibres and 
is illustrated in fig. 2 which shows the change of resistance plotted ae 
field strength for a fibre whose orientation corresponds to the meget 2 
the ¢ curve. The fibre had a zero-field resistance of 170 ohms. It will be 
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seen that for field strengths less than approximately 350 oersteds the change 
in resistance is negative and the minimum value occurs in the neighbour- 
hood of 180 oersteds. The decrease in resistance is extremely small but, 


Higseds 


Isometric drawing of the holder. AB is the axis of the fibre and the ¢-axis. 


Fig. 2. 
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Variation of resistance with field-strength for 170 ohm fibre. 


Electrical Conductivity of Bismuth Fibres 773: 


by virtue of the comparatively high resistance of the specimen and the- 
sensitivity of the arrangement, it is easily detected and is readily repro- 
ducible. 

It should be emphasized that this phenomenon is not characteristic of 
bismuth fibres per seq, and in fact only a fraction of the order of 10 per 
cent of the fibres exhibited this effect. For this minority the effect was. 
quickly and convincingly demonstrated. Having balanced the potentio- 
meter in the absence of a magnetic field, a suitable permanent magnet is 
brought up to the specimen and the galvanometer observed. 

All the fibres which showed the negative change in resistance did so. 
when the orientation corresponded to the minimum of the ¢ curve. With 
one exception, which will be referred to later, no such effect was observed. 


Biss 3: 


{200 
H Coersteds) 
Variation of resistance with field-strength for 6 ohm fibre. 


when the orientation corresponded to the maximum of the ¢ curve. The 
significance of these two extreme positions lies in the fact that they 
represent simple orientations of the main cleavage plane with respect to the 
field (Stierstadt 1933, Donovan and Conn loc. cit.). No negative change 
in resistance was observed in any fibre when its axis was parallel to the 
field (9-orientation). 
Considerable variation was apparent in the order of magnitude asso- 
ciated with the anomaly. In particular, the range of field strength over 
which the resistance change was negative varied widely from one specimen 
to another; the limits of the critical or “change-over” field strength were. 
roughly 300 to 1000 oersteds. No useful purpose will be served ah 
including a large number of experimental curves, but a further yaaa e 
may be given (fig. 3) illustrating a wider range of negative resistance change.. 
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The critical field strength in this case is approximately 930 oersteds. 
The fibre had a zero-field resistance of 6 ohms and again was set at the 
minimum of the ¢ curve. 

The magneto-resistive behaviour of the fibres which did not exhibit the 
effect in question was quite consistent, and is illustrated in fig. 4, which 
was obtained with a representative “‘ orthodox ” specimen. The change 
of resistance for weak fields is shown for orientations corresponding to the 
maximum and minimum of the ¢ curve 

Finally, we give curves (fig 5) obtained with a fibre which was unique 
in its behaviour and which provided the solitary exception previously 
mentioned. This fibre had a zero-field resistance of 46 ohms and was 
remarkable in two respects: (a) the decrease in resistance was about 
fifty times larger than that obtained in other cases and (b) the decrease in 


H(oersteds) 


“* Orthodox ”’ magneto-resistance effect in weak fields. The upper and lower 
curves represent the maximum and minimum respectively of the ¢ cycle. 


resistance was obtained with orientations corresponding to the maximum 
as well as the minimum of the ¢ curve. The critical field strengths in 
these two cases are approximately 1300 and 1700 oersteds respectively, 
and these are considerably greater than any other observed values. 


§ 3. PossrBLE ERRORS. 


In view of the small magnitude of the decrease in resistance, and the 
nature of the phenomenon, suspicion is at first to be directed towards 
possible sources of experimental error. The various effects, however, 
which might conceivably give rise to spurious results can be discussed, if 
not with finality, at least with sufficient cogency to enable the ‘““‘Umkeh- 
reffekt ” to be regarded as genuine. 


Hlectrical Conductivity of Bismuth Fibres 7715 


The possible sources of error which are relevant in magneto-resistance 
measurements have been enumerated by Thompson (1938). Of these 
the most serious is undoubtedly due to the presence of the Ree eet 
effect. If the electrodes are attached to the specimen in such a way that 
they do not lie on a line of current flow, then, when the magnetic field is 
applied, a component of the Hall e.m.f. will appear across the terminals 
Such a misalignment is illustrated, in exaggerated form, in fig. 6 (a) whieh 
represents a section through a specimen perpendicular to the magnetic 
field. The Hall e.m.f. reverses, of course, when the field is reversed and the 
maj ority of previous workers have used mean values obtained from the two 
‘directions of the field. Under these circumstances the resistance change 
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Variation of resistance with field-strength for 46 ohm fibre. 


must necessarily be a symmetrical function of the field and no ‘“ Umkeh- 
reffekt ’’ could have been detected. Moreover, it should be pointed out 
that the decrease in resistance is associated essentially with weak fields and 
that, in the past, more attention has been paid to investigations in strong 
fields, or at low temperatures, where the resistance changes are very large. 

As mentioned in the previous paper already referred to, specimens in 
the form of fibres possess the distinctive advantage of combining the struc- 
ture of a single crystal with a resistance value suitably high for accurate 
“measurements. In addition, the shape of the specimens is such as to 
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minimize the influence of transverse effects, and in this respect fibres 
are superior to the previous types of single crystal, such as triangular 
prisms (Stierstadt 1933) and plates (Thompson 1938). Two other reasons. 
may be adduced for believing that the Hall effect referred to above is not 
responsible for the observed negative resistance change. Firstly, the 
electrodes are not “ localized ”’ at points on the specimen, as would be the 
case for instance with fine wires ; the method of construction ensures that 
contact is made over the whole cross-section, as shown in fig. 6 (b). Thus 
any effects due to the type of misalignment depicted in fig. 6 (a) are elimin- 
ated. Secondly, the effect is not universal and it is difficult to see how such. 
a general explanation as the Hall effect can be valid when about 90 per 
cent of the fibres show no indication of the phenomenon. 
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(2) Possible misalignment of specimen. Magnetic field 
perpendicular to plane of paper. 
(6) Type of mounting used for fibres. 


The evidence indicates an explanation in terms of the crystai structure 
of the specimen, and suggests that the ‘‘ anomalous ”’ fibres are charac- 
terized by departure from the general structure discussed in the previous. 
paper. Whereas in the majority of fibres the main cleavage plane is parallel 
to the fibre axis, it seems reasonable to suppose that in the “ anomalous ” 
cases this no longer applies, and the main cleavage plane is inclined at a 
small angle to the axis. Small variations in this angle might well account 
for the different values of the ‘‘ change-over ”’ field strength observed. 

It is natural to seek confirmation of such conclusions by examining the 
X-ray diffraction of these fibres. A preliminary investigation undertaken 
through the kind co-operation of Mr. R. B. Shepherd has served to clarify 
the problem but it cannot be said as yet to have yielded clear evidence. 
The inclination of the cleavage plane to the axis, upon which our inter- 
pretation is based, is sufficiently small to make discrimination difficult 
by X-ray methods and, bearing in mind that the fibres are of the order of 
10~* cm. in diameter, it would appear that this corroboration, which is 
clearly desirable, is likely to require a separate investigation in itself. 
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§4. OUTLINE oF Previous Work. 


Before any discussion of the foregoing experimental results is attempted, 
it is necessary to give a brief survey of the existing literature. The 
theoretical approach may be directed along one of two lines: the starting 
point may be either the classical equations of conduction in an anisotropic 
medium, or the quantum mechanical formulation of transport phenomena. 
The latter method of treatment has not, hitherto, proved particularly 
satisfactory in the case of bismuth, as the effect of the crystal symmetry 
has not been rigorously investigated. 

The negative change in resistance would be readily explicable if the Hall 
effect contained a longitudinal as well as a transverse component, i. ¢. if an 
e.m.f. were produced parallel to the current flow. This must not be con- 
fused with the spurious Hall effect due to a bad experimental arrangement 
considered above. Such a longitudinal Hall effect would give rise to a 
resistance change whose sign depends on the direction of the magnetic 
field. The transverse character of the Hall effect was postulated by Voigt 
(1928) but the consequences of removing the restriction were examined by 
Kohler (1934) following some experimental results on the Hall effect 
(Verleger 1932) which could not be reconciled with Voigt’s treatment. In 
view of the considerable experimental evidence supporting the Voigt 
theory, Shoenberg (1935) was led to suggest that Verleger had incorrectly 
described his experiments and that in fact the discrepancy was not 
genuine. Shoenberg dealt with some of the implications of Kohler’s work 
and endeavoured to obtain experimental evidence by means of a magneto- 
resistance experiment with a bismuth crystal. No indication, however, of 
a longitudinal Hall effect was forthcoming. 

A more detailed theoretical investigation of the magneto-resistance 
effect, based on Kohler’s work, was carried out by Meixner (1939). Starting 
with the general expression for the resistance, in terms of odd as well as even 
powers of the field strength, Meixner showed that an “ Umkehreffekt ” 
could occur in anisotropic substances, for certain orientations of the cur- 
rent, field and crystal axes. In bismuth, when the conditions for max- 
imum ‘“‘Umkehreffekt”’ are fulfilled, the linear term should out-weigh the 
quadratic term for field strengths less than about 6700 oersteds and the 
change of resistance should be negative. 

In a later paper Meixner (1941) made a comprehensive examination of 
electrical transport phenomena, from the quantum mechanical standpoint. 
This treatment, which was quite general, led to the result that neither the 
electrical nor the thermal conductivity was dependent upon the direction 
of the magnetic field. Any ‘“‘ Umkehreffekt ” was thus inexplicable on 
this basis, but, in surveying the experimental results on the Hall effect in 
bismuth, Meixner admitted that the question must be regarded as still 
open. Some significance is to be attached to the fact that an ‘‘ Umkeh- 
reffekt ” of thermal conductivity, also excluded by Meixner’s theory, has 
been observed in bismuth by Griineisen and Gielessen (1936). It would 
appear that the quantum mechanical formulation necessarily involves 
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assumptions that render it invalid in the case of bismuth. In a further 
application of the electron theory of metals, Kohler (1941) suggested that a 
possible explanation of these discrepancies might lie in the neglect of the 
quantization of the electron paths in the magnetic field. 

The only previous experimental evidence of this effect for bismuth has 
been reported by Casimir and Gerritsen (1941). Their measurement 
showed that the resistance was actually a function of the direction of the 
field, for a single crystal at 77:3°K. Casimir and Gerritsen concluded, 
however, that the cause was spurious, since a similar “‘ Umkehreffekt ” 
was found with a crystal orientation which was precluded on theoretical 
grounds. No reference was made to Shoenberg’s paper and, as no details. 
of the experimental arrangement were given, the interpretation of the 
results cannot profitably be discussed. 


§ 5. Discussion oF RESULTS. 

We examine now the implications of the relevant theoretical results. 
(Shoenberg 1935) for various special cases and conclude by advancing a 
tentative explanation to account for the results in the second section of 
this paper. 

The most general expression for the Hall effect is 

E =R jz; Ji, 158 . . . . . . . (1) 
where E, J, H are vectors denoting the electric field responsible for the Hall 
e.m.f., the current density and the magnetic field respectively. The R,,, 
are governed by the symmetry of the crystal and in the case of bismuth it 
has been shown by Kohler (1934) that there are only four independent 
constants : 

Rosi=—Ry32=Ry, | 
Rsi2=—R3.1:=Re, 
Rygg=—Roi3=Rs, 

Run= —Ri22= —Rogi9=—Roo1=Ry. 


If there is no longitudinal Hall effect, then R ,,,=—R,,;, and it follows. 
from (ii) that 


(ii) 


Ry=Re and Rie=0) 2s Se Pe eG} 
In the subsequent application to the bismuth crystal, we shall use a 
rectangular cartesian system of axes with OZ parallel to the principal 
(trigonal) axis and OX parallel to one of the binary axes. The main 
cleavage plane is thus the (x, y) plane. In the experiments which concern 
us, the magnetic field is invariably perpendicular to the current flow. For: 
this case Shoenberg has shown that a longitudinal Hall effect, which is. 
‘oto in equation (i), produces a change in the specific resistance given 
ue 
Ap=H[(R,—R,) cos « sin « sin y 


—R, sin? «(sin y sin 3f--cos y cos « cos 3 L)|, Sea aaiivy 
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Here « is the angle between the principal axis and the current direction, ¢ 
is the angle between the binary axis and the projection of the current 
vector on the main cleavage plane, and y is the angle between the magnetic 
field and the projection of the principal axis on the plane normal to the 
current vector. The significance of these angles can be seen from fig. 7 (a). 
Expression (iv) represents a magneto-resistive effect which is linear in 
H, and the change in the observed resistivity when the field is reversed 
should be 2Ap. Ifthe current flow is perpendicular or parallel to the main 
cleavage plane («=0° or 90°) then the effect vanishes, as is also the case 
in polycrystalline material. 

In Shoenberg’s own experiment «=45°, ¢=17° and y was varied by 
rotating the magnet. Reversal of the magnetic field produced no detect- 
able change in resistance for any value of y and the essentially transverse 
nature of the Hall effect was regarded as established. However, only one 
value of magnetic field (13,700 oersteds) was used and it is probable 
that this was too high for the effect to appear. The effect is exceedingly 
small and is clearly more easily detected in weak fields where the orthodox 
magnetic resistance is not large. 

In the orientation discussed by Meixner (1939) the magnetic field is 
parallel to one of the binary axes and the current is restricted to the (y, z) 
plane. Thus y=€=90° and equation (iv) reduces to 


Apj/H=(R,—R,) cos « sina +-R,sin?« . . . . (¥) 


For «=45° (see fig. 7 (b)) this expression has a maximum value in agree- 
ment with Meixner. An identical arrangement, with «=45°, was used by 
Casimir and Gerritsen (1941), and produced an apparent ‘‘ Umkehreffekt”’; 
however, such an effect was also obtained when the magnetic field was 
perpendicular to the binary axis, which disagrees with the theory. 

In discussing our own experimental results, the significance of the ¢ 
resistance curves must be bornein mind. Stierstadt (1933) has shown that 
in P, (L) crystals (current flow perpendicular to a binary axis) the resis- 
tance in moderately low fields is controlled by the orientation of the (111) 
plane, i. e. the main cleavage plane. When this is parallel to the field the 
resistance has a maximum value, and when it is perpendicular to the field 
the resistance has a minimum value. Further, in P, (||) crystals, that is, 
current flow parallel to binary axis, the controlling influence on the 
resistance is exerted by the (111) plane, i. e. the subsidiary cleavage planes. 
The behaviour of these planes is opposite to that of the main cleavage 
plane. For example, if one of the subsidiary cleavage planes is parallel to 
the field the resistance has a minimum value. In this case the inclination 
of the main cleavage plane to the field is approximately 71° and therefore 
it is clear that in any P, type crystal the minimum of the ¢ curve represents 
a state of affairs in which the main cleavage plane is roughly perpendicular 
to the field; at any rate the departure must be less than 20°. 

Returning to fig. 7 (b), we see that if the current vector J is parallel to the 
y-axis we have a normal P, (L) crystal with a main cleavage plane 
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parallel to the field, 7. c. at the maximum of the ¢ curve. IfJ is restricted 
to the (a, y) plane, any type of P, crystal can be represented and as «90° 
in every case no “ Umkehreffekt ” is possible. As it has already been 
established in the previous paper that in these bismuth fibres the main 
cleavage plane is parallel to the axis and hence parallel to the current flow, 
we see why the majority of fibres do not exhibit this effect. 

The explanation of the anomalous behaviour of the minority described 
in § 2 would seem to lie in the exceptional nature of the crystal structure of 
these fibres. In drawing the fibre the main cleavage plane has become 
tilted with respect to the normal orientation, with the result that the 
‘current flow is no longer parallel to the main cleavage plane and « is no 
Jonger 90°. This displacement must necessarily be quite small and of the 
order of a few degrees. 

We consider first the orientation corresponding to the maximum of the 
6 curve. In P, (L) crystals the main cleavage plane is parallel to the 
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‘Orientation of current J and magnetic field H with respect to crystal 
axes, OZ parallel to principal axis ; OX parallel to binary axis. 


(a) General case ; OM is projection of J on (x, y) plane ; 
ON is projection of z-axis on plane normal to J. 


(5) Special case with J in (y, z) plane and H parallel to OX. J, 
and J, correspond to values of 45° and approximately 85° 
respectively for a. 


field and, as shown in fig. 7 (6), the current vector is displaced slightly 
‘from the y-axis although still remaining in the (y, z) plane. For this case, 
therefore, equation (v) will apply. In P, (||) crystals, and intermediate 
‘cases, both H and J will be inclined to the main cleavage plane and no 
‘corresponding simplification of equation (iv) can be made. 

Next we consider the orientation corresponding to the minimum of the 
¢ curve. In true P, (L) crystals the main cleavage plane would be 
perpendicular to the field; in imperfect P, (L) crystals, that is to say, 
crystals possessing a fault in accordance with the conjecture made above, 
this will no longer by the case, but the principal axis will be in the plane 
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defined by H and J. It may be instructive to consider this orientation, 
which is shown in fig. 8 with the crystal fault intentionally exaggerated. 
If the principal axis lies in the plane of H and J, then y=0 and this restric- 
tion alone causes (iv) to reduce to 


ApJ =—Rjsin’acosacos3¢ . . . . . . (vi) 


The * Umkehreffekt ” is controlled by «, i.e. the magnitude of the crystal 
fault, and by ¢, which is determined by the position of the binary axis. In 
crystals of the quasi-P, (1) type, the binary axis is normal to the current* 
and ¢=90°. Hence the resistance change should be normal. For all 
other orientations of the binary axis expression (vi) has a finite value. 

In crystals of the quasi-P, (||) type the position is complicated by the 
fact that the principal axis is inclined at an angle of approximately 19° 
to the plane containing H and J. For intermediate P, crystals the angle 
will have some value between 0 and 19° and, since «, y and @ are all 
finite for these cases, equation (iv) must be used as it stands. 


Fig. 8. 
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MAIN 
“—~ CLEAVAGE 
PLANE 


PRINCIPAL 
AXIS 


Crystal fault (exaggerated) in which main cleavage plane is not perpendicular 
to magnetic field. The principal axis, current J and field H all lie in 
one plane. 


CONCLUSION. 


The nature of the specimens we have used confers two very real advan- 
tages in exploring the reality of such an “‘ Umkehreffekt.” The prepara- 
tion of crystals of such a relatively high resistance enables one to examine 
with care and precision any departures from the normal behaviour. This 


benefit is to some extent obtained at the expense of exactness in specifying 
Oe a a aa 


* Tn this case fig. 8 represents a section in the (y, z) plane. 
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to which crystalline category any given fibre belongs. Furthermore the 
size and the method of mounting would seem to eliminate the possibility 
that such observed effects were due to a component of transverse Hall 
coefficient. 

It is difficult to resist the conclusion that the evidence here presented 
confirms the existence in these specimens of a longitudinal component of 
the Hall coefficient. Because there is at the moment an inevitable lack of 
precision in our specification of the category of any given fibre, we can go 
no further than the above qualitative interpretation setting out the various 
possibilities which can give rise to an “‘ Umkehreffekt.”” With reference 
to the order of magnitude, it is relevant to quote Casimir and Gerritsen 
(loc. cit.): ‘“‘A linear term would lead to a decrease of resistance at low 
fields but ....this decrease would only amount to 210-4p).” The 
results illustrated in figs. 2 and 3 are in striking agreement with this 
prediction. 
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ABSTRACT. 


Theories of the origin of the dipole components of the Earth’s main 
magnetic field are of two types; distributed theories attribute it to a 
fundamental property of rotating matter, core theories to current systems 
within the core. The variation with depth below the surface of the 
horizontal field intensity is different for the two theories. Experimental 
values are given for this variation, obtained by measurements in coal 
mines, which are near to the values predicted by a core theory, and 
significantly different from those predicted by a distributed theory. The 
magnetic effects of the sedimentary rocks and magnetic anomalies due to 
the basement rocks are shown to have negligible results on the 
measurements. 


Iv has been suggested, notably by Blackett (1947 and 1949), that the 
axial dipole component of the geomagnetic field might be a fundamental 
property of any rotating spherical mass ; the non-dipole part of the field 
(about 15 per cent) arising from subsidiary processes such as induction. 
The other possible theory of, the geomagnetic field suggests that it is 
due to the flow of electric currents maintained by thermo-electromotive 
forces or by dynamo action ; these must be confined to the deep interior 
of the earth, where the electrical conductivity seems likely to be 
appreciably higher than that of the surface rocks. It is convenient to 
call the former a Distributed Theory and the latter a Core Theory. 
Bullard suggested that it might be possible to decide between these 
two theories by measurements of the change with depth inside the earth 
of the horizontal and vertical components of the geomagnetic field. 
Preliminary measurements were made in mines in South Africa (Hales 
and Gough 1947) and in Lancashire (reported by Chapman 1948 a). 
Much more extensive measurements have now been made in the coal 
mines of this country and methods of establishing the significance of the 
results, which were lacking in the earlier work, have been developed. 


a ee a a ee =e Se Ra le 
* Formerly at the Physical Laboratories of the University of Manchester. 
+ Communicated by the Authors. 
3H 2 


784 §.K. Runcorn, A. C. Benson and A. F. Moore on the Experimental 


Were the distribution of the Earth’s field exactly that of a dipole, the 
field components would vary radially according to an inverse cube law 
except in regions where the field vector is rotational. For a core theory 
such a region (of electric current flow) is inaccessibly deep within the 
Earth. The simplest interpretation of a distributed theory—though 
perhaps not the only one—is that the whole Earth is a region of virtual 
current flow, having magnetic but not electricaleffects. Runcorn (1948) and 
Chapman (1948 b) calculated, for the two groups of theories, the differences 
(AZ, JH) between the vertical and horizontal field components at .a 
small depth ‘“‘d’’ below the Earth’s surface (Z, and H, respectively) 
and the vertical and horizontal components at the surface (Z) and Hy 
respectively). They showed that though for a Core Theory, if R ise 
the radius of the Earth, 


AZ=Z,—Zoy=+ po i 
3:0.d 
and 4H=H,—H,)=+ ape He 


a Distributed Theory, assuming the density of the surface rocks is 2-7, 
gives 


ALT, Tee ae Vis 
AH=H,—H,=— = i Ee 


The essential problem of the experiment is to find conditions in which 
measured differences of the geomagnetic components between a point 
on the Earth’s surface and a point underground may be attributed to 
the main field and not to anomalies. For though the geomagnetic field 
is mainly a dipole field, the gradients of its components over small distances. 
at the surface.may be very different from the dipole gradients. - Thus. 
in England, while the surface gradients to be expected from the dipole 
field are of the order of 5 gammas per kilometre (1 gamma=10-* gauss), 
the observed gradients in some regions may be of the order of hundreds. 
of gammas per kilometre. These anomalous gradients seem to arise from 
surface features of the igneous basement rocks or from horizontal changes. 
in their magnetite content; the sedimentary layers above usually 
producing little effect. 

Therefore in regions where geological data indicates that the probable 
minimum depth at which extensive bodies of igneous rock are likely to 
occur is greater than 8000 feet, as in South Lancashire and Yorkshire, 
the gradients of the field components over distances of a few thousand 
feet, both along the Earth’s surface and vertically, should be sensibly 
linear. If this is so, then the effect of anomalies arising from the basement 
rocks on the measured values of 4Z and 4H can be estimated by a surface. 
survey and allowed for, if necessary. 
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For at any point where there is no current density (real or virtual) 
the gradients of the geomagnetic field along the 3 axes of a rectangular 
coordinate system are related by 4 equations, because 


curl H=0 and div. B=0. 
Such an equation is 


dH  dZ 
dz dx’ 


where z is a vertical distance and a a distance along a tangent to a line 
of longitude at the point. 

The measurement of the surface gradients may prove difficult as the 
values of the field at points on the surface are likely to be perturbed by 
the field of igneous boulders in the drift clay, thus some scatter of the 
points from a straight line has to be expected. Also the measurements 
may be made impossible owing to the vicinity of built up areas. Thus 
the preliminary measurements at the Parsonage Colliery at Leigh in 
Lancashire are of less interest because no adequate surface survey could 
be made in the region. 

It is difficult to find useful sites in mines as steel arches and rails are 
used in almost all passages and even if these are cleared some iron is 
usually buried. Passages also have to be a sufficient distance from the 
large concentrations of steel near the shafts, so that no gradient is 
detectable along the passage. Only then can the geomagnetic field in 
the passage be assumed to be free from man made disturbances. Only 
a very few of the cleared sites proved to be free of buried iron. 

Measurements of the vertical and horizontal components are made to 
an accuracy of 1 gamma by the quasi-absolute B.M.Z. (magnetometric 
zero balance) and Q.H.M. (quartz horizontal magnetometer) respectively 
which were developed by Dr. La Cour at the Danish Meteorological 
Institute. The temperature coefficients of the instruments must be 
known accurately and frequently checked in the laboratory, because their 
values are of the order of 10 gammas per degree centigrade and the 
difference in temperature between the surface and underground sites 
may be as much as 30° centigrade. For a determination of JH and 4Z 
measurements of H and Z are taken on the surface at a fixed point, the 
instruments are then taken to the underground site, the measurements 
repeated, and finally a further series is taken at the surface site. 
Throughout the experiment, instruments are left on the surface to record 
the temporal change in the geomagnetic field so that a correction for 
this may be applied to the results. The surface survey was made by 
H and Z variometers of the Schmidt balance type over a grid of 4:5 kms. 
around the surface site and from these measurements the surface gradients 
were computed and the value of the field at the surface site was corrected 
for any local anomaly. 
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The results of the measurements are shown in Tables I-V. ‘Table I- 
shows the position of the underground and surface sites where measure- 
ments were made, the depth of the former and the height above 
Ordnance Datum of the latter. The distances of these sites from the 
shafts and the steel materials and buildings associated with them are 
also given. Only in the case of the Nook surface site is there any 
interference, due to railway-sidings, and the correction for this is deter- 
mined from the surface survey of the area. This Nook surface site was. 
also used as the surface site for the Astley Green measurements, the 
necessary correction again being found from the surface survey of the area. 


TaBLeE I. 
Latitude and , Distance ,,, 
longitude and National Depth of eal from ee 
National grid grid under- OM.a8 shafts of h pis ¢ 
Colliery coordinates coordinates ground net ys under. **# fa 2. 
of under- of surface site Ba a ground ae ae 
ground site (feet) a °) site “ty 5 
site (reer) (yards) (varus) 
Hickleton 53° 32’ 42” N. 
Main 
(near PUT 36 We 
Doncaster, 
Yorks.). 44684056 44704065 2520 210 515 1200 
Nook 53° 29’ 50" N. 
(near 90997 
eich, 28a. iO" VV: 
Lancs.). 36924004 36914009 2930 100 560 410 
Astley 53° 28’ 40" N. 
Green 
(near 2° 267 48” W. 
Leigh, 
Lancs.). 37053980 36914009 2859 64 1900 410 


Table II. shows the results of the surface survey of the Lancashire 
and Yorkshire Areas. Values of H and Z on a kilometer grid covering a 
square with sides 4:5 km. from the surface base were used to compute, 
by a least squares method, the field gradients along each grid line. The 
gradients along the lines of latitude and longitude were then averaged. 
Figs. 1-4 show plots of H and Z against distances along lines of longitude 
(x) and along lines of latitude (y). The drawn straight line and the values 
given on the graph were derived by the method of least squares. The 
survey in South Lancashire was less satisfactory than that in Yorkshire, 
as small hamlets prevented measurements being made near certain of 
the grid points and the East-West gradient is rather large. However 
the corrections involved in 4H and AZ are small. 
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TABLE II. 
Surface Gradients of Geomagnetic Elements. 


ee 


Theoretical Experimental Theoretical eae care 
Colter value of value of value of Me ae ‘ 
eis gradient of Z gradient of Z_ gradient of H Peas 
/km jae /km gradient of H. 
y/km. y/km. y/km. Pict 
Variation with Latitude (Positive to North). 
Hickleton +5-0 +5-0-+0°3 —4-0 —3-4+0-5 
Nook and 
Astley Green +5-0 +4-0 —4-0 —6-0 
Variation with Longitude (Positive to East). 
Hickleton 0 +1-6+0-5 0 +1-2+0-6 
Nook and 
Astley Green 0 +3:2 0 —8-0 


a 
70 8-0 9:0 10:0 


—74 4 =! 1 
F 0: f 40 5:0 60 
10 20 30 Northing (hms) Ngan 


Hickleton Surface Survey. 


H 
Variation of H with latitude ae al —3-4+0-5y{km. 


e 
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Gammas 


Fig. 2. 


1 =e So 1 fe ee te 
FO 20 S004 Oe OG Onn CS 0 ees: Cee 00 
SOUTH Northing (kms) NoRTH 


Hickleton Surface Survey. 


Variation of z with latitude = =+5-0+0-3y}km. 


Fig. 3 
+10 
) 
~10 
SSS el! _—— = = 
West 4:5 65 a5 10-5 East (kms) se 


Hickleton Surface Survey. 


Variation of z with longitude a =+1-6+0-5y|km. 
Y 
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, The effects on the measurements of the induced magnetism of the 
sedimentary rocks was examined. From_ boreholes and shafts in 
Lancashire and Yorkshire representative collections of strata samples 
were obtained and the magnetic susceptibilities of the specimens were 
determined with an astatic magnetometer. Table III. gives the average 
susceptibilities for the different types of rock. 


Fig. 4. 
+10 


Gammas 
° 


a 2 
West 4:5 GS (kms) 85 10-5 East 


Hickleton Surface Survey. 


Variation of H with longitude oss =+1-2+0-6y| km. 


dy 
Tasxe III. 
(K is volume magnetic susceptibility.) 
Mudstones, 
shales and Sandstones Fireclays Tronstones 
Locality siltstones 
Number Number Number Number 


of Average . of Average of Average of Average 
samples K.10® samples K.10® samples K.10® samples K.10® 


Lancashire 40 23-5 21 14:5 a 19-4 14 56-7 
Yorkshire 60 23-4 a 11-7 6 14-1 4 69-7 


Calculations were made of the magnetic effects of the strata at Hickleton, 
Assuming that the various beds have uniform susceptibility and thickness. 
the vertical field components due to each bed at the surface and the 
underground observation points could be determined; if the horizontal 
extent of the bed were known. The effects could then be summed for 
the whole succession, the difference in the underground and surface values 
being a measure of the anomaly in 4Z. The horizontal extent of each 
bed is not known, but if it is assumed to be that for which its effect is a 
maximum, we get an anomaly which cannot possibly be exceeded but 
which is certain to be very greatly overestimated. A more realistic 
estimate is obtained by computing an average value of the anomalies 
which would be produced by beds of all sizes up to a certain arbitrary 
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maximum. The “maximum” and “average’’ anomalies in 4Z are 
shown in Table IV. The susceptibilities of the strata from the surface. 
down to 2700 ft. have been measured and the susceptibilities of the strata. 
beneath these down to the base of the Coal Measures at 4000 ft. have_ 
been estimated. The effect of all the sediments below the Coal Measures. 
is negligible. 


TaBeE IV. 
Vertical Component of Magnetic Fields due to Strata 
(gammas) 
Thickness of 
Strata over Maximum Effect Average Effect 
which the effect 
is summed Under- Under- 
(feet from surface) At ground  Differ- At ground __ Differ- 
surface 2367 ence surface 2367 ence 
feet feet 
(2725 + = 16 +30 utd 45 aL 12 
2728-4000 0 44 44 0 44 +4 
Total -F18 +6 
TABLE V. 


(Positive Values denote Increase in Component with Depth). 


Theoretical Experimental Theoretical Experiment al 


AZ AZ AH AH 4H 
Both theories Core Distributed 
yy ye ae Y te 
Hickleton +16 +11+83 +6 —13 +2242 
Nook +19 -Fli-+3 +8 —16 + 7+2 


Astley Green +18 +18+3 +7 =15 +10+3 


Table V. shows the expected change with depth of the horizontal and 
vertical components on both the distributed and core theories together- 
with the final experimental values. The measured values have been 
corrected so that they represent the difference between the components. 
in the mine and at a point immediately above on the surface. The values. 
of the geomagnetic components to be used at the surface are deduced 
from the surface survey and those underground from the survey of the 
passage in which the site is situated. It will be seen that the measure- 
mets are near to those predicted for a Core Theory and do not support. 
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a Distributed Theory. Only in the case of the 4H values at Hickleton: 
do the results deviate significantly (i.e. by twice the standard deviation), 
from the prediction of a Core Theory. 

Twenty sites in various mines have been magnetically surveyed but 
only six have been found to be free from buried iron. The result of the 
experiments in three of these sites has been described above. The 
experiment is proceeding in two more mines in the Doncaster area in 
order to obtain more 4H measurements and a site in a colliery near 
Sunderland, Co. Durham is also being used. The detailed experimental 
results will be published elsewhere. 


The experiments were done mainly from the Physics Dept. of 
Manchester University but are being continued from the Dept. of Geodesy 
and Geophysics of the University of Cambridge. We wish to thank 
Professor P. M.S. Blackett for his support and interest and Dr. Bruckshaw, 
the Danish Meteorological Institute, the Air Ministry, the Royal Greenwich 
Observatory, Dr. Veldkamp and Professor Whetton for the loan of 
instruments at various stages in the work. ‘We are also grateful to the 
National Coal Board and its officials for providing facilities for the: 
experiments and to many undergraduates of Manchester University for: 
their assistance in the measurements. 
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ABSTRACT. 


A systematic treatment is given of the energies and widths of domain 
boundaries in ferromagnetics ; it covers cubic crystals, with positive or 
negative magneto-crystalline anisotropy coefficients, and crystals with 
a single easy direction of magnetization. This extends the investigations 
of Bloch (1932), Landau and Lifshitz (1935), Lifshitz (1944) and Néel 
(1944) on cobalt and iron, and enables estimates to be made of the 
energies and widths of domain boundaries in nickel. The effect of 
magneto-elastic anisotropy energy is considered. Theoretical values of 
the energies and widths of the more important types of boundary in 
nickel, iron and cobalt are calculated. 


$1. INTRODUCTION. 


A FERROMAGNETIC material is, in general, divided into elementary 
domains in each of which the spontaneous magnetization vector, I), 
is constant, except possibly for small fluctuations in direction or intensity 
which may be neglected in the present calculations. The magnitude 
of I, is the same in all the domains and is equal, to a good approximation, 
to the saturation intensity of the material at the temperature considered ; 
the direction of I, varies from domain to domain. The transition region 
between contiguous domains is known as a boundary or wall and in it 
I, rotates from its direction in one domain to that in the other. 
Estimates of the energy associated with domain boundaries and of 
their effective thickness were first made by Bloch in 1932. Since then the 
treatment has been refined and further developed by Landau and 
Lifshitz (1935), Néel (1944) and Lifshitz (1944). Landau and Lifshitz 
considered a crystal, such as hexagonal cobalt, with a single easy 
direction of magnetization ; Néel and Lifshitz considered body-centred 
cubic iron, for which the easy directions lie along the cube edges, 
corresponding to a positive magneto-crystalline anisotropy coefficient. 
The aim of this paper is to obtain more general formulae from which 
results may be obtained for all the particular cases of interest, and to 
apply similar methods to the more difficult problem of nickel, which has 
not previously been considered. For nickel, which has a face-centred 


* Communicated by Professor E. C. Stoner, F.R.S. 


Energies and Widths of Domain Boundaries in Ferromagnetics — 793. 


structure, the easy directions lie along the cube diagonals, corresponding 
to a negative magneto-erystalline anisotropy coefficient. Most ferro. 
magnetic materials are either cubic, with positive or negative anisotropy 
coefficients, or hexagonal, and these are all covered by the general 
formule obtained for the relative energies and widths of different types 
of boundary. Estimates of absolute values, such as are made for nickel. 
iron and cobalt, and could be made similarly for materials for which the 
basic experimental data are available, are much more uncertain. 

In §2 the method of approach is indicated ; §3 deals with the energies 
and §4 with the widths of the domain boundaries neglecting the effect 
of magnetostriction. This is considered in $5, and in §6 the theoretical 
results for nickel, iron and cobalt are collected together. 

A survey of the background to the subject of this paper is given by 
Stoner (1948, 1950; see in particular 1950, §3). Reference may also. 
be made to a review article by Kittel (1949). 


§2. MeTHoD. 


In an unstrained ferromagnetic material in zero applied field each 
domain is magnetized along one of the crystallographically and 
energetically equivalent easiest directions. In domains in cobalt the. 
magnetization vector, I,, can only lie along the hexagonal axis and thus. 
only 180° boundaries (¢. e. boundaries in which I, rotates through 180°) 
can occur. In iron the domains can be magnetized along any one of the 
six equivalent cube edges, of type [100], and 90° and 180° boundaries. 
are possible. Domains in nickel can be magnetized along a cube diagonal, 
of type [111], giving rise to 70-53°, 109-47° and 180° boundaries. 

Néel (1944) emphasizes the importance of the condition that I,, the 
normal component of I,, should have the same value not only on each 
side of the boundary, but also in it. This corresponds to the absence of 
free poles (7.e. a zero value of div I) in the transition region. The 
presence of internal free poles would in general correspond to a free energy 
very much greater than that for some other arrangement satisfying the 
condition of zero divergence. It is convenient to specify the direction 
of magnetization at any point in the boundary by (9, 4), where @ is the 
angle between I, and the normal to the boundary, and ¢ the angle 
between the projection of I, on the boundary plane and some selected 
zero direction in this plane. In particular the initial and final directions 
of I, in the adjacent domains may be denoted by (4, ¢;) and (6, ¢5), 
constancy of I,, being equivalent to constancy of 0. 

The free energy per unit volume of the boundary may be expressed 
in the form 

eee De a eee Le 


where F,, is anisotropy and F’, exchange energy. The field energy is 
zero since the applied field and div I are both zero (¢f. Stoner 1950, 
equations (3.14), (3.15)). The energy and effective width of the domain 
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boundary are determined from the condition that F, integrated over the 
range in which ¢ changes from ¢, to ¢5, shall be a minimum. The 
anisotropy energy, F,,, may be expressed in the form 


| ond ied WS er ee ee ee) 5) 


Bha= Bet B'S tip] ee SL ae (2.3) 


where F, is magneto-crystalline and F,, magneto-elastic anisotropy 
energy; f,, f, and f,, are expressions for reduced energies adjusted, by 
subtracting appropriate constant terms, to make each zero for I, along 
(0, d,) and (0, d,); f, fp’ stand for the appropriate magneto-crystalline 
-and magneto-elastic anisotropy energy density coefficients respectively. 
The effect of the term F,,, will not be considered until $5. 


“OF, 


‘General expressions for boundary energy and width. 


It is convenient to consider the boundary to be in the yz plane with its 
normal along Ow and to take Oy as the zero direction for measurement of ¢. 
If 5% is the angle between the spins in adjacent layers 64—sin 654, and 
the exchange energy, which is proportional to 1—cos 65, may be written 
(Néel 1944, p. 7; cf. Stoner 1950, equation (3.1)) 


Feo" sin Oldgida ye. a. nal ien eres mer eee 


where « is proportional to the exchange energy for parallel spins, and 

a is a lattice constant. With the Weiss treatment of ferromagnetism « 

would be proportional to $N,,,1?, N,, being the molecular field coefficient. 
The energy, y, per unit area of the boundary is then given by 


y= [{ua? sin? O(dd/dx)?+Bf idx. . . . . . (2.5) 


It follows from a simple variational treatment that, for y to have a 
minimum value, 
Bi eat sin® 6(ddide)?,i wy ate Ly a tee) a ec) 


v.€. the exchange and anisotropy contributions to the energy are equal 
(see Néel 1944, pp. 6-8; cf. Becker and Doring 1939, p. 191; Stoner: 
1950, equations (3.1)—(3.6)). Equation (2.6) may be written 


(dajdd)=by sin Ofs*, 2 2°... . . (2.7) 
giving, by substitution in (2.5), 
“Po 
vivo=2 sind) fids, 2. 2. 8) 


where b, is a unit of length normal to the boundary, and y, a unit for the 
energy per unit area of the boundary (cf. Stoner 1950, equation (3.7)) 
given by 

bo=(aa?/B)t; yo=(aa2B)t. 2 2 2... (2.9) 
Equations (2.8) and (2.7) will be used to calculate domain boundary 
energies and widths. 
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Magneto-crystalline anisotropy energy. 


~ For a cubic crystal the magneto-crystalline energy, F,, apart from 
constant terms, may be represented, with sufficient approximation, by 
(ef. Stoner 1950, equation (2.1)) 


B= K,(ogag+ogat+atad). 2 2. 2. (2.10) 

where «1, %, «3 are the direction cosines of I, relative to the cube edges 

of the crystal. The reduced energy f, is a function of 6 and ¢ and may be 
written 


BFA, $)= (Fea, Penns PO OD (2.11) 


where (F,), 4 is the value of F, corresponding to I, along the direction 
(9, 6); (F.)min that corresponding to I, along an easy direction; and 
for cubic crystals, 8=|K,|. For iron K, is positive and thus 


Ea ames (te be (919) 


Fig. 2.1. 


Cubic crystal, showing easy directions in iron (edges) and in nickel (diagonals). 


For nickel K, is negative giving 
(F)mnin=Ka/3= — 8/3. Se eg IS) 


The form of the function f, can be determined for boundary planes of 
all possible orientations, and in the most general case a very complicated 
expression is obtained, but only certain boundaries satisty the condition 
imposed upon T,. All 90° boundaries in iron are similar to those in 
which I, rotates from [100] (i.e. OA in fig. 2.1) to [010] (OB); and [110] 
or [Ill] (i.e. any direction in the plane OCGF) are the only eetanenrs 
which make equal angles with these cube edges. In nickel all 109-47 
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boundaries are similar to those in which I, rotates from [111] (HD) to 
[111] (HE) and, again, only [11/] or [117] make equal angles with these cube 
diagonals. All 70-53° boundaries in nickel are similar to those in which I, 
rotates from [111] (HE) to [111] (HA); and [hk0] (i. e. any direction in the 
plane parallel to OAFB drawn through H) are the only directions making 
equal angles with these initial and final easy directions. All 180° 


boundaries in iron are similar to those formed by a rotation of I, from © 


[001] to [001] and in nickel from [111] to [111], the normals to the 
boundaries thus lying in the (001) and (111) planes respectively. Thus, 
for 90° boundaries in iron possible boundary planes are parallel to the 
{111} planes ({11l} being understood, in the usual crystallographic 
notation, to mean any plane of the type (11/)) ; for 109-47° boundaries in 
nickel to the {1i/} planes; and for 70-53° boundaries in nickel to the 
{hk0} planes. It is convenient to consider {001} to be {11/} planes with 
I—co; the parameters h, k and J can each take any value between 
—oand-+co. The {001}, {110} and {111} planes are particular examples 
of the. {11} planes; {100} and {110} of either the {117} or {hk0} planes. 
In these cases f, assumes particularly simple forms and it is then possible 
to derive analytical expressions for the boundary energies and widths. 
The normals to all 180° boundaries in iron and nickel lie in the {001} 
and {111} planes respectively. 

For hexagonal cobalt, the magneto-crystalline energy may be repre- 
sented, with sufficient approximation, by 


B26f.(6, d)=Ki sin’d, ee 2 Ae) 


where B=K3. The angle @ is equal to 90° and hence ¢ is the angle 
between I, and the hexagonal axis. This relation also holds, if the 
natural magneto-crystalline effect can be neglected, for materials of 
positive or negative isotropic magnetostriction under tension or 
compression respectively. 


Magneto-crystalline anisotropy energy of cubic crystals as a function of Oand ¢. 


If the normal to the boundary is the [111] direction it is possible to 
determine the direction cosines of I, with respect to the axes [111], [110], 
[viz]. A simple transformation then gives the direction cosines with 


respect to the [100], [010], [001] axes defined by the cube edges, and it 
follows from (2.10) that 


je ny [4(202-+1)+-4(2-+1)(2—4) sin? 9—(322+4)(2—4) sint 6 
4(1?—1){6-+ (I? 4) sin? 6} sin? 6 sin? d 
+ 4(/?-+-3)(?—1) sin* @ sin’ 6 
— 2°? 1 cos 6 sin 6 sin 4{2(/2—1) 


, —(5?+4) sin? 0-+-2(/?+-5) sin? @ sin? di]. . | (2.15), 


Le ee 
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This relation holds for {11/} boundary planes in any crystal whose 
magneto-crystalline energy may be represented by (2.10). In particular, 
I=0o, 0, 1 correspond to boundary planes whose normals are [001], 
{110], [111] respectively. The energy expression is now in a convenient 
form for use in calculating boundary energies and widths. 


§ 3. BouNDARY ENERGIES. 


The energy of a domain boundary when magneto-crystalline anisotropy 
and exchange energy alone are considered, given by replacing LaaVAl, 
in (2.8), is 


: pa 
y/yo=2 sin ay fo dd. So a eee ae) 
Nickel. #or a boundary plane in nickel whose normal is [111] the values 


of 6 and 43, ¢, corresponding to I, along the [111], [111] directions 
respectively, can be determined in terms of /:, 


6=cos—? 1/{3(2+2)}, OR 77 2 RE aS ee (3.2) 
d3=cos! {(?+2)/(2+3)}?, 0<dy<7/2 ; | 

(3.3) 
by=7— $s- 


The domain boundaries in nickel which are of particular importance 
are now considered, and the results are summarized in Table 3.1. 
A [001], 70-53° boundary is understood to mean a 70-53° boundary the 
normal to whose plane lies in the direction [001], about which I, rotates 
from the initial to the final easy direction. The values of 0, 4, and ¢, 
are obtained by substituting the appropriate value of / in (3.2) and (3.3), 
¢, and ¢, are found from 4, and ¢,, and f? is then determined by 
substitution in (2.15) and using (2.11) and (2.13). The relative energies, 
(y/vo), of the boundaries are then calculated using (3.1). 

The [001], 109-47° boundary is effectively two successive [001], 70°53° 
boundaries since I, must rotate through [111] (é.e. HC in fig. 2.1) in 
passing from [111] (HD) to [111] (HE). Its energy is, therefore, twice 
that of the [001], 70:-53° boundaries. It is impossible to have [111], 
70:53° boundaries since the angles between [111] and any two easy 
directions containing an angle of 70-53° (e.g. HD and HC in fig. 2.1) 
are not equal. 

180° Boundaries. All 180° boundaries are similar to those formed by 
a rotation of I, from [111] to [111]. A new system of axes Ow'y’z’ which 
has directions [111], [110], [112] (corresponding to /=1 in the [11/], [110], 
[712] axes of § 2) with respect to the system, Oxyz, formed by the cube 
edges of the crystal, may be set up. A particular boundary may be 
specified by w if I, rotates through OP, which lies in the y’z’ plane, and 
where w is the angle between OP and Oy’ (see fig. 3.1). The boundary 
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thus formed igs in the plane whose normal is [3* sin w#—cos w, 
—3! sin w— cosw,2cosw]. The angle w can have any value and, in 
particular, when w+30°=n60°, where n is any integer, I, rotates through 
an easy direction intermediate between [111] and [1IT]. It is thus 
convenient to write * 


Ae tcl 30S: all. wi Cae ee) 


and the normal to the boundary plane as [sin A—3? cos A, —2 sin A; 
sin A+ 3! cos A]. The angle A is that between the projections of I, and 
the easy direction [111] on the (111) plane. 


TABLE 3.1. 


Reduced energies of 70-53° and 109-47° domain boundaries in nickel. 

6, angle between I, and Oz, the normal to the boundary. 

41, $9, angles between projection of I, on yz plane and Oy, corresponding to 
initial and final easy directions respectively. 

f., reduced magneto-crystalline energy. 

y, boundary energy per unit area. 

Yo=(xa*B)* (see (2.9). 

n=cos—1(2/3)3(==35-26°) : 2n=angle between [111] and [I11]- 


Easy 
Type of directions; $4; eyes 
boundary initial, Y ie Be 22 f 7 33(y/yo) (y/vo) 
‘ final. : 
See lla 0 ke 
001], 70-53° [! ‘} ree Ll gagaie sere 
[001] rin 71 al. 332 sind cosd /3 0-5443 
ee) AGO) 0 eet <\ 
(O11 100-472, eee alee \ 92982 sin bc 95/2 aeies 
[O01] / 74. 37 —1 e 3°? sin d cosh /3 1-0887 
aoe 111] ‘A hs 
1107, 70-532 (LUE 4g ig ante ae ae : 
Lar htee [LE iees puma OT cael renee rnti RTE Waa fae i 
[110], 109-47° LULL}, 37 is 32(3 sin? d—1) 24!+2/2—y 1-3680 
[111]. SPE cies a 
Pena s\n ee 
[111], 109-47 rs ‘ 2 ete 272(1—gin 34)! 2152/34 -1.2903 


The magneto-crystalline energy may be represented by 


F,=(K,/12){4 cost 6+3 sint 6+ 2°” cos ¢ sin? d sin «(3-4 sin? w)}, (3.5) 


* This angle, following the usage of Néel (1944), is denoted by A. It is not 
to be confused with a magnetostriction coefficient, which is always written 
in this paper with a directional suffix, e. g., Arg9- : 
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by putting J=1 and replacing 6 by ¢ and ¢ by w in (2.15). The 
magnetization vector, I), must rotate in the plane of the boundary and 
sin @ in (3.1) is, therefore, equal to unity. From (2.11), (2.13), (3.1) and 
(3.5) it follows that 


ylyo=3* F (1+-p sin ¢ cos +7 cos ¢)i sin¢dd¢, . . (3.6) 
“0 
where 
p=—2°" sin 3 | - 7 
eee yi, 3.7) 
=2% cos 3A | 


Digests 


Coordinate axes used in describing 180° boundaries in nickel. 
OP, direction in y’z’ plane through which I, rotates. 
w, angle between OP and Oy’. 
¢, angle between I, and Oz’. 


It is possible, without any difficulty, to evaluate the integral analytically 
when p=+2°? since (sin? ¢-+2°” sin ¢ cos +8 cos* ¢) is a perfect square. 
Similarly when p=0 the integral can be evaluated analytically in terms 
of an inverse trigonometric hyperbolic function. For a general value of p 
the integration is slightly more troublesome. 

312 
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It is possible to make the mixed term of (3.6) equal to zero by writing 
Bei ee a 


where y is a new variable and « is a function of p only, and choosing « 
such that : 
tan 2«=p/7. PMN eenn neh ee) eS) 


Then 
y/yp= 3 3(7 sec2x)? {| 


1 —K ; 
y2—gin? y)? sin xk cos yd 
x X&X 


135 


pu—k 


a (r?—sin? y)? cos « sin xdx} ; #(3-10) 


=—K 
The first integral is evaluated by substituting w=(sin x)/r and the second 
by substituting v—(cos x)/s and, after some manipulation, the energy 
of a boundary whose normal is [sin A—3? cos A, —2 sin A, sin A+-3* cos A] 
may be written A 


ylyp= 3-H 22+(7 sec 2x)'(r3t sin-1 t+s%p sinh-! p)}, . . (3.11) 
where 
1472=9 cos 2k+7 | t=(sin x)/r i 
@ Gees ia |e (3.12) 
14s?=9 cos 2x—7 | p=(cos k)/s. | 


When A=0°, 120°, 240°, p27; 760°, 180°, 300°, p= — 2°, and these 
six particular cases correspond to boundary planes of the (110) type 
and to a rotation of I, through an easy direction intermediate between 
[111] and [III]. Thus the energy of these boundaries is equal to the 
sum of the energies of [110], 70-53° and [110], 109-47° boundaries 
(cf. Table 3.1) and by substitution in (3.11) 


vive ot (2? + a/2—— 29) 1-202. 5) a (813) 


‘This is the minimum value for the energy of 180° boundaries in nickel. 
When A=30°, p=0, the normal to the boundary is [112] and the energy 
as given by (3.11) is a maximum, 

ylyp=3-4{292-+47-4 sinh-! (74)} =2-0040. . . . (3.14) 


Iron. For a boundary plane whose normal is [11/] the values of 
6 and ¢5, ¢4, corresponding to I, along the [100], [010] directions are 
given by 


'é@=cos1(/?+2)-4 Oc Oe, hae team salt 
d3—cos 1— {(0?42)/2(2+1)4, 7<43<3n/2; ] 
gy=cos t{(24+2)/2(1?+1)}* , 32/2<¢,<2r. | 


‘The domain boundary energies may be calculated in a way similar to 
that used for nickel, using (3.15), (3.16) and (2.12) instead of (3.2), (3.3) 
and (2.13). The results are summarized in Table 3.2. 


(3.16) 
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MABE #323 


Reduced energies of 90° domain boundaries in iron. 
In all cases the initial and final easy directions are [100] and [010] 
respectively. (See Table 3.1 for symbols.) 


Type of $y, ; 

boundary” ie ai rIYo ro 
° —3n/4, ‘ fs ae 

[001], 90° a/2 ==5/4 2(1—2 sin? ¢) 1 1-0000 

[110], 90° =/4 ~G psingis—3 sin? p)*? 1+2-323-45 sinh-1(3/5)? 11-7274. 
O__ fs —57/6, 93 2-3/2 * 4 95/23 { > 

[111], 90°x/2—y —alG. (2 3-32(1-+-sin3¢)* 25/3 1-1852 


Fig. 3.2. 


Coordinate axes used in describing 180° boundaries in iron, 
OP, direction in «xy (y’z’) plane through which I, rotates. 
w, angle between OP and Oy’. 
A, angle between OP and Oy. 
dé, angle between I, and Oz (O2’). 
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180° Boundaries. A rotation of I, from [001] to [001] through OP, 
which is in the zy plane (see fig. 3.2) is considered by substituting /= co 
and replacing 6 by ¢ and ¢ by w in (2.15). By writing 


N=04-45° Se ae 


(cf. equation (3.4)) the magneto-crystalline anisotropy energy may be 
expressed by 
f=sin? d—psintd, .. 7 5) Soa) 
where 
p=(7—Ccos 4w)/8, 


(3.19) 
=(7-+cos 4A)/8 (cf. Néel 1944, p. 11). 


The angle A is that between [010] and the projection of I, on the (001) 
plane. 

The normal to the boundary specified by 4 is [—cos A, —sin A, 0]. By 
substitution in (3.1) it follows that 


ylyo=2{1+ (1—p)p sinh“! (u/1—p)*}. . . (3.20) 


When A=n7/2, where n is any integer, ~=1, and the boundary can be 
considered to be the sum of two [001], 90° type boundaries. Its energy 
as calculated from (3.20) is twice that of a [001], 90° boundary and is 
given by 


Vga - Se Beare (S20) 


This is the minimum value for the energy of 180° boundaries in iron. 
If A=n7/2+45°, corresponding to »=3/4, the energy is a maximum and 
equal to 


yive=o-oe sinh (3t)==2:'7603 2 ea (322) 


for boundaries in the {110} planes. 

The relative energies of domain boundaries in iron have already been 
calculated by Néel (1944), whose method is essentially the same as that 
employed above. 


Cobalt. Only 180° boundaries occur in hexagonal cobalt and F, is 
zero along the easy directions. Thus, from (2.14), 


Jo=sine gd oe ey 
giving, from (3.1), writing sin 0=1, 
Y/Yo=?2 sin p dp 


a RR Cs ee me hey 


The energies of all boundaries in cobalt are the same. 
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§4. Bounpary WIpTus. 


If exchange and magneto-crystalline anisotropy energy alone are 
considered the distance x, measured normal to the boundary, can be 
obtained as a function of ¢, the angle between the projection of I, on the 
boundary plane and the zero direction, by replacing f, by f, in 
equation (2.7), giving 


(a Co \==Demin 7 a key wean ut (4.1) 


This section will be devoted to a consideration of the variation of x with co) 
for the various boundaries considered in §3, and of the resultant 
‘boundary widths ”’. 


Fig. 4.1. 


Variation with distance normal to the boundary plane of the angle, ¢’, through 
which I, has rotated, for [110], 70-53° and [110], 109-47° domain boundaries 
in nickel. 

The broken line is an asymptote common to the two curves. €=2/bo. 
x, distance normal to boundary. 


For ¢>¢4, or ¢5, f,>0,|xz|—, and if the boundary is to satisfy the 
minimum energy condition it has no sharply defined width, as is shown 
in figs. 4.1, 4.2 and 4.3. It will be convenient to introduce a non- 
dimensional parameter (cf. Stoner 1950, equation (3.8)), 

Cath. |. eee tre) (4. ) 
for measurement of length normal to the boundary. All ¢, € curves for 
70-53° and 109-47° nickel, 90° iron and 180° cobalt boundaries are 


symmetrical about their points of inflexion which lie at d=(4,+¢,)/2. 
For such boundaries the effective width, b, may be conveniently taken as 


(b/b)) =(b2—$1)(4E/44) nfexchw eet Se (4.3) 
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where the notation is self-explanatory. All ¢, € curves for 180° boundaries 
in nickel and iron are of the form indicated in figs. 4.2 and 4.3 respectively, 
and the effective width, b, may be taken as 


(b[bp) == (Eo—£.) + bq(dEldb)g-+(m—b_)(dE |b)...» (4-4) 


Fig. 4.2. 


Variation with € of ¢ for 180° boundaries in nickel. The numbers on the 
curves indicate the values of X; the curve marked 0° (with m.-s.) is that for 
A=0° when magnetostriction is taken into account (see § 5). 
£=2/b). «x, distance normal to boundary. 
¢, angle between I, and [111]. 


Fig. 4.3. 


0 aA! 
C0. 8. 6 9 4= RanO qe eeciemrG ane Clear 


Variation with € of ¢ for 180° boundaries in iron. The numbers on the curves 
indicate the values of A; the curve marked 0° (with m.-s.) is that for A=0° 
when magnetostriction is taken into account (see $ 5). 

€=a/by. x, distance normal to boundary. ¢, angle between I, and [001]. 


Se ae 
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where (¢,, &,), (¢., €,) (see figs. 4.2 and 4.3) are the first and third points. 
of inflexion. It will_be-noted that all 180° boundaries in iron, and 
180° boundaries in nickel when A=n60°+-30° (where A is defined as the 
angle between the projections of I, and the easy direction [111] on the 
(111) plane (see (3.4)), are symmetrical about ¢=7/2. The expressions (4.1). 
and (4.2) are used to find & as a function of 4, the constant of integration 
being omitted, and (4.3), (4.4) are used to determine the boundary 
widths. The results for 70:53° and 109-47° boundaries in nickel are 
summarized in Table 4.1 and those for 90° boundaries in iron in Table 4.2. 
The values of 0, J, and ¢, and f,! are given in Tables 3.1 and 3.2 and are- 
not repeated here. 
TABLE 4.1. 
Reduced widths of 70-53° and 109-47° domain boundaries in nickel. 
=clbe 
E oe normal to boundary. 
b, effective boundary width. 


by, unit for measurement of x and b; by=(aa?/B)* (see (2.9)). 
n=cos1(2/3)?(=35-26°). 


Type of boundary (2/3)2é (2/3)#(6/b,) (b/bo) 
[O01], 70-53° In | tan ¢ | 7 3°8476 
[O01], 109-47° oe) ee) 
[110], 70-53° 2 tanh—! (2? tan ¢) 2/2 4.2642 
[110], 109-47° —2 coth1 (23 tan ¢) 22(a7—2n) 3°3093 
[111], 109-47° In | tan (7/8—3¢/4) | 7 3-8476 

TABLE 4.2. 


Reduced widths of 90° domain boundaries in iron. 
(See Table 4.1 for symbols.) 


Type of boundary E (b/b5) (b/d) 
[001], 90° In| tan (6+7/A4) | 7 3-1416: 
[110], 90° —sinh-1 { (8/5)? cot 4} (8/5)4ar 3:9738: 
[111], 90° —In| tan (37/8—3¢/4) | 7 3-1416 


It may be noted, in connection with the expressions for € in the Tables, 
that the inverse hyperbolic functions may be expressed as logarithmic 
functions as follows : 


sinh~! 4=In{x-+ (x?-+-1)?} ; 
2 tanh-! x=In{(1-+2x)/(1—x)}, 2 
2 coth-1 w=In{(x+1)/(x—1)}, # 
2 tanh! (tan ¢)=/n | tan (6+7/4) | - 
Values of these functions are tabulated by Comrie (1949). , 
If the magnetostriction term, /,,, 1s not introduced, the [001], ee 
boundary in nickel can only be regarded as two successive 70°53 


boundaries, and as these would necessarily have an infinite distance 
between them, the resultant boundary width would be infinite. This. 


case will be considered in § 5. 


(4.5) 


1s, 
tt 
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, In fig. 4.1 the 4, € curves for [110], 70-53° and [110], 109-47° boundaries 
in nickel have been plotted. All 4, € curves for 70-53° and 109-47° nickel 
and 90° iron boundaries are similar to these. 


180° Boundaries in nickel. The parameter € is expressed as a function 
of 4 by using (4.1), (4.2), (3.5), (2.11) and (2.13) and may be written 

&=3'2f[(1+ p cot +8 cot? d)-? cosec? ddd eating MELA: a4) 

= — (3/2)! sinh~1 {(16 cot ¢+p)/(32—p?)*} l 


(4.7) 
=— (3/2)! sinh~! {cosec 3A(23? cot 6+ cos 3A)}. | 


The ¢, € curves corresponding to A=1° and 30°, where A is the angle 
between the projections of I, and the easy direction [111] on the (111) 
plane, are plotted in fig. 4.2. The curve for A=30° is a limiting case; and 
as X approaches 0° the central portion broadens and becomes more 
nearly parallel to the € axis until, when A=0°, the boundary width is 
infinite and the wall can be regarded as two infinitely distant walls of 
the form shown in fig. 4.1. When 


A=’ 0°,,(5/6,)—=0o 3 } 

A= 1°, (b/bp) = 15-4029 ; 

A=30°, (b/by) =2(3/2)#{sinh-? (6#)-+ (7/2)# sin? (2/72)} 
=7-9161. 


(4.8) 


180° Boundaries in iron. For iron € may be written, by using (4.1), 
(4.2) and (3.18), in the form 


é== [sino sind) x" date ne 8) a eee oe eo) 
=—sinh—! {(cot ¢)/(l—y)*}=—sinh— (2 cosec 2A cot d). (4.10) 
Curves, similar to those for 180° nickel boundaries, are plotted in fig. 4.3, 
that for A=45° (A being the angle between the projection of I, on the 
(001) plane and the easy direction [010]), corresponding to a normal in 
the direction [110], being a limiting case. In the other limiting case 


A=0° and the boundary can be considered to be two infinitely distant 
boundaries of the [001], 90° type. When 


A==07; (b/by) = © ; 


A=1°, (b/by)=12-6258 ; 
A=22:5°, | (b/bo)=2{sinh-? (6) + (7/2) sin-} (2/7#)} 
= 64635 ; Sree 


A=45°, — (b/by) =2{sinh-} (2#)-++3# sin-1 (2/3)#} 
=5:6018. | 
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Boundaries in cobalt. For cobalt € may be written, by using (4.1), 
(4.2) and (3.23), in the form 


Ge] COSCO OU et enw  . e412) 
Nie COE alge | ea ae er a) 


This curve is similar to those for 70-53° and 109-47° nickel and 90° iron 
boundaries. The reduced boundary width, 6/b,, is equal to 7 and is the 
same for all the boundaries. 


§5. Errect oF MAGNETOSTRICTION. 


The introduction of a magneto-elastic anisotropy energy term, F,,, 
has important consequences. This section will be devoted to a considera- 
tion of its effect on 180° boundaries, and particularly to those which have 
been considered in §§3 and 4 to be the sums of two infinitely distant 
boundaries. They no longer have infinite “ boundary widths ”’ since the 
effective result of the introduction of this term is an attraction between 
neighbouring 109-47° and 70-53° nickel boundaries and between adjacent 
90° boundaries in iron (see Stoner 1950, § 3 (i); Kittel 1949, pp. 563-4). 


The strain tensor and the energy term B'f,,. The magneto-elastic energy, 
F,,, of a cubic crystal is the sum of magnetostrictive and elastic parts 
and may be written (Becker and Doéring 1939, p. 146) 


Bin={—3CgAy99(atA 1 +0GAoe+ a3 Acq) 
—3C A431 12 («1% 2A 12+ %9%3A53+%3%,A5;)} 
+{$C,(Ayi+Ag2t+Ags)?+C2(At +Aget Ags) 
+20,(A7,+A5,+A51)}, 


where the A;, are components of the strain tensor, and Ajo9, Aq are 
magnetostriction coefficients. Comparison of (5.1) with Bitter (1937, 
p. 172, equation (138)) shows that the constants C,, C,, C, are given 
in terms of the elastic coefficients C11, C12, C4, by 


Ci=Cy2 ; Cy= (C41 —649)/2 5 Cy=Cy,. - . |. (5-2) 


The stable configuration of the crystal may be determined by minimizing 
F,,, with respect to the A,,. It follows that 


A jp=(3/2)Aq 11% i es A j= (3/2)A yoo{0G —Cy/(83C,+2C2)}, - (5.3) 

tk; 
from which it is evident that the components, A;;, of the strain tensor 
are the same for domains magnetized in opposite directions. It may be 


shown further (cf. Lifshitz 1944, pp. 338-9) that they remain constant 
throughout the 180° boundary, 7. e. the deformation of the material of 
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the boundary is effectively the same as that in the adjacent domains. 
The A,,, are found by substitution in (5.3) of the direction cosines corre- 
sponding to the easy directions in the domains. 

The excess magneto-elastic energy of the boundary, f’f,,, with respect 
to the domains can now be calculated by substitution in (5.1) of the 
appropriate values of [a,«%3], corresponding to the directions of I, in 
the boundary and in the domains. (See Lifshitz 1944, p. 339; Kittel 
1949, p. 564.) 


Nickel. The directions of magnetization in adjacent domains are 
[111] and [111] (cf. §§ 3 and 4) and the strain tensor has components 


Ajg=Agg—=Agy=Aq44/2 ; Ay y=Age—=Agg—=C oA 199/(83C1+2C,). (5.4) 


For a rotation of I, through OP (see fig. 3.1) the values of «,, « and a 


in the boundary are 
a= 6-#{2! cos 6—(3? sin A+ cos A) sin ¢$}, 


%& = 6-#{2? cos +2 cos A sin d}, ea ee 8 ashi 
a,=6-}{2? cos ¢4-(3* sin A—cos A) sin J}, | 


giving, for the excess magneto-elastic energy, an expression, 


Bfin=(9/2)CgAqas? sin? 6, (56) 


which is independent of the angle A. The reduced anisotropy energy 
may be written, from (2.3), (2.11), (2.13), (3.5), (3.7) and (5.6), in the form 


fa=((L+7)/12}{1+ p’ sin ¢ cos d+o cos? 6} sin? 6,  . (5.7). 
where 
T= (540303 41")/ | K,|; ep =p/(1-+7) ; o=7/(l-+-7). .' (5.8) 


Substitution in (2.8), with sin 6 equal to unity, and comparison with 
(3.6) gives 


yiyo={(l+7)/3}* Ik (1+p’ sin ¢ cos ¢+0 cos? $)? sin 6 dd, | 
0 
={(1-+7)/3} {(6+1)!+(o sec 2«)#(r3t, sin-! t, Pe eel i!) 
+SiDy sinh? p,)}, 
where 
2ori=(o+2) cos 2x-+o 


208?—=(o+2) cos 2x—o 


t3=(sin k)/r, | 


> s=ri—l; j (5.10), 


p,=(cos K)/8y, 


and tan 2xk=p’/o=p/7. 


Widths of Domain Boundaries in Ferromagnetics 809 


The elastic coefficients of single nickel crystals have been measured by © 
Bozorth et al. (1949). -Values of the magnetostriction coefficients are 
given by Becker and Déoring (1939, pp. 277-80), and of the magneto- 


crystalline anisotropy coefficient, K,, by Stoner (1950, § 2). The constants 
required are 


6, ;=2°50 70 ers cm, A109 = — 46 X 10-8, 
Cy.>= 1-60 X 10” erg cm.-3, Aqq=— 25 x 10-6, 
£,,—1-185 x 10% erg em.-3, K,=—0-50 x 105 erg cm.-3, | 
: (5.11) 
and the value of 7 calculated from them is 0-80. 
When i 
p=t22, y/yo=3-2?[(8-+7)!24 {(9+7)/9} sin-} (9-+-7)-¥2 
+ (27/9) sinh- (8/7)3], 
aS See. ie | i air OM T0515) 
p=0, —— y/yo={(1+-7)/3}**{(o+1)¥*-+074? sinh (01)}, 
ee ee Take oe hee ee 4 eee e139) 


These energies are greater than those calculated in §3 where the effect 
of magnetostriction is neglected. 

The parameter € may now be expressed as a function of ¢, using (5.7), 
by 


&=—2{3/(8+7)}? sinh-1 [{2(8+7) cot b+ p}/{4(8+7)(L+7)—p?}1?]. 


(5.14) 

When 
p=+25?, b/by=791; | Bae 
p=0, b/bp=4:45. j , 


The boundary width no longer approaches infinity as p approaches 4235/2 
and the introduction of a magneto-elastic energy term produces effectively 
an attraction between [110], 70-53° and [110], 109-47° boundaries. 

Iron. The directions of I, in adjacent domains are [001] and [001] 
(cf. §§ 3 and 4) and the strain tensor has components (cf. Lifshitz 1944, 
pp. 338-9; Kittel 1949, p. 564) 


Ayp=Ag3=Agi=0 5 Ayy=Age= 3A 4 99C1/2(8C,+2C2) > 
Agg=3A199(C1+Cs)/(83C1 +203). - ees). 10) 


For a POURHOn through OP (see fig. 3.2) the direction of I, in the boundary 
is [— sin A sin 4, cos A sin d, cos ¢] and, hence, 


B'fim=(9!2)C2A 907 Sin? $, Peter we ess (Onk 1) 


810 B. A. Lilley on the Energies and 


which is independent of the angle A. The reduced anisotropy energy 
may be written in the form 


ji (lea) (sin? p—p sin ), See reo. 8S (5.18) 
where 
T= (ICaAteg-) 2a he LE) pe en On LoDy 
It follows that 
viyo= 2 +7’ )HL+ (1p!) sinh? {y'/(1—w’}4, (5.20) 


and that. 
€é=—(1+7’)-} sinh! {(1—p’)-* cot ¢} ; 
blbp=2(1-+-7’)-? [sinh-? {(2u/—1)/(1—p’)}# +2’? sin-? (2u’)-4]. 


The constants ¢,,, Cy2, C4q are given by Kittel (1949, p. 556), Arog, Arai 
by Becker and Doring (1939, p. 280), and K, by Stoner (1950, § 2) as 
follows : 


5.21) 
ie 


C,,= 2/4) x10 Mere ems) Agog = 19°5.410-4, 
Cyo= 1-46 10Mergem-=, A, ,, = — 18-8 x 10-6, . (5.22) 
C112 10 Mere cms. eda 1 0° ere icin: 


The value of C, obtained from these is 0-48x10" erg cm-* giving a 
value of 


7 =0-0017. 
When 
A=0,_ ‘} 
y/yo=2{(1-+7')#+-7’ sinh-1 (7’-*)} 
eM) 
a2O2 70) gat en Mors, e-paper Ba ee eee 
b/bo —2(1-+7’)-#[sinh-! {(l—7’)/7’}4 
+2(L-p7')# sin“? ((1+7')/2}4] 
=10°87 . . . - . . fe * 4 - ‘ (5.24) 
A=457% 


—— | neater) (-b4e89 sinh-1 {3/(1+47')}4] 


2:77, eae Tee, 
bby =2(1+-7')-# [sinh-? {(2—47')/(1+ 47’) } 
+ {3#/(1+-7’)} sin-t {2(1-+-7')/3}] 

650/04 2) 4 a: Ee een 
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Thus the inclusion of a term representing the magneto-elastic energy of 
the crystal has little effect on the value of the domain boundary energy 
(e.g. 2-02 as compared with 2-00 for a [001], 180° boundary ; for other 
values see Table 6.1), but results, as in the case of nickel, in an important 
modification of the results obtained, in § 4, for the boundary widths. 

In this section the numerical values of 7 and +’ are given to two: 
significant figures ; those of y/y) and 6/by are given to three figures. The 
accuracy of these values is, however, determined by the accuracy with 
which the constants concerned (given in (5.11), (5.22)) have been 
measured. 

Cobalt. The effect of magnetostriction in cobalt is not important since- 
only 180° boundaries are possible and none of these has an infinite width 
even if magneto-elastic energy is neglected. Moreover, the introduction 
of this further term would have very little effect on the boundary energy 
since the magnetostriction coefficients of cobalt are smaller than those 
for iron, and the magneto-crystalline anisotropy coefficient, 8, is much. 
greater (by a factor of about 10) than that for iron. 


§6. COLLECTED RESULTS AND CONCLUSION. 


In conclusion, it is appropriate to tabulate the results obtained. This 
is done in Tables 6.1 and 6.2, in which ‘“‘ with m.-s.”’ is to be understood 
to mean “with magnetostriction taken into account’. In cobalt the 
hexagonal axis lies in the domain boundaries. 


ABLE GL le 


Reduced energies, y/y), of domain boundaries in ferromagnetics. 
y, boundary energy per unit area. 
For estimated values of yp see text. 


Type of boundary Direction of normal to boundary 
[O01] [110] [111] [112] 
70-53° 0-5443 0-4611 ; me = 
109-47° 1-0887 1-3680 +290: 
Nis 180° = 1-8292 a 2.0040 
180° (with m.-s.) < 2-19 = 2.27 
90° 1-0000 1-7274 1-1852 
Fe< 180° 2-0000 2-7603 — 
180° (with m.-s.) 2-02 2-77 — 
Co 180° 4-0000 


For 180° boundaries in both nickel and iron the angle A is that between 
the projections of I, and an easy direction on a plane whose normal is 
the initial direction of I, (see §3). In Table 6.2 the values of 6/b, for: 
\—1° are included in brackets below the infinite values for A=0°. 
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TABLE 6.2. 
Reduced widths, 6/b), of domain boundaries in ferromagnetics. 6, effective 
boundary width. 
For estimated values of by see text. 
Type of boundary Direction of normal to: boundary 
[001] [110] [111] [112] 
70-53° 3-8476 4-2642 — 
109-47° oe) 3-3093 3°8476 
Ni< 180° — oe) — 79161 
(15-4029) 
180° (with m.-s.) a 7-91 . — 4-45 
90° 3°1416 3:9738 3-1416 
Fe2 180° oe) 5-6018 
(12-6258) 
180° (with m.-s.) 10-87 5-59 
Co 180° 3-1416 


The effect of introducing a magneto-elastic energy term is most 
important in nickel, less important in iron, and of least importance in 
cobalt. Values of y, and by are given by Stoner (1950, table 3.2), who 
also discusses the uncertainty of these estimates due to uncertainties 
in the detailed treatment of exchange interaction effects. For nickel, 
iron and cobalt y ), measured in ergcm™, is 0-14, 0-62 and 2-05; 
bo, measured in 10-* cm, is 260, 130 and 50. (Néel (1944) obtains for iron 
values of 0-7 and 165 for y, and by respectively ; Kittel (1949), using the 
method employed by Lifshitz (1944), in which the exchange energy is 
calculated from the coefficient of T?? in Bloch’s law of approach to 
saturation at low temperatures, obtains corresponding values of 0-9 
and 230.) Thus the energies, y erg cm~*, of 180° boundaries in nickel, 
iron and cobalt (taking A—0° for nickel and iron, and without neglecting 
magnetostriction) are 0-31, 1:25 and 8-20; the widths, b, in 10-8 cm 
are 2060, 1410 and 160 respectively. These comparisons have a bearing 
on the character and relative ease of production of powder patterns and 
on the sizes of domains in the three ferromagnetics (see Stoner 1950). 
For nickel, iron and cobalt, the boundary widths, at room temperature, 
are, in general, small compared with the linear dimensions of the domains. 
If the boundary widths are of the same order of magnitude as the domain 
size the ferromagnetic can no longer be regarded as being magnetized 
along some easy direction over most of its volume. This effect may be 
important in many alloys, and also in the ferromagnetic elements 
themselves at higher temperatures. 

I am indebted to Professor E. C. Stoner, F.R.S., for suggesting the 
problem and for his guidance and constant encouragement; to Dr. P. 
Rhodes for helpful discussions and constructive criticism; and to the 
Department of Scientific and Industrial Research for a grant. 
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ABSTRACT. 


Fundamental criticism is made of two recent papers by Meltzer on 
electrical noise on general and particular grounds. Correction of the 
errors in the analysis also provides an opportunity for a discussion of the 
fundamental relation between a stochastic equation and the corresponding 
frequency spectrum of a fluctuating variable. An example taken from 
the field of classical Brownian movement clarifies the concept of the 
observable velocity of a particle. 


Two recent papers by Meltzer (1949, 1950) in this Journal on electrical 
current fluctuations (noise) present much cause for criticism. The first, 
(A), implies that a proof has been given of a fundamental physical 
identity between thermal and shot noise. But in brief the paper simply 
reiterates the well-known fact that any purely random phenomenon, 
which may be idealized as a Markoff process, obeys a certain rather 
simple statistical law. Membership of the same statistical class, however, 
is remote from an essential physical identity of two phenomena, and the 
suggestion that a serious proof of this has been given is at least rather 
misleading. The second paper, (B), however, takes up the problem of 
current noise in metallic resistors and uses analysis which is invalid 
in fundamental, but elementary, respects; the ensuing results and 
comparison with the experimental work of Bernamont are consequently 
fallacious. 

We first consider the specific errors of the second paper and then take 
the opportunity of a more general discussion (relevant to both papers) 
of the statistics and frequency spectrum of fluctuations. An example 
taken from the field of classical Brownian movement clarifies the concept 
of the observable velocity of a particle. 

It is simplest in considering (B) to present afresh in correct form 
Meltzer’s analysis, adhering as far as possible to his notation. Consider 
some cross-section of a conductor and let us observe the charge crossing 
this during a set of time intervals of (arbitrary) length, 7 +. Let n be the 
total number of charges e;, of magnitude e, crossing in any particular 
member of the set. ‘This number will of course not be constant, because 


* Communicated by the Author. 
+ We prefer to use 7 in place of Meltzer’s ¢. 
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of the very presence of fluctuations, but will vary from interval. Let Ny 
be the number crossing in the “ positive” direction and M. in the 
“negative ” direction, and let the respective probabilities be p and q. 
Thus denoting the mean values of the three statistical variables 
n, Mm and ny by N, N, and N, we have: 


. N,=pN ; N,=qN 
and the mean current, I, is given by: 
eN 
I=(p—q) —. Sanremo © ghee. come L) 


Then, following Meltzer, we set : 


Q= ‘eh €;=(n4—N,)e (2) 
Q= (i, —i,)e=Ne(p—g)=Ir 
and Q?=2e?+ L ee, (3) 
i ij 
Q?==7 6? +n ce? +n4(n,— 1)e?+-n,(ng—1)e? —2nn,e? (4) 
=(N,+N,)e?+ (nN, )e?+ (n2-N, Je? —2nynge” 
Now let 
m,=N,+6,; m2=N,+5, 
. W=N248? ; n2=N3+83 ; mm, =NiN, +5152 
and Q?=Ne?2+ (N, —N,)?e2-+ (52-++ 52 25,8, (Ni +N,))e2 . (5a) 


=Ne?+N?(p g)re2-+ (A? INVES Oe Sk es OO) 


if we set 4=6,—5,. 

This should be compared with Meltzer (B), equation (5). Meltzer has 

throughout made no distinction between average and mean square 
values, etc., of the relevant statistical variables, which appears a rather 
serious error in a paper intended to deal with fluctuations. 
The achievement of further real progress depends now very naturally 
- on the assumption of some physical model for the system considered. 
The only assumption of a general nature that we can essay here is that of 
an overall equilibrium distribution : 


yap den=d\h Ps ix. 50 as ce Mm OLED) 
or more specifically to assume : 
SN, ; S=N,; 58-0 ©... . , (68) 
both yielding : 
MPS NezENAp— gyre? cw we wt 


2 
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and using (1): 
@aNet Pe 9 Sal eee 


Using finally Einstein’s classic result (Einstein (1906)) that for a 
resistor, R, in thermal equilibrium at temperature T, the effective 
fluctuation in charge (Ze? after equation (3)) is given by (2k7'/R)r, we 
may write : 


This is of course just Meltzer’s starting point, (B) equation (1), from 
which no amount of algebra by itself can possibly derive any essentially 
new expression for the fluctuation involved. By the very fact of writing 
down this equation, Meltzer has implicitly assumed a complete knowledge 
of the statistics of the fluctuation in the presence of the term in 7, 
representative of a purely random component, and the term in 7”, repre- 
sentative of a pure “ d.c.”’ contribution. The simple addition of the two 
terms embodies immediately the assumption that these are entirely 
independent—or physically, to sum up, that the passage of a direct 
current I has had no effect whatsoever on the fluctuations to be observed. 

In particular, then, all equations subsequent to (5) in Meltzer (B) are 
spurious and consequently the suggested comparison with the experiments. 
of Bernamont (1937) is invalid. 

Let us remark that it may prove possible to set up a phenomenological 
theory of current noise in thin metallic films and semiconductors 
proceeding from (our) equations (5) by making appropriate assumptions. 
about the last term (in either form) in any particular physical model. 
In any case, however, great care will be necessary in formulating and. 
interpreting an analysis of this type since the current observed at the 
terminals of a resistance, say, is composed of the contributions of all 
electrons in the body and a fluctuating current may well be induced by 
motion of remote charge without electrons actually crossing locally a 
specific cross-section of the body. It is im this sense that a rigid 
interpretation of N may be entirely misleading and, in particular, 
Meltzer’s “ equivalent shot-effect current ”* defined by I,=kT/Re in (A) 
has no very direct physical significance. 

A general feature which perhaps is not sufficiently widely appreciated 
is the complete, immediate equivalence of a statistical equation such as (9) 
above to the statement of the current fluctuations in terms of a frequency 
spectrum. If we are provided with the auto-correlation function 


y(t) .y(t-er) of a stationary statistical variable, y(t), then it should now 
be well known that the theorem of Wiener and Khintchine provides an 
immediate transformation to the corresponding (power) frequency 
spectrum. In addition, however, as we have shown (MacDonald (1949)) 
a knowledge of the average mean-square deviation of a variable over any 
arbitrary time-interval 7 also leads directly to the frequency spectrum : 


” 
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this concept is sometimes more natural than the auto-correlation function 
which deals with the instantaneous values of the variable concerned. 
Thus equation (9) is in fact precisely equivalent to writing : . 


B= 5 kT 4f+L . 3(/) Pa a SLO 


(where 5(f) is a Dirac function “at” zero frequency, representative of 
the direct current, I, such that jel f)df=1) 

An interesting application of these ideas arises in the original extended 
equations of Brownian movement due to Langevin (1908) and others 
(e.g. Furth (1920)). 

If x be the displacement of a free Brownian particle, of mass m, in a 
medium of viscous resistance A at temperature 0, then in time 7 : 


Foe Seem ey jy) ewheres T==m/A; 9; «.-. {11) 


thus the mean square deviation of the velocity averaged over the 
interval 7 is 


tae =f Fa am, Us ee) 
Ta 


7? A |r 


We can then immediately, as pointed out above, derive the frequency 


spectrum of the true “instantaneous ”’ velocity : 
=> 4k? df 
2— — ,—___~___}, Are ee eee |) 
f= Sl redTH = 
By integration over all f from 0 to «~, we have 
Se lg, ee eR 
2rAT" 2 m 


in agreement with the classical equipartition formula of Einstein, while 
(13), on the other hand, shows immediately that, unless our measuring 
instruments respond to, or resolve, events of natural frequency 
considerably greater than 1/T, we shall not in fact observe the “ true ” 
Brownian movement velocity as given by (14) (cf. Ramsay (1892), Exner 
(1900), Svedberg (1906)). If, however, the ee resolution time of 
our instrument is, say, T, then (13) enables us to determine the value of 
$v, Which should in fact be observed. Thus: 
4hO pts df 
om (1+ (2nfT)?) 

2 2nT a 
= = (arctan): ee eee Ld) 


TM 1 


obs 


80 Sbs =) 


awe 4h /'T 
or iol< !,, SS = (mr): Me Ms eer es sin LO’ 0) 
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In examining the early literature to find measurements for comparison. 
with equations (15), I was unable to find entirely satisfactory data, 
generally because no very precise record of T, was kept. However, those 
of F. M. Exner (loc. cit.) enable a good comparison to be made. Since 
T=m/d and A=6rya, for a small spherical particle, of radius a, observed 
in a medium of viscosity 7, we have : 


8 Te ie a sss Tire (16) 


It appears that Exner followed visually the track of the particle 
continuously, in so far as he was able. It therefore appears reasonable 
to take T;~0-1 second, with an upper limit certainly of 1 second. 
Now, 

= 3-5< 10 cm. 

7 = 10? poises. 

G== 293" K., (20°C. }. 
Whence y/(5v2,,,)2v5 X 10-4 cm./sec. (with a lower limit of 1-6 1074 
em,/sec,). The velocity observed by Exner at this temperature was in 
fact 3-2 x 10-4 em./sec., and the agreement is thus very satisfactory. 

This analysis appears also of interest because it has been stated in the 

past that the velocity of a Brownian particle as observed visually has 
only a subjective significance. 


I should like to thank Professor M. H. L. Pryce and Mr. D. F. Gibbs 
for helpful discussion. 
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_ ABSTRACT. 


_ The frequency of production, and the angular and energy distributions 
of single cosmic ray protons produced in photographic emulsions exposed 
at 3650 m. are given. It is shown that, after corrections for protons 
emitted from stars have been made, the results agree with those given 
in a previous paper (Lattimore 1949). 

The rates of production of fast protons in paraffin wax, carbon, 
aluminium, and lead are deduced from measurements of the number 
of protons observed to cross 100u C.2 emulsions surrounded completely 
by these materials. From the results for paraffin wax and carbon, the 
rate of production in hydrogen is also deduced. 


INTRODUCTION. 


In a previous paper (Lattimore 1949) it has been shown that the number 
of proton tracks observed crossing 100u C.2 emulsion is several times 
greater than the number to be expected if these tracks were due solely 
to protons from cosmic ray stars in the surrounding material. From the 
number of such tracks, and the number of protons observed to end in the 
emulsion, it was deduced that 62 protons were produced/c.c/day in 
emulsion at 3650 m., and that their energy of production was 30 MeV. 
or more. In the calculation of the number of protons crossing the 
emulsion, star protons of energies greater than 30 MeV. were ignored. 
The number of such star protons was therefore included in the value of 
63/c./c./day. 

It must be emphasized that at the time of these previous experiments, 
emulsions sensitive to tracks of minimum or near minimum ionization 
were not available. The C.2 emulsion was not of sufficient sensitivity 
for protons with energies of 30 MeV. or greater to give recognizable tracks. 
Now, however, emulsions sensitive to minimum ionization tracks are 
available, and it is possible to study the actual beginnings of the single 
proton tracks. The present paper is concerned with the origin of these 
single protons and their angular and energy distributions. The relative 
numbers of protons produced in various materials is also given. 


* Communicated by Sir George Thomson, F.R.S. 
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PART A: 
1. Frequency of Production. 

Ilford G.5. 2004 emulsion on glass were used. The emulsion was 
searched twice, using a X95 oil immersion objective. A total of 89 fast 
protons, other than those produced by stars, were observed to start 
within the emulsion. Of these, 74 occurred singly. Five cases were 
observed where two fast protons were produced together. Of these, two 
occur near the surface, and may be coincidences. One case was observed 
of what may be a proton of energy greater than 100 MeV. either producing 
a proton of 60 MeV. and vanishing itself, or losing considerable energy 
in passing through a nucleus, without producing visible tracks. Two 
cases were found where a fast proton apparently produces a second proton, 
and escapes itself. The value of 25+-2/c.c./day is obtained for the rate of 
production in emulsion of these fast protons at 3650 m. In addition to 
these protons, we now observe many fast protons emitted from stars. 
These were not observed in the C.2 emulsions, in general, because of the 
low sensitivity. In Table I. is shown the number of protons of energy 
greater than 30 MeV. associated with stars. 


TABLE I. 
Mean No. of Protons of No. of Protons of 
Star Size Energy greater than Energy greater than 
30 MeV. per star 30 MeV. produced/c.c./day 
2, 3, or 4 prongs al 20-5 
4 prongs 1-26 7:5 


This gives a total of 28 protons of energy greater than 30 MeV. 
produced/c.c./day in emulsion at 3650 m. from stars. 

Considering, therefore, both the protons produced singly, or in processes 
other than normal star production, and those protons of energy greater 
than 30 MeV. produced in stars, we find a total rate of production of 
53/c.c./day in‘emulsion at 3650 m. These are protons actually observed 
to start in the emulsion. 


2. Angular and. Energy Distribution. 


The angles between the projection of the proton tracks on the surface 
of the emulsion and the vertical were measured by means of a protractor 
attached to one of the eyepieces of the microscope. The angles of dip 
were also measured, by means of a graduated fine focusing adjustment. 
This was calibrated by measuring the thickness of the emulsion, which 
was assumed to be 200 microns. The lengths of the projections of the 
tracks on the surface of the emulsion were measured, and the number of 
grains along the tracks. From these four measurements we could 
calculate the true angle between the proton and the vertical, and also 
true grain density. Not all of the tracks found were suitable for measure- 
ment, either because of shear in the emulsion, or because they dipped too 
steeply for us to perform a reliable grain count. 

A grain density-energy curve for protons was constructed by measuring 
grain densities along p mesons which stopped in the emulsion, and using 


No. per 20 MeV. interval. 
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__ Energy in MeV. 


the relation dN/dR=e?f(8), where dN/dR is the grain density, ““e”’ the 


charge of the particle, and f(8) some function of the velocity of the 
particle. 


Bags lw: 


No. per 20 MeV. interval. 
a 
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Fig 16. 
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Energy in MeV. 

The energy distribution of those protons produced other than in stars 
is shown in fig. la. It has a maximum between 40 and 60 MeV., and 
falls to zero at 120 MeV. The energy distribution for both protons from 
stars and protons produced singly is shown in ng b. 
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The distribution shown in fig. 1a falls to zero at an energy of 20 MeV. 
Single protons are produced with energies below 20 MeV. (see previous. 
paper), but most of these have very low energies (5 MeV.) and it is probable 
that these are true single prong stars. The cut-off at 20 MeV. is therefore 
somewhat arbitrary. 

Fig 2a. 


28 
24 


20 


No. per 20° in'erval. 


Oe 


In fig. 2a is shown the angular distribution of the single protons, of 
energy greater than 20 MeV., at production, and in fig. 26, the angular 
distribution of all fast protons produced, i.e. both those produced in 
stars and those in other processes. In these figures @ is the angle between 
the track and the vertical. 


3. Comparison with Previous Results. 

In a previous paper (Lattimore 1949) it was deduced that 62 fast 
protons were produced/c.c./day in emulsion at 3650m. This value was 
deduced from the number of proton tracks observed to cross the emulsion, 
and from an estimate of the sensitivity of the emulsion used. It is clearly 
difficult to assess this latter quantity with any accuracy, and it is probable 
that the sensitivity was underestimated in these earlier experiments. 

In fig. 3 is shown the angular distribution of all fast protons produced. 
N(#) represents the number of protons emitted/unit solid angle at an 
angle of 6° or (180—6)° to the vertical. The full line is the distribution 
obtained in the previous paper (fig 46). The agreement is reasonably 
good, although the statistics are poor. 
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No. per 10° interval. 
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Panes GB. 
Variation of Frequency of Production with Material. 

For this part of the experiment, it was necessary to revert to the 
original method ‘of counting, single protons which crossed the emulsion, 
and which were therefore produced in the surrounding material. 
100 C.2 emulsions on film were used and these were surrounded on all 
sides by the various materials and exposed horizontally, emulsion 
uppermost, at 3650 m.C.2 emulsions were used in preference to G.5 
because the latter records all ionizing particles, including fast mesons 
and electrons. It was found that the number of protons crossing the 
emulsion was independent of the surrounding material, in agreement 
with results of Perkins (1947) and Morand and Cuér (1949). 

In Table II. are shown the number of protons observed to cross the 
emulsion with the different surrounding materials. Also shown are the 
calculated rate of production in the various materials. For the calculation 
of the rate of production, it is necessary to know the sensitivity of the 
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emulsion, by which we mean the energy of the fastest proton that can 
be observed to cross the emulsion (see previous paper). This was taken 
as 35 MeV., a value which brings the rate of production given in our 
previous paper (Lattimore 1949) into agreement with the rates found 
in emulsion in the present experiment. We can then find the rates of 


Fig. 3. 


30 


20 


°5 10 20 30 40 SO 60 70 80 90 


Ge, 


production in the various materials by using the formula given in the 
previous paper. The ranges were calculated from the equation 


5 
dE /da=4-15p(Z/A) u,/B? | toe pe —F| 


(Bethe 1937), where p, Z, and A are respectively the density, atomic 
number and atomic weight of the stopping material, ,., the mass of the 
electron, 6 the velocity of a proton of energy E, and a the range of the 
proton in gms./em.’. [is the ionization potential of the stopping material 
(see Table II.). This equation was integrated graphically. Since the 
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calculated absolute values for the rate of production depend upon the 
value taken for the sensitivity of the emulsion used, they are perhaps 
open to doubt, but the relative values should be approximately correct. 


TABLE II. 
No. crossing N 
Material emulsion I(eV) a ae ie 
Arbitrary units STRICTED: 
Pb. 38-712 1150 6-5 40-5 
Al. 37-442 Wes 9-65 40-7 
C. 40-4--2 81 15-1+1-1 
CH. 42-7412 18-6 +1-4 
H. 13-5 39-6+12 


Discussion of Errors. The principal error involved, other than purely 
statistical errors, is introduced by taking the emulsion thickness as. 
200 wu. Thethickness was unfortunately not measured before development. 
Measurements on two further batches give errors of 6 per cent and 8 per cent 
in the thickness. Any error will affect the angular distributions only 
slightly by alterig the angle of dip, but affects the number of single 
tracks produced/c.c./day directly. 


CONCLUSIONS. 


This work has confirmed our earlier conclusions given in the previous 
paper, that protons are produced singly, and that they are mainly 
directed downwards. 

It would seem reasonable to suppose that these protons are produced 
by the same primaries as produce stars, since the number both of stars. 
and single tracks vary in the same way under ice (Harding e¢ al. 1949) 
and under air (Lattimore unpublished). It follows that any theory of 
single proton production must take account of star production also. 
Attempts are being made to produce such a theory, and any results wilh 
be published later. 
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ABSTRACT. 


A technique is described for controlling a Wilson cloud chamber by 
-an internal counter whose cathode consists of a hexagonal array of rods 
and whose sensitive volume is not separated from the main cloud chamber 
‘volume by any wall and in fact forms part of the usable chamber volume. 
The counter has no stopping power beyond that of the gas unless the 
incident particle chances to hit one of the rods. Events in which heavily 
ionizing particles are produced in the cloud chamber and events in which 
low energy particles traverse the chamber can thus be photographed. 
Discrimination between single fast particles or events of high ionization 
-can be made at will. Altogether 26 slow particles (excluding electrons) 
have been photographed—corresponding to a vertical intensity of 
9:3 x 10-5 slow particles cm." min.~! sterad.—1 at sea level. Details of a 
proposed arrangement for photographing nuclear stars in the gas of the 
chamber are given. The technique is applicable to the study of the 
discharge phenomena. of counters and some preliminary photographs of 
these phenomena are presented. The counter is quite stable and the 
technique of operation is not unduly difficult. 


$1. INTRODUCTION. 


Ix many cloud chamber experiments it would be an advantage if the 
chamber could be triggered on the occurrence of a heavily ionizing event. 
A proportional counter operated within the chamber offers a means of 
detecting such events, but for the observations of processes involving or 


caused by very low energy particles the absorption in the wall of the | 


counter is an important factor. The possibility of employing a counter 
in which this absorption was reduced to a minimum by means of a cathode 
assembly consisting of rods was investigated by Hodson and Page in 1948 
with a view to detecting and photographing nuclear stars occurring in 
the gas of a cloud chamber||.. The idea of using a cathode of this form 
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| After this initial work, attention was drawn to a proposal of Cohen to use a 
cloud chamber as an ionization chamber for this purpose (Cohen 1948). The 
success of this method has since been reported (Cohen 1949a, b) 
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was obtained from the work of McCusker (1942) who had used such an 
arrangement for Geiger-Miiller counters in glass envelopes. After 
experimenting with various rod assemblies McCusker showed that a 
hexagonal arrangement of rods was the most satisfactory. Hodson 
and Loria (1949) succeeded in controlling a Wilson cloud chamber by 
an internal counter having a cathode of this type but whose sensitive 
volume was not separated from the main chamber volume by any 
envelope and in fact formed part of the usable chamber volume. 

This paper gives details and progress of this technique, together with 
some preliminary results obtained. 


§2. EXPERIMENTAL DEralIzs. 
(a) Constructional details of the counter : AM 
The counter at present employed is 25 cm. in length and 4 cm. in 
diameter and is illustrated in Pl. XXV. .The cathode rods are of steel, 
2 mm. in diameter, supported at each end by a brass framework, the whole 


Rig, 


Diagram of wire mounting. 
Anode wire. 

Shield: 

Rubber gasket. 
Glass cylinder. 
Threaded support. 


EOD 


cathode assembly being chromium plated. It is mounted across the 
middle of a cloud chamber of 28cm. diameter and 9 cm. depth. The 
anode wire is attached to two metal supports (for details see fig. 1) passing 
through holes in the glass cylinder of the chamber, the holes then being 
sealed with small rubber gaskets. The compressibility of the gaskets 
also provides a useful means for the final tensioning of the wire. Both 
tungsten and platinum wires (diameter 0-1 mm. and 0-17 mm.) have been 
used successfully, but, for reasons to be given later, platinum wires are 
preferable. In future applications it is proposed to use cathode rods of 
smaller diameter with surfaces blackened to reduce the amount of light 


scattered. 
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(b) Chamber filling : 

In normal cloud chamber operation advantage is taken of the readiness: 
with which electrons produced by ionizing particles attach themselves to 
oxygen molecules in the gas to form heavy negative ions. In counters, 
however, the production of gas multiplication depends on the electrons 
remaining freely mobile, and electron attachment must be avoided. For 
this reason precautions must be taken to exclude oxygen and water 
vapour. The cloud chamber is therefore filled with commercial argon 
(stated to be of 99-8 per cent purity and oxygen-free) and saturated with 
absolute ethyl alcohol vapour—a mixture which gives good cloud chamber 
tracks and at the same time is a satisfactory filling for the counter. As 
it is not possible to evacuate the chamber, the air initially present is. 
removed by partial extraction with argon : argon is admitted toa pressure 
of 1-5 atmospheres and the pressure is then released to one atmosphere. 
This procedure is repeated at least ten times and the chamber then 
finally filled with argon to the required pressure of about 118 cm.Hg in 
the compressed position. 

If there is no appreciable negative ion formation the electrons produced 
along the path of an ionizing particle are rapidly swept away by the 
general clearing field applied across the chamber. The tracks observed 
when the cloud chamber is expanded thus arise from positive-ion 
condensation only. 

It is found that if no further precautions are taken the pulse size of 
the counter decreases steadily over a period of several days after the 
cloud chamber has been filled—presumably due to the slow release of 
oxygen and water vapour from the walls and velvet background of the 
chamber and to the slow diffusion of air into the chamber through the 
rubber diaphragm. The increasing amount of electron attachment which 
thus occurs results in an undesirable doubling of the tracks, showing that 
condensation is now occurring on both positive and negative ions. (A 
similar difficulty has been experienced by Cohen (1949 b) with his. 
combination cloud- and ion-chamber). ‘The status quo can be restored. 
temporarily by reflushing the chamber with argon and for the present 
purposes it has proved sufficient to do this once a week. However, a 
method of continuously removing oxygen and water vapour is essential 
for a permanent and reproducible pulse size in the counter, and we are 
now investigating various possible methods. 


(c) Operation of the counter : 

The counter will operate satisfactorily either in the proportional region 
or in the Geiger region. With the pressures used (118 cm. of mercury) 
and a wire diameter of 0:17 mm., suitable proportional operation is 
achieved over the range 1950-2300 volts, while the threshold potential for 
the Geiger region is about 2500 volts. 

The cathode rods of the counter are earthed and the anode wire is 
connected to a positive high voltage supply through a 2 megohm resistance. 
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It will be seen later that the most suitable method for detecting particles 
inside a cloud chamber is to operate in the proportional region, using a 
fairly low gas multiplication and high external (valve) amplification. In 
this mode of operation the negative pulses obtained from the wire are 
amplified by a high gain linear amplifier (gain variable between 500 
and 3-10°) incorporating a biased diode discriminator at the output stage, 
which, in conjunction with the gain setting, determines the minimum size 
of pulse that will be accepted. 

Particular care is necessary to avoid spurious pulses at the input stage 
of the amplifier. Any appreciable leakage from high voltage points 
associated with the anode wire will give rise to spurious pulses, but by 
mounting the anode wire directly on the glass cylinder of the chamber 
so that the leakage path is made as long as possible and by using poly- 
styrene supports for all high voltage components, leakage is readily 
reduced to a negligible level. A silicone varnish (trade name “ Aquatect ”’) 
has proved very effective in reducing residual leakage across the surface: 
of the glass cylinder. These initial precautions being taken, the counter 
has been operated continuously for many weeks without any trouble due 
to spurious pulses. 

For these initial investigations the counter has been operated in 
coincidence with two Geiger counters placed above the chamber in a 
vertical telescope arrangement. The output pulse from a threefold Rossi 
coincidence circuit is used to expand the cloud chamber, and within a 
few microseconds (q. v.) of a coincidence being detected the high voltage 
on the wire is removed electronically. The voltage is thus removed before 
the positive ions along the path of the incident particle have had time 
to move any appreciable distance and the track is therefore visible and 
undistorted even within the sensitive volume of the counter (Photograph (a), 
Pl. XXVI.). The small “bead” of ionization observed on the wire near 
to the track of the incident particle consists of the positive ions formed by 
the gas multiplication which occurs near the wire during the discharge 
of the counter. By using high values of the amplifier gain and quenching 
rapidly it is, in fact, possible to remove the voltage on the wire before the 
process of gas multiplication has reached completion—the “bead” 
observed is then correspondingly smaller and may even be too small to 
be visible on the cloud chamber photograph. However, it is often 
desirable to know which particle discharged the counter and this is 
possible if a track can be associated with a “ bead” on the wire. For 
this reason the removal of the anode voltage is usually delayed by approxi- 
mately 20 microseconds; this allows the discharge to develop until it 
is quenched by a space-charge reduction of the field near the wire (Stever 
1942), and a “ bead ” of reasonable size is observed. 

The columns of ionization marked W in Photographs (a) and (6), 
Pl. XXVI., and also visible in the subsequent photographs arise as a result 
of surface irregularities of the wire. Ifthe wire is not smooth or if particles 
of dust adhere to it an ‘electric wind”’ or corona phenomenon produces 
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continuous streams of positive ions moving away from the wire and 
occurring at the same points in successive photographs. This phenomenon 
does not interfere directly with the normal operation of the counter but 
the streamers are rather undesirable from the point of view of clear 
cloud chamber pictures. Fortunately it is possible to avoid most of them 
by cleaning the wire carefully. The initial cleaning consists merely of 
washing the wire with alcohol with a fine brush. Platinum wires are 
preferable to tungsten as they have a smoother surface and may be safely 
glowed in air after closing the chamber as the final stage of the cleaning 
process. The “electric winds” have been investigated in some detail 


Fig. 2. 
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and the exact mechanism of their production will be discussed in a 
subsequent paper. They are more pronounced at higher voltages on 
the wire. 

The general circuit arrangement used for removing the anode voltage 
of the counter after a comcidence has been detected is shown 
diagrammatically in fig. 2. The valve V is normally biased beyond anode 
current cut-off so that the full positive potential is applied to the wire of 
the counter. When a coincidence is detected a pulse, delayed by about 
20 microseconds, triggers a single-shot multivibrator. This produces a 
square positive pulse (~1/3 second long) which is applied to the grid of 
the valve V. Valve V then becomes fully conducting and the anode 
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voltage of the counter is rapidly removed. A few milliseconds after the 
cloud chamber is triggered the control mechanism of the chamber closes 
a relay contact K and maintains the grid of valve V at earth potential 
until the cleaning cycles of the chamber have been completed. The 
anode voltage of the counter is therefore approximately zero during this 
period (fig. 3) and the production of fresh ionization by discharges of the 
counter and by the “electric wind ” phenomenon during the cleaning 
expansions is thereby avoided. . 

Operation of the counter in the manner described above, i.e. in the 
proportional region with removal of the anode voltage a few microseconds 
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after triggering the cloud chamber, will hereafter be referred to simply 
as ‘the normal mode of operation ’’. Other modes of operation will be 


described in §3. 


§3. Ta DiscHARGE PHENOMENA OF PROPORTIONAL AND GEIGER-MULLER 
CoUNTERS. 


The technique of operating an open counter in the gas of a ee 
chamber is one that is obviously applicable to the study of disc arge 
phenomena in counters. When a counter Is operated in the en eee 
the discharge is known to spread along the wire from the ee wl : 
particle passed through the counter (Stever 1942). Severa eta ions 
of the velocity of this discharge spread have been oe e ie Cee 
mechanisms for this process have been suggested. ee a ; son 
(1948) have published a useful review of ots work and it ne oa 
discussed more recently by Balakrishnan and Craggs (1950). special - 
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interest also is the transition between the proportional region of @ 
counter and the Geiger region. The present technique is being applied. 
to the study of these problems and an account of this work will be 
published separately. However, a few photographs will be given here to 
indicate the possibilities of the method. 

If the counter is operated in the proportional region and the cloud 
chamber is triggered without removing the voltage from the anode wire 
the positive ions formed near the wire in the discharge of the counter 
move outwards towards the cathode until condensation begins. Instead. 
of forming a ‘“‘ bead ” near the wire as in Photograph A, the positive ions. 
now appear as a “ring” of ions in Photograph (6) (Pl. XXVI.).  Photo- 
graphs (a) and (b) demonstrate directly the extreme localization of the 
discharge in the proportional region. 

Photograph (c) (Pl. X XVI.) was taken a ta rather higher anode voltage 
than Photographs (a) and (6) and is typical of photographs obtained in 
the transitional region between the proportional and the Geiger regions. 
It was obtained by removing the anode voltage about 20 microseconds 
after triggering the chamber and may be compared with Photograph (a). 
At the higher voltage it is seen that the discharge spreads along the wire 
in a series of quite characteristic and fairly regularly spaced ~ jumps ”’, 
and the positive ion space charge consists of a number of ~ beads ” 
along the wire. It has been observed that as the anode voltage is further 
increased the discharge becomes progressively more continuous along the 
length of the wire. 

From the previous description of the normal mode of operation it is. 
clear that the anode voltage is removed only when expansion of the 
chamber is about to take place. But bursts of ionization are being 
produced continually in the counter and although the initial radial 
movement of the “bead” is very rapid the time for the complete 
collection of the positive ions on the cathode rods is several milliseconds, 
owing to the high gas pressure used in the chamber. There is thus a 
probability that besides the selected particle and its associated “‘ bead ” 
of ionization there will be present tracks and “rings” of ionization 
due to events which occurred within a certain interval before the principal 
track. Examples ofthismay beseen on Photographs (e)and(f), Pl. XX VII. 
At low voltages, where the discharge is very localized, the “ rings” so 
formed are narrow and not very intense. At voltages sufficiently high to 
produce spreading of the discharge along the wire the effect of a previous 
count becomes more noticeable, due to the formation of extended “ ring ” 
systems (¢. g. Photograph (d), Pl. XX VI.) or even continuous “ cylinders ”” 
of positive ions (e. g. Photograph (c), Pl. XXVI.). These are undesirable in 
photographs intended to show interesting particle events in a cloud 
chamber, and from this point of view the most suitable operating region 
of the counter for particle detection is the proportional region. By 
keeping the gas multiplication fairly low any single “ rings ”’ present do 
not seriously obscure any part of the cloud chamber. 
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§ 4. Discussion or Rzsutts. 


We now describe a series of photographs which demonstrate the 
potentialities of the internal counter technique for detecting particles 
inside a cloud chamber. The internal counter was used in the normal 
mode of operation, with the threefold telescope arrangement described 
in §2 and with the discriminator biased towards events of high total 
ionization. 

The photographs may be divided into three types : 

Type (a). Photographs showing the internal counter being triggered 
by single heavily ionizing particles. 

Type (6). Photographs of local showers and air showers in which two 
or more fast particles of minimum ionization traverse the internal counter 
simultaneously. 

Type :c). Photographs of single fast particles which produce slow 
knock-on electrons in the gas or in the cathode rods of the counter. 

Type (a). Slow Particles. 

Photographs (e) and (f), Pl. XXVII., show slow particles, probably 
protons, triggering the internal counter. The slow particle in (f) appears 
to originate in a nuclear star in the glass cylinder of the cloud chamber. 

Photograph (g), Pl. X X VIT., shows an event which is almost certainly the 
decay of a x-meson in the glass cylinder at the bottom of the chamber. 
The decay electron re-enters the chamber, travels upwards at a rather 
steep angle, unfortunately misses the counter and passes out of the 
illuminated region of the chamber. Both tracks project back to a point 
in the cylinder. 

Photograph. (h), Pl. XX VII., shows a slow particle which stops within the 
sensitive volume of the counter. It is out of the illuminated region at , 
the top of the chamber. From the observed scattering in the gas the 
particle is very probably a slow meson. No decay electron is visible in 
the chamber and unfortunately the end of the track is obscured by the 
large bead of ionization near the wire. For a given quenching delay, the 
size of the bead round the wire increases with the amount of initial 
ionization in the counter. As the operating conditions of the counter 
which affect the size of the bead (in particular the anode voltage and the 
gas purity) are not necessarily the same in all the photographs, too strict 
a comparison of bead sizes is not justified. However, this bead is 
appreciably larger than any obtained under the same conditions. Mutual 
repulsion of the positive ions in the dense space charge before condensation 
occurs is more likely to be the cause of the unusually large diameter of 
the bead, rather than their motion under the influence of the applied 
field in the few microseconds before it is removed. The spherical nature 
of the bead supports this view. 

In Photograph (i), Pl. XXVIIL., a slow electron is seen to discharge the 
counter. This shows that it is possible to trigger a cloud chamber by 
events which involve electrons or f-particles having energies too low for 
detection by conventional techniques. 
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Altogether 26 slow particles (excluding electrons) of ionization greater 
than about three times minimum have been photographed in a total 
sensitive time of 727 hours—corresponding to a vertical intensity of 
about 9:3 10-5 em.-2 min.-! sterad.-1, that is, about 0-01 per cent of 
the total vertical cosmic ray flux at sea level. In three cases there are 
definite indications that the heavily ionizing particles are mesons (e. g. 
Photographs (g) and (h), Pl. XXVII.). Inthe remaining cases, however, 
the particles cannot be identified with certainty. On account of the 
material above the internal counter (chamber wall, external counters, etc.), 
the slow mesons observed would have incident momenta in the range 
(9-2-10-1) x 107 eV./c., while the protons would lie within the range 
(4:3-6-0) x 108 eV./c. 

It is difficult to compare the above rate with the results of other 
workers obtained under different experimental conditions, but 
Montgomery (1949) in discussing the rather indefinite evidence available 
has suggested a value of the order 610-4 cm. min.! for the omni- 
directional intensity of slow particles at sea level. Anderson and 
Neddermeyer (1936) and Widhalm (1940) have shown that the heavily 
ionizing particles are much less concentrated about the vertical direction 
than is the total ionizing component, so the unidirectional intensity is 
probably of the order 10~4 particles cm.~? min.~! sterad.~1, a value which 
is in satisfactory agreement with our observed rate. 

Rossi, Sands and Sard (1947) have measured the slow meson intensity 
at sea level in the momentum range (6—12)x 10% eV./c. and obtained a 
value 5-4x10-> slow mesons cm.-? min.! gm.-! (range in air). This. 
corresponds to an omnidirectional intensity of 14:0 x 10—> cm. min.-1 
in our momentum range. The unidirectional intensity in this momentum 
range is therefore probably of the order 5 x 10-> slow mesons cm.~? min.~! 
sterad.~!, since the angular distribution of slow mesons is likely to be such 
that there is a general concentration about the vertical. (On the other 
hand, the angular distribution of slow protons may. be expected to be 
much more isotropic, since many are secondary particles emitted from 
nuclear evaporations in local material (Anderson and Neddermeyer, 1936). 
Rochester and Bound (1940) measured the number of slow protons in the 
range (3-9-5-5)x10%eV./c., and from their results Rossi (1948) has 
deduced a rate of (3+1)x10~ slow protons cm.-? min.—! sterad.—!. 
These results for the intensities of slow mesons and slow protons taken 
together give a total slow particle rate which is again in fair agreement 
with the present results for comparable momentum ranges. It seems 
reasonable to suppose that the numbers of slow mesons and slow protons 
observed with the present experimental arrangement are of a similar 
order. 


Type (b) Showers. 


The internal counter obviously cannot distinguish between one slow 
particle having an ionization three times minimum (say) and three 
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simultaneous particles of minimum ionization. The threefold telescope 
arrangement is therefore sensitive to showers, either incident from the 
air or locally produced. Examples of events of this type are shown in 
Photographs(j)and(k)ofP].X XVIII. Itwillbenoticed that when particles. 
pass simultaneously through the sensitive volume of the counter, each 
produces a local discharge near the wire. A particle traversing the 
counter some short time before the cloud chamber was triggered would be 
associated with a “ring” of positive ions, more or less expanded (e. g. 
Photograph (f), Pl. XXVII.), while a particle arriving after the counter 
voltage was removed would produce no discharge. This ability to 
determine whether two or more particles are contemporary decreases 
the effective resolving time of the chamber and could be employed with 
advantage in certain experiments. A preliminary estimate of the time 
interval required for the positive ions to move outwards from a “ bead ”’ 
to a recognizable “ring” gives the effective resolving time of the order 
100 microseconds, whereas the resolving time of a normal chamber of 
this size is of the order 1/50 second. 


Type (c). Knock-on electrons. 


Photograph (/) of Pl. XXVIII. is typical of many photographs obtained 
with the experimental arrangement described : a fast cosmic ray particle 
produces a slow knock-on electron in the gas. The ionization of the fast 
particle alone would not have given a sufficiently large pulse to pass the 
discriminator, and it is clearly the slow electron which contributed most 
to the large pulse obtained. Slow knock-on electrons in the gas or 
cathode rods of the counter occur fairly frequently, and must be 
responsible for many of the large pulses obtained from the normal 
enclosed type of proportional counter. (This conclusion was previously 
reached by Cook (1949) in some recent work using proportional counters). 
In experiments using conventional proportional counters there is often 
an undesirable uncertainty in knowing exactly what events triggered the 
counter. This is avoided in the present technique. 


§ 5. Tue InrerNAL CounTER TECHNIQUE FoR DetTectinc NUCLEAR 
DISINTEGRATIONS. 


The internal counter need not, of course, be operated in coincidence: 
with an arrangement of external Geiger counters as in the present 
work. For example the cloud chamber may be triggered by single pulses. 
from the internal counter alone. With the discriminator set to record 
only events of high ionization a short run on these lines produced 
photographs mainly of contamination a-particles. This arrangement would 
undoubtedly be also sensitive to nuclear disintegrations occurring in 
the gas of the cloud chamber. A rough calculation based on photographic 
plate evidence (George and Jason 1949) gives the rate of nuclear stars in 
the gas of the cloud chamber as four per day. A comparison of this with the 
measured «-particle rate of the counter, gives an expected yield of only 
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one star photograph per 500 photographs of contamination o-particles. 
(Cohen has reported a similar difficulty with his combination cloud- and 
ion-chamber). To bias the counter towards events of higher ionization 
than that of a contamination «-particle would of course mean that many 
stars would not be recorded, especially those whose origins lay outside 
the sensitive volume of the counter. 

However, an arrangement consisting of two or more open counters 
operated in coincidence inside a cloud chamber would be sensitive to 
stars but not to contamination «-particles provided the counters were 
separated by a distance greater than the range of such «-particles 
in the gas (~5cms.). At the same time there would be no necessity 
to bias towards events of very high total ionization. It is hoped to employ 
this arrangement in the near future. 


§ 6. CONCLUSIONS. 


The technique of controlling a Wilson cloud chamber by an internal 
counter of the type described has the following features : 

(a) The counter may be used in the proportional region and set to 
record either single particles of minimum ionization or heavily ionizing 
events. 

(6) Unless a particle happens to hit one of the cathode rods the counter 
has no stopping power beyond that of the gas filling. Slow particles 
traversing the cloud chamber and events producing slow particles within 
it can therefore be photographed. 

(c) The counter may be operated in such a way as to obtain tracks 
throughout the sensitive volume of the counter, and thus the useful 
volume of the cloud chamber is not reduced. This is an important 
feature, especially since it is proposed to put two or more of these counters 
inside a cloud chamber to detect nuclear disintegrations in the gas. 

(d) It is possible to observe the exact event which triggered the 
counter. / 

(e) The technique is one which is readily applicable to the study of the 
discharge phenomena in counters. 

(f) The internal counter affords a means of testing the contemporary 
nature of particles which pass through it, and the effective resolving time 
of the cloud chamber is thereby greatly reduced. 

(g) The technique is not difficult and, given the few elementary 
precautions described, the counter is quite stable. 
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SUMMARY. 


In order that momentum and energy be conserved in the process of 
n—p-decay at rest, it is necessary to assume the emission of a single 
additional non-ionizing particle in a direction opposite to the initial line 
of motion of the w-meson. There exists strong experimental evidence 
for the conclusion that the neutral particle has small or negligible rest- 
mass. Hence. we may suppose that it is either (a), some light neutral 
particle possibly identical with the neutrino of B-decay, or, (b), a photon. . 

Some 250 suitable cases of 7—y-decay in Ilford G5 Nuclear Research 
plates have been examined with a view to testing the hypothesis (6). A 
careful search has been made in the requisite direction, for associated 
pairs of electrons. In a total scan-length of 38 cms. no such case was 
observed. Statistical analysis shows that, if indeed a photon were emitted,, 
the probability that it would fail to produce a pair under the conditions. 
of the experiment was 4x 1073. 

Since this result would appear to exclude the possibility of photon 
emission, it is inferred that the 7-meson normally decays in accordance 
with (a). 

It is shown that the present experiment has much in common with 
those devised to demonstrate the existence of the neutrino of B-decay. 


§1. THe Decay oF z-MEsons. 


In the p-decay of z-mesons as observed in photographic emulsions,, 
it is known that the distribution in range of the p-particles is consistent 
with the assumption that they are always emitted with constant kinetic: 
energy, the departures of the individual values of the range from the: 
mean being due to “straggling ’”’. It follows from this result that the 
disintegration is a two-body process, momentum, energy and charge being: 
conserved by the emission of a single neutral particle, which leaves no, 
recognizable trace in its passage through the emulsion. 


* On leave of absence from University College, Cork. 
t Communicated by Professor C. F. Powell, F.R.S. 
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A study of many examples of the process as it ordinarily occurs in the 
plates has shown that the spatial distribution of the trajectories of the 
#-meson is isotropic with respect to the final direction of motion of the 
disintegrating 7-meson (Lattes et al:, 1947, Bradner 1949). Taken in 
conjunction with the constancy in the range of the p-meson, this result: 
proves that the 7-mesons are commonly moving with very small velocities 
at the instant of decay. If it were otherwise, the range of the ,.-mesons 
would vary widely according to their direction of emission with respect 
to the line of motion of the parent 7-particles. The neutral particle formed 
in the decay will therefore be projected in a direction opposite to the 
initial line of motion of the u-meson. 


Mass oF THE NEUTRAL PARTICLE. 


By application of the conservation laws it may be shown that the 
rest-mass (%)) of the neutral particle bears a ratio (P) to that of the 
f-meson given by 

Bi Ore 206 al) a eet) 
where Q=m,/m,, and By is the value of (1—f?)-¥? proper to the 4-meson 
at the instant of its formation. The following Table calculated from (1) 


gives, for certain pairs of values of Q and B,,, some real values of the mass 
of the neutral particle expressed in terms of the mass of the electron, m,. 


TaBLeE I. n,/m,. 


Bu ‘ 
Q 1-036 | 1-037 1-038 | 1-039 | 1-040 


1-325 21-2 18:3 14-7 10-0 — 


1-321 18-7 15-3 10-8 1-2 — 


feu ee |) 11.6 42 = i 
1-315 | 142 Gioia sé a 
1-313 | 12-3 6-0 _ - = | 
1-312 | 11-3 3-4 xo az. 7 
1-311 | 10-1 i | ia s | %5 


The most recent measurements (Berkeley Group 1950) give 
m,=276+6m,; m,=210+4m, ; Q=1-315. They are in good agreement 
with those found for cosmic-ray mesons in this laboratory (Franzinetti 
1950). The observed range of the ~-meson in Ilford G5 emulsion is 
595+10 microns. Using the established range-energy relationship 
(Fowler 1950), and setting m,=210 m,, we find B,=1-038. 
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It will be clear from Table I. that an exact determination of n/m, in 
the region n, > 0 requires very precise knowledge of the values both of 
Q and of Bu. Nevertheless, the accuracy of the experimental measure- 
ments seems to be such that we may conclude with confidence that the 
neutral particle has a mass not greater than a few electron masses. Since 
there is no certain evidence for the existence of neutral particles with a 
mass of this order, it appears preferable to suppose that it has negligible 
rest mass. This assumption is consistent with all present experimental 
information and will be made in what follows. 

In order to conserve momentum the neutral particle, if of negligible 
rest-mass, must have energy—30-1 MeV. Such a particle might be 
(a) the neutrino of B-decay, (6) some light neutral particle other than the 
neutrino or (c) a photon. 

If it is the neutrino of f-decay, its interaction with matter would 
.appear to be negligible (Nahmias 1935, Ellis and Wooster 1927), and thus 
it would be expected to escape detection in the photographic emulsion. 
Evidence based on the shape of the f-ray spectrum of tritium indicates 
that its rest-mass is probably less than 1/30 m, (Nielsen 1949). 

If the neutral particle is analogous to the neutrino, but not identical 
‘with it, it is impossible to state, a priori, in what manner it will interact 
with matter, save in respect of its energy loss by ionization which will 
‘be negligible. 

On the other hand, if the particle is a photon, its interaction with matter, 
though small, will not be negligible. For a photon of 30 MeV., the most 
important interaction will lead to the creation of electron pairs. Later, 
it will be shown that the cross-section for the production of a single Compton 
electron is only 1/7 that for pair production. For this reason we shall be 
justified in disregarding such of these as may be produced under the 
conditions of the experiment. It has been the aim of the present investiga- 
tion to detect pairs of electrons, if they exist, by scanning under high power 
in a direction opposite to that of the initial path of the »-meson. The 
observations would appear to exclude the hypothesis that the neutral 
radiation emitted in the decay of the «-meson consists of photons, and 
‘strong evidence is thus provided for the existence of some form of 
“neutrino”. That this is identical with the neutrino of B-decay is not 
certain, but the two particles will be shown to have several properties 
in common. 


§ 2. EXPERIMENTAL METHOD, 
(a) Geometry of Problem. 


The search for electron pairs associated with 7-decay was carried 
out in the batch of Ilford G5 electron-sensitive plates referred to in 
previous communications from this laboratory (Camerini et al., 1950). 
About 400 cases of z-u-decay were examined. On each of these, 
measurements were made of the position in the emulsion of the origin 
of the p-meson and of its apparent angle of dip using the first 25u of track. 
Tracks of which the dip was greater than 15° were then rejected. 
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The emulsion of these plates was 400, thick before processing, and. 
after processing, 1604. It has been shown (Fowler 1950) that for regions. 
of the plates more than 1 cm. from the edges, this shrinkage takes place: 
only in a direction at right angles to the plane of the emulsion and is a 
uniform linear transformation in this coordinate. We define the. 
shrinkage factor t,,/t,—S—400/160 for the plates concerned. Here, and 
subsequently, the subscripts w and p will refer respectively to the- 
unprocessed and processed emulsions. 

Fig. 1 represents a vertical section of the emulsion. The path of the 
neutral particle through the unprocessed emulsion will be OE=t, O being 
the origin of the 4-meson, and E the point at which the particle leaves. 
the emulsion. The angle of dip, d=arc tan PQ/25, PQ being measured. 


Fig. 1. 


in microns. We have the following relationships valid for processed. 
and unprocessed emulsions : 


tan #,=S tan #3, ie Tet ine! tes oe Se t(2) 
y=d, cot d,=d, cot #, eye tatess cor 1h)! 
MA SAAC A) i eee enemas a9 


In the actual experiment the most convenient distance to measure 
is y. For each track, ten measurements of dip were made by means of 
an eyepiece scale, and the graduated focusing drum of the focusing scale 
of the instrument. The distance d, was measured for the same setting 
of the instrument. From these data, the value of y was computed using 
equation (3). The path-length thus estimated is less than the actual 
path ¢,, in the unprocessed emulsion, but it may be seen from (4) that the 
defect of y from ¢, is small for the values of } included in the measurements, 
and further, it may be allowed for. 


(b) Method of scanning and allowance for multiple EE: 

The scanning of the plate was carried out under magnifying power: 
objective x 45, with binocular eyepieces x 8. A hair-line in one eyepiece 
was set so that its direction was parallel to the y-axis of motion of the 
microscope stage. The plate was then rotated until the hair-line appeared 
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to lie evenly along the mean-line of the grains in the first 25. of the 
p-meson track. The field of view extended to 1004 on either side of the 
hair-line. 

The path followed by the neutral particle will be tangent to the 
-meson track at the point of origin. In general, because of multiple 
scattering, the directions of the mean-line and of the tangent will not 
coincide. The angle between these directions will be normally distributed 
about the direction of the mean-line and the variance of the distribution 
may be computed from the momentum of the scattered particle and the 
composition of the emulsion (Rossi and Greisen 1941). Calculation shows 
negligible probability. that the projection of the tangent lies outside a 
sector of angle 2° on either side of the mean-line, with centre on the 
origin of the p-meson. For values of y<2500p it may be seen that the 
whole of this sector lies within the 200, strip defined by the field of view. 
For longer scans, 1 cm.>y>0-25 cm., two additional 200 u strips were 
‘scanned on either side of the central strip. 

The estimate of y is, in similar manner, uncertain because of the effect 
of multiple scattering. Referring to fig. 1 it will be clear that this 
estimate, say y,, depends on a measurement of the dip of the projection 
on a vertical plane, of the chord of the first 25 of u-meson track. A 
correct estimate, say y,, would be based on the dip of the tangent at the 
origin. It may be shown that the angle between chord and tangent has 
a normal distribution about the direction of the chord, and has variance 
(i?) =4¢¢7) (Rossi and Greisen 1941). <¢?) is the mean square value 
‘of the projection of the angle between the tangents at the end-points 
of the segment. The value of (<?))!? for the conditions of the experi- 
ment was calculated to be 1-8°. 

For any individual track, y, may be greater or less than the observable y,. 
For tracks of sensible dip, because of the small half-width of the 
distribution of 7, this difference is small or negligible. It is otherwise 
for tracks of small dip, and especially for “ flat” tracks (#,-~0), and 
must be taken into account. This was done by finding the distribution 
function for y in terms of d, and & (fig. 1) assuming normal distribution 
of 4. The mean value <y), of y was obtained for this by graphical 
integration. For tracks of measurable dip the scan was continued to 
Ye>Y- (fig. 1) such that the mean value of y, computed as above, was equal 
to y,. For “ flat” tracks, which by their nature make large contribution 
to the length scanned, it was found possible to express (y) in terms 
of a closed expression involving the tabulated functions Erf(x) and Ei(a). 


§ 3. EXPERIMENTAL RESULTS. 


The total number of z-~ decay events selected as suitable was 
253. The sum of all the effective scan-lengths Xy, computed as above 
was 36cms. The corresponding true path-length in the emulsion 
+t, was 38cms. 35 pairs of particles, probably electron pairs, were 
observed to have their points of origin within the scan strip-width 
(but not necessarily within the sector). The total area of the scan was 
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1-32 cm’. This yields a density of 26 pairs per cm2. of these plates. Other 
workers in this laboratory (Carlson et al., 1950) report a density of 
24 pairs/em®?. The agreement is satisfactory. Neither figure has been 
corrected for failure to observe the small proportion of very steeply 
dipping pairs which are bound to escape detection. 

For 24 cms. of this scan length each field was scanned twice and for 
the remaining 12 cms. the field was scanned once only. The corresponding 
areas of scan were calculated to be, respectively, 1-0 cm?. and 0-32 cm2. 
The number of pairs observed in the first series was 28 and, in the second 
series, 7. It may be concluded from these figures that there was no 
significant reduction in the efficiency of detecting pairs introduced by the 
change to single scanning. The apparent disproportion between 
scan-length and scan-area arises from the facts that the ratio scan- 
area/scan-length increases with the cube of the scan-length, and that the 
incidence of long scans was slightly smaller in the second series. 


Theoretical Cross-Sections. 


It has been shown (Heitler 1944) that, for high energies, both 
electrons of a pair are emitted in the forward direction. The solid angle 
Q into which the electrons are emitted is of the order 2~m,c?/T, where 
T is the total energy of the pair. For 30 MeV. photons Q~1/60. 

The cross-section for pair-production o,,;, by 30 MeV. photons 


=7Z? x 5-71 x 10-28 em?, and for Compton Effect 


—3>x 10-2Z x 6-57 X 10-25 em?2. 


Fcomp 


From data supplied by the manufacturers of the G5 plates we may 


calculate Z?=451, Z=13:1. The number of atoms per c.c. of emulsion 
= 8-37 102. We obtain from this 0,,;,/¢c¢omp=7, and a probability of 
pair production per cm. path in the emulsion=0-152. Hence, if the 
neutral particle is a photon, we should expect to observe, on the average, 
5:77 electron pairs “‘ associated ” with the cases of 7--decay examined. 

The criteria of ‘‘ association”? are somewhat complicated by the 
multiple scattering of the »-meson, but for present purposes, may be 
formulated with sufficient accuracy as follows : 

(1) The origin of such a pair must be found within a solid angle con- 
taining OE (fig. 1) of such magnitude as to allow, in any individual case, 
for uncertainty as to the direction of the tangent at the origin of the 
p-meson track. The projection on the surface of the emulsion of this 
solid angle is the sector already referred to. 

(2) The pair must lie within a small calculable solid angle with OE as 
axis since the electrons are projected in the forward direction, 

(3) It must be a pair of electrons of energy 30 MeV. to within the 
experimental error of measurement. oi, 

36 pairs were examined in respect of these criteria. Al, save one, 
were rejected because of violation of (1) and (2). One pair remains which 
satisfied (1) and (2) but did not satisfy (3). This event will be fully 
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described in the appendix to this paper where it will be shown that it is: 
most unlikely to be a pair of electrons. It has not been found possible- 
to furnish a precise identification of its nature. 

It is extremely improbable that any “associated” pairs escaped - 
detection. Firstly, such pairs, if produced, would be more favourably 
situated for observation than the background pairs. Their dip will be 
small on account of the exclusion of steeply dipping u-mesons. Further- 
more, even if their points of origin were undetected in one field, their 
direction would be such that they would be readily observed as converging 
pairs of tracks in the succeeding two or three fields. Secondly, as has. 
already been stated, the superficial density of background pairs is in 
satisfactory agreement with that found by other workers using plates. 
chosen from the same series. 

There remains to be considered a statistical possibility namely that the 
neutral particles are indeed photons, but that they failed to be materialized. 
in a scan-distance in which 5-6 pairs were to be expected. Each scan- 
length y; microns corresponding to the jth 7» event, constitutes a 
“trial? in which the probability of producing a pair is 1-52 x 10~*y,. 
The, probability of failure in this trial is (l1—1-5 x 10~*y,) and the proba- 

n 
bility of failure in all n trials Pp= I (l—1:5x10™y,). A sufficient 


approximation to P, will be Poisson’s Function of Rare Statistical 

Frequency applicable to the case where all trials have equal probability. 
. This yields Pp= e“=«-*77=4 10-3. Rejecting this supposition on the 
» grounds of its extreme improbability, we may conclude with reasonable 
confidence that the neutral particle, whatever be its nature, is not a 
photon. 

§ 4. INTERPRETATION OF RESULTS. 
(a) Spin of the 7 and u-mesons. 

Further support for the reasonableness of this conclusion appears from 
considerations of spin conservation. Experiment has shown that the 
p-meson decays into an electron and more than one neutral particle. 
(Leighton et al. 1949, Davies et al. 1949). Evidence based on the average 
energy of the electron spectrum favours the emission of two neutral 
particles, presumably identical. The results of these workers would 
appear to be confirmed by a more extended series of measurements on 
which the writer is at present engaged. Counter experiments (Hincks and 
Pontecorvo, 1949, 1950 ; Piccioni, 1949) appear to exclude the possibility 
that these are photons. Taken in conjunction with the results of the 
present work, these findings indicate that the decay chain z-p-e takes 
place as follows : 

7 > [L-yy 

> e+2yp. 
Whether the spin of rv, is integral of half-integral, it follows that that of the 
p-meson is half-integral. This result is consistent with evidence based 
on the distribution curve of burst frequency (Christy and Kusaka, 1941). 
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If v, has spin integral or zero, it follows that the -meson has half-integral 
spin. This is in contradiction to the role of the 7-meson in the theory of 
nuclear force fields which envisages interactions of the type p>ntnt ; 
N—>p-+-m, an assumption supported by experimental evidence for the 
creation of 7-mesons by nucleon-nucleon collisions*. We infer that v, 18 
a fermion (a particle of spin 4), a conclusion not contradicted by the 
results of the present experiment. 

Similarly, we may conclude that if the competing decay process 
7 > €+v; exists, it is likely to lead to the production of monoergic electrons 
unless an odd number >1 of neutral particles is emitted. The theoretical 
cross-section for this process is of the order 1/100 that for the chain 
m-y-e (Ruderman and Finkelstein 1949). 


(b) Haistence of the Neutrino. 


It appears that the characteristics of the decay of 7-mesons may have 
a close bearing on this problem. It would be tempting to identify, 
v;, vy, and v, with the neutrino of B-decay. In common with this particle, 
it would seem that they are fermions of negligible rest-mass. However, 
the present state of theory is such that perhaps it would be unwise to press 
the analogy between z-decay and the £-decay of radioactive elements. 

It is a necessary consequence of the assumed properties of the 
“neutrino ”’, that such particles have a very weak interaction with matter. 
Attempts have therefore been made to prove the existence of the particles 
indirectly. It has commenly been considered that satisfactory evidence 
would consist in a demonstration that the conservation of momentum, 
as well as that of energy, requires the assumption of the emission of a 
neutral particle in the £-decay of a nucleus, or in the process of K-capture. 

Several such experiments have been described involving the f-ray 
spectra of C1°* (Crane and Halpern 1938, 1939) Kr®* (Jacobsen and 
Kofoed-Hansen, 1945) and on the formation of Li’ by K-capture from 
Be’ (Allen 1942), the magnitude of the recoil momentum has been shown 
to be consistent with the neutrino hypothesis. These experiments, 
though brilliantly conceived, are attended by great difficulties, because 
the disparity between the masses of the electron and the nucleus has the 
consequence that the energy of recoil of the latter is very small. 

The present experiment has points of similarity with that of Allen in 
which the recoil momentum of the Li? was shown to be due entirely to 
the emission of the neutrino. The recoil energy was only about 50 eV. 
and was measured electrostatically. The present method is similar 
in principle, but the p-recoil has energy 4-1 MeV., is visible in the emulsion 
and its momentum is measured by use of the well established range-energy 
relationship. If subsequent developments should establish the fact 
that of the neutral particle is identical with the B-decay neutrino, the 


* Note added in proof—Additional direct experimental evidence for the con- 
clusion that the 7-meson is not a fermion comes from the recent work of 
Panofsky, Aanodt and York (Phys. Rev., 78, 825, 1950) on the y-radiation 
associated with the absorption of 7—mesons in hydrogen under high pressure. 
This would appear to demonstrate the occurrence of the reaction am +p—>n+y. 
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present experiment could be regarded as supplementing and confirming 
the B-ray experiments. In view, however, of the anomalous nature of the 
particle apparently associated with a case of 7-y-e decay, already referred 
to, and described in detail in the Appendix, it would seem at present 
unwarranted to conclude that, the neutral particles emitted in meson 
decay are to be identified with the “neutrino” postulated by Fermi and 
by Pauli in order to explain the observed phenomena of B-decay. 


APPENDIX. 


In the course of the search for electron pairs which could be attributed 
to y-rays liberated by the decay of 7-mesons, an event was observed of 
which photo-micrographs are reproduced in Pl. XXIX. Although no 
definite conclusion can be drawn from it, it displays a number of remarkable 
features, and it was therefore found desirable to describe it in detail : 


Fig. 2. 
? 


4,=-1-4° 


4-3° 26° 


1500 » from the point of arrest, O, of a 7-particle in the emulsion, and 
in the line of motion of the .-meson produced by the decay of the z-particle, 
and which subsequently decayed itself into an electron of 36 MeV., a pair 
of particles p, and p, appear to originate at a point P. OP is indis- 
tinguishable from the initial line of motion of the u-meson. Further, 
OP is inclined to the plane containing the initial directions of motion of the 
two particles p, and p, at an angle less than 1°. These relationships, the 
angles of dip (A), positive if towards glass, of tracks p, and p, and the 
lengths (R) of the tracks are shown in fig. 2, ; 
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- Although the ranges of the two particles in the emulsion are relatively 
short, it has been found possible to make a determination of the scattering 
parameter o (Fowler 1950), for the tracks are heavily scattered and short 
“ cell-lengths ” of 20 1 can therefore be employed. The values of g thus 
obtained are :—for p,, %1o9,—1-63-L0-25 ; and for p,, 100. 1-60-35. 
The grain-density in the two tracks are 16-3-43-0 and 20-4 grains per 
504, respectively, values which may be compared with those due to 
particles with minimum ionization recorded in the same plates ; namely, 
11-5 grains per 504. The uncertainties quoted are standard deviations. 
_ Ifthe tracks are due to particles with the elementary electronic charge, 
the observed values of the scattering parameter correspond to particles 
of kinetic energy ~ 10 MeV., and in view of the considerable limits of 
error due to statistical fluctuations, it would therefore be reasonable, 
from the point of view of the kinetic energy of the particles, to attribute 
them to electrons produced by a 30 MeV. y-ray. Other features of the 
observations, however, make it very difficult to accept this explanation. 

Firstly, the grain-density in the tracks is exceptionally high if they 
are due to particles moving at relativistic velocities, as would be the case, 
over most of the trajectory, if the particles were electrons. Here again, 
the statistical fluctuations, due to the relatively restricted lengths of the 
trajectories is rather large, and the evidence is not decisive. Observations, 
however, on more than 1000 tracks, of length +3000 p, of which the grain 
density has the minimum value, have been made in this laboratory in 
connection with other investigations. These observations were analyzed 
to determine the frequency with which the grain-density in the tracks 
of particles of high energy, measured over a length of 500, gave an 
apparent grain-density of 16 or greater. Inno case was a value observed 
as high as the mean for the two tracks under discussion. 

A second feature of the event is that the two trajectories both appear 
to terminate in the emulsion without any indication of an increase in 
the specific ionization due to a gradual reduction in velocity. Such an 
effect might be explained as due to the annihilation of a positive electron, 
and to a concomitant large angle deviation in the track of an electron. 
Since both of the particles cannot be positive electrons, the apparent 
endings of the tracks have been scrutinized. They occur in the middle 
of the emulsion where there are no abnormalities of any obvious kind. 

A third feature of the observations is the angle between the pair, 
viz. 9:5°. This angle is very large if the pair is attributed to the 
materialization of a 30 MeV. y-ray, but not impossibly so. Because of 
the occurrence of these three unusual features in the same event, it is 
reasonable to reject the view that it is due to a pair of electrons. In this 
case, we have to account for a phenomenon in the line of motion of a 
neutral particle emitted in the decay of a 7-meson with exceptional 
features not previously observed. 

The values of grain-density and scattering constant % correspond to 
masses of the order 30m, for both tracks, but, as is well known, 
statistical fluctuations may lead to serious error in an estimate of mass 
based on observations on a single relatively short track. 
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If we accept the evidence of the grain-density and scattering of the 
track as not being seriously in error, it appears to be very difficult to 
account for the phenomena in terms of charged particles of which the 
existence is established. It is also difficult to find a consistent 
explanation of the observations in terms of known particles which have, 
by chance, passed through the emulsion during development. Two 
protons, for example, originating in a nuclear explosion during develop- 
ment would produce thin tracks but they should show evidence of an 
increase in grain-density towards the end of their range, a feature which 
is not observed. It appears best, however, to defer an attempt to identify 
the event until a second obser'vation of the same kind allows us to exclude 
the possibility that it is to be attributed to a rare but trivial occurrence. 
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LXXVI. CORRESPONDENCE 


An Estimate for the Position of the Lowest Dipole-level of a Nucleus. 
By B. F. Touscuex, 
eo men of Natural Philosophy, The University, Glasgow *. 
[Received June 20, 1950.] 
* Communicated by the Author. 


THE completeness property of the set of nuclear eigenfunctions can be 
used to obtain an inequality for the position of the lowest dipole level 
of a nucleus, i.e. for the lowest level that can be reached from the 
groundstate by a dipole transition. 


FEENBERG (1935) has applied sum rules to obtain a lower limit for the 
energy of a hypothetical excited state of He’. His method can be adapted 
to apply to any nuclear system, though the actual evaluation of the 
limiting value for the excitation energy can only be carried out on the 
basis of a strongly simplified nuclear model. By considering two extreme 
cases it will be shown that the limiting value of the excitation energy 
does not depend critically on the model used for its derivation. In this 
paper the method will be outlined for dipole transitions though, of course, 
a similar procedure may be applied to quadripole and higher transitions. 
The estimates obtained for quadripole transitions give a much higher 
limiting energy than in the dipole-case (140A~2 MeV. as compared with 
50A~2/8 MeV. for a dipole) and are therefore probably less useful. 
The basic idea is to evaluate the sums 


| eer Pert we clin. ast. cated) 
and e 
Heer ests te 5 52 | ai (2) 

b 


by applying the completeness of the set of nuclear eigenfunctions, 
%q, is the electric dipole moment (after the separation of the motion of 
the centre of gravity) of the transition a > 6b between the states a and b 
’ with energies E, and E, of a nucleus. If a=O denotes the groundstate 
and b=1 the first Reo state with x), 40, then all pee E,—E, are positive 


and one has 
E,—E, laud. . . . . . . . . (3) 


In a one particle picture in which the nucleus consists of independent 
neutrons and protons the sums (1) and (2) can be readily evaluated : 


R2Z h? 
xX =— 5 2(1+I) and H=>5 Z4(1+D), 


where [=1—2Z/A is the isotopic ratio, M is the nucleon mass and R the 
radius of the nucleus. The factor 1/5 corresponds to the assumption 
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that the nucleons are uniformly distributed over a sphere of radius R. 
The factor 4(1-+1) takes account of the motion of the centre of gravity : 
it vanishes for a system consisting of particles of one and the same specific 
charge. Taking RA~'?=1-5 x 10-% em. one obtains the estimate 


Eee MeV, eee 


a 
for the position of the first dipole level of a nucleus. 

It is well known that the independent particle model gives a poor 
description of the sum-over-cross-sections for nuclear photo effects and 
recent calculations by Lévinger and Bethe (1950) on the yield of photo 
neutrons have shown a strong dependence of the results on the assumed 
exchange character of nuclear forces. This can readily be understood 
if one considers the tendency of group formation inherent in a force-picture 
with saturation properties. This group-formation will result in strong 
correlations between the coordinates of different particles (in particular 
between two neutrons and two protons forming an intranuclear «-particle), 
so that considerable changes may be expected in the values of H and X. 

To see how the estimate (4) is affected by different assumptions about 
nuclear forces one can consider an independent particle model in which 
the nucleus is assumed to consist of «Z free protons, }Z(1—«) «-particles 
and A—Z—Z(l—«) free neutrons. The parameter « can then be 
interpreted as the degree of dissociation of the «-particles in the nucleus. 
A model of this type with a reasonable value of « (Wergelands (1941) 
results on Be® suggest a value «~0-25) should be expected to give a fair 
rendering of some of the effects of exchange forces: the formation of 
a-particles as relatively stable sub-units and the resultant correlations 
between the coordinates of the nucleons. In this model the sums (1) 
and (2) can again be evaluated in an elementary manner : 


xats {I+-«+3I?(1—«)} 
RZ, » Syery eee! 
= Gq t+} 


In the case of complete dissociation of the «-particles this gives, of course, 
the result already obtained for the free neutron-proton model. The 
case =O gives X=—O and H=O if A=2Z. In this case all the particles 
must have the same specific charge so that both the cross-section and the 
average energy-loss of dipole transitions must vanish. From (5) one 
obtains the refined estimate for the position of the lowest dipole level : 

5h? 1 ; 
SMRA I A2319(1 2a)(lae |) a 


and it is immediately obvious that for moderate isotopic ratios this 
expression is nearly independent of «. Even for very heavy nuclei the 
difference between the values of (6) for «=O and «=1 respectively is not 


E,—E,y<= 
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more than 40 per cent. This small variation seems to indicate that the 
estimate of the position of the first dipole level is better than the theory 
from which it has been derived and suggests to try to determine a similar 
limit for quadripole levels by means of the free neutron-proton model. 
Again an elementary calculation similar to that leading to the inequality (4) 
gives 

(F,—E oq <n MeV. PAs Mane een 7) 
This limit is not very useful for light nucléi since there it corresponds 
to an energy far in excess of the “ionization-energy ”’ of a neutron. 
For A>70 it should fall into the discrete part (<8 MeV.) of the nuclear 
spectrum so that one concludes that every nucleus with A>70 should 
have at least one quadripole level in the discrete spectrum. The 
corresponding consideration for dipole levels gives a critical atomic 
number of A~16 above which every nucleus should have at least one 
discrete dipole level. Feenberg’s calculations seem to indicate that the 
actual level usually occurs very near to its theoretical limit. O18 with a 
dipole level at 7 MeV. and a limiting energy of approximately 8 MeV. 
is an example for this behaviour. 
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Examples on the Calculation of Energy States by the Uncertainty Relations. 


By Hsry P. Sou, 
National Chekiang University *. 


[Received June 14, 1950.] 


AuruoucH the Uncertainty Principle forms one of the most fundamental 
concepts in quantum mechanics, the possibility of making quantitative 
predictions based on this principle is little explored. In this letter we 
shall give some examples illustrating the possibilities in this direction. 


I. Hydrogen atom. 
The energy of the system satisfies the inequality 
1 ze 
— 2a — 
ee 2m Bs ean 
since the contribution of angular motion is not included in the right side 
of the above expression. The minimum energy is given by the “spread” 
(mean square deviation) of P, and that Oley Wwe. : 
1 5 ere" 
APs— Tr’ 


* Communicated by Professor M. Born, F.R.S. 


E 


min Im 
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where AP, is the minimum measurable momentum and Ar is the minimum 
measurable distance of the electron from any other particle, here the 
nucleus, The minimum of P, is now postulated to be connected with that 
of r by the relation JP,4r>h*. Hence 


h2 1 ze? 


Emin Sm Ze? Ar 


The last expression can be minimized by finding a suitable value of Ar, 
- which may be varied by changing, for instance, the wavelength of the 
light used in the measurement. Using this minimum value, Ar=h?/ze"m, 
we get finally, H>—z%e4m/2h?. 


II. Relativistic hydrogen atom. 


For the relativistic hydrogen atom one has 
24 2p2]s 2e 
BH >[mpct+c?Po]t— aE —m,C?. 


The argument in the case of this problem is the same as in the last 
example, 2. e. 


ze h 
E>[mct+c? AP2}}— UE —m¢?> | mictee oar | — — —m,C?. 


Varying Ar and substituting in this formula, we obtain 


E>m,c[(1—«?z7)?—1], — 


Ill. Hydrogen atom inchuding the motion of the nucleus. 
In this case 


ia sore 
E> oy et a, arene 

The minimum energy is given by the spread of Pp, P.and r,,. But the 

minimum of 7,,, the distance of the electron from the nucleus, cannot be 

smaller than either of the minimum measurable distance of the electron 

from any other particle, 47, or of the minimum measurable distance of 

the nucleus from any other particle, AR. Hence if 4r>AR, then 


(eo Ah he ol ze" 
b> 2 See 
‘ 2M AR?’ Im dP Ary 


* In the case of two dimensional motion this e i i 
xpression can. be rigorousl 
derived from the wave property of the electron. We shall, Kowoverr fie it 


as a postulate. In dealing with more complicated atoms we shall 
AP, dr >nh for “ higher quantum states - Se an 


—————— 


<< 
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Now the first term is independent of the other terms and has its minimum 


value when AR has the greatest value, which by our assumption is at 
most equal to dr. Hence 


The same lower limit is obtained if A4r< AR is assumed. 
In the three examples given above the lower limit of energy obtained 


from uncertainty relation agrees in every case with that obtained by 
other methods. 


IV. Helium-like atoms. 
For helium-like atoms 
E> ss (P2 + P2)—2e? (= =f ~) 4: 7 :; 
The minimum of the first four terms is given by the spread of P.., P,., 
r,,andr,. The minimum of the last term is given by the largest value of 
*12 Which must be less than 7,+7,. The last limit requires that the three 


_ particles should lie on a straight line which is hardly realizable. However, 
the inequality 


E h? 1 1 i dee a 1 ra e 
> om \ Are > Te dr, * dr,) * dr,+ Ar, 
is certainly true. It sets a lower limit to the actual minimum. Mini- 
mizing the above expression by varying Ar, and Ar, we get 
E>—2(z—4)?Rh. 


In the following Table the limits thus calculated are given together with 
the experimental minima for the first few atoms. 


jal= He Li Be B C 
Den O58: —45:12  —66-12 
E (cal.) ia eee G19 1512 28-12 45-12 
E (exp.) —1-05 Ba 8790 14.664, 2756  —44.06% —64-8 


The discrepancy between the theoretical and the experimental values is 
to be expected from the crudeness of our estimate. The significant 
point is that all the calculated values lie lower than the empirical ones. 
The agreement can be much improved by considering the uncertainty in 
the electronic distributions. Hither by averaging over all directions, 
remembering that r,, must be greater than either 4, or 4,,, or by 
assuming correlated normal distributions for the electrons, we obtain, 
practically in all cases, results in reasonably good agreement within 
1/2 to 1/5 per cent, with the experimental values. 


SER. 7, VOL. 4I, NO. 319.—AUGUST 1950 3N 


854 Correspondence 


Other complicated systems, such as Li-like and Be-like atoms, have 
been treated in the same way with reasonably good results, by assuring 
Pauli’s principle and using the relation 4P,4dr>2h for the third and 
fourth particles (see footnote). It is, however, needless to go into details 
until the physical ideas presented here, for instance, the argument in the 
third example, can be formulated on a mathematical basis capable of 
extension to many-body problems. 


I wish to thank Mr. W. T. Soh for some of the earlier calculations on the 
Berylium system and Dr, T. M. Hu for much of the later calculations. 


December 18, 1949. 
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LXXVII. Notices of New Books and Periodicals received. 


Radiative Transfer. By S. CHANDRASEKHAR. (Oxford : The Clarendon Press.) 
[Pp. 393.] Price 35s. 


Tuts book presents a systematic treatment of transfer in a plane parallel 
scattering medium, a subject which has been investigated principally by 
astrophysicists, but acquires new applications in the scattering of neutrons: 
indeed, other applications will assuredly be found in chemistry and applied 
physics for the quantitative results now available. The main line of attack 
consists in using certain principles of invariance (as that the diffuse scattering 
from a semi-infinite plane parallel atmosphere should be invariant to the 
addition of a layer of arbitrary optical thickness) to formulate integral equa- 
tions : a method introduced by Ambarzumian in 1943 and greatly developed 
by Chandrasekhar himself. Astonishingly, the resulting non-linear inhomo- 
geneous integral equations are then solved. Those who would recoil from the 
more difficult problems which are treated with success here, or who take no 
interest in radiative transfer, may nevertheless study with profit the application 
of Gaussian quadrature (and its generalization) to the approximate (and some- 
times exact) solution of certain simpler integro-differential equations. 

It is presumed,that this will remain for long the standard text on its subject. 
But even should its techniques be superseded, the easy and illuminating style, 
and the continuous courtesy to the reader—in keeping him informed of where 
and how an argument is leading ; in converting the plain “ list of references ”’ 
at the end of each chapter into concise and instructive ‘‘ bibliographical notes ” ; 
in providing an index of definitions in addition to that of names and subjects— 
will allow it to stand as a model of the way in which a text in mathematical 
physics should be written. F.C. F. 


Absorption Spectrophotometry. By G. F. Lorutan, M.A., F. Inst.P. (A Hilger 

Publication.) Price 26s. net. ° 
THis book caters for the reader with a scientific background, commencing 
research in absorption spectrophotometry or requiring an informed survey of 
the practice and possibilities of the subject presented within reasonable compass. 
It succeeds admirably in these aims, and at the same time will be valued by 
experienced workers in the field as a balanced account of up-to-date practice. 

The book is divided into three parts :—The first is introductory and assesses 
conditions affecting accuracy of measurement, and here, as throughout the book, 
the mathematics is an absolute minimum. : 

The second part, and perhaps the most valuable, presents a most interesting 
and varied collection of the types of problem arising in modern practice with 
appropriate methods of attack. Few readers will be familiar with the whole 
field surveyed, but specialists and non-specialists alike will profit by this 
comprehensiveness. os 

The third part discusses apparatus and techniques. The descriptions of 
apparatus are most effective, but those of actual techniques somewhat com- 
pressed. 

The book is made up from good paper, is well illustrated with 70 figures and. 
photographs, and also lists some 300 references, well indexed, and half of 
them within the 1944-48 period. This ample documentation mitigates those 
accounts which are necessarily brief in order to keep the book to its present 
aii reviewer finds a score of minor errors; a more careful layout of the 
ion of the figures with reference to the text and a greater consistency in 


osit 
k However, this is a good work and pleasant to read. 


lettering would not be amiss. 
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Cosmic Rays. By L. JAnossy. 2nd Edition. [Pp. 454, Plates 10.] (Oxford : 
Clarendon Press, 1950.) Price 40s. net. 


Tue first edition of this work was published only two years ago and was reviewed 
in this Journal (39, 826, 1948) as authoritatively written by an investigator who 
has himself made numerous contributions to the subject. The need for a new 
edition at so early a date brings with it difficulties for the author, because 
important contributions to’ the subject are continuing to appear in rapid 
succession through the application of the photographic emulsion technique. 
It was probably inevitable at this stage that the attempt to bring the work up 
to date had to be made by what one may term tinkering and repair rather 
than by re-writing. But when a third edition is planned, changes in construction 
and in emphasis will undoubtedly be imperative. A.M.T. 


History of Physics. By Max Von Lavz, translated by R. E. Oesper. 
[Pp. 150.] (Academic Press, New York.) Price $2.30. 


In a foreword to his book, Professor von Laue says that the first draft was a 
war product written in 1943. ‘ Precisely because of this circumstance the 
thought of the culture that is the common property of all nations, and which 
was then being so despicably mistreated, was ever present throughout the period 
of composition”’. A very readable book by an eminent physicist. 


Manual of Organic Chemistry for Advanced Students, Vol. I. By G. Dyson. 
[Pp. 984.] (Longmans Green.) Price 63s. 


Tuts is the first part of a three volume treatise and deals with the compounds 
of carbon, hydrogen, oxygen and the halogens. Its scope may be indicated 
by saying that it is more comprehensive than a text-book, without being a 
collection of specialized monographs. Compiled by one who has made very 
important contributions in the field of notation in organic chemistry it may be 
looked upon as authoritative. Readers of this journal should find the volume 
useful as a work of reference. 


BOOKS RECEIVED 


The Association of Research Institutes of Tohuku University Japan have 
started to issue a series of Science Reports consisting of collected papers on 
work carried out by them in Physics Chemistry and Metallurgy. Volume I. 
Numbers 1 and 2 (May and August 1949) have been received. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


